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Abstract 

Passive  multiple-input  multiple-output  (MIMO)  radar  is  a  sensor  network  comprised 
of  multiple  distributed  receivers  that  detects  and  localizes  targets  using  the  emis¬ 
sions  from  multiple  non-cooperative  radio  frequency  transmitters.  This  dissertation 
advances  the  theory  of  centralized  passive  MIMO  radar  (PMR)  detection  by  propos¬ 
ing  two  novel  generalized  likelihood  ratio  test  (GLRT)  detectors.  The  first  addresses 
detection  in  PMR  networks  without  direct-path  signals.  The  second  addresses  de¬ 
tection  in  PMR  networks  with  direct-path  signals.  The  probability  distributions  of 
both  test  statistics  are  investigated  using  recent  results  from  random  matrix  theory. 
Equivalence  is  established  between  PMR  networks  without  direct-path  signals  and 
passive  source  localization  (PSL)  networks.  Comparison  of  both  detectors  with  a 
centralized  GLRT  for  active  MIMO  radar  (AMR)  detection  reveals  that  PMR  may 
be  interpreted  as  the  link  between  AMR  and  PSL  sensor  networks.  In  particular,  un¬ 
der  high  direct-path-to-noise  ratio  (DNR)  conditions,  PMR  sensitivity  and  ambiguity 
approaches  that  of  AMR.  Linder  low-DNR  conditions,  PMR  sensitivity  and  ambigu¬ 
ity  approaches  that  of  PSL.  At  intermediate  DNRs,  PMR  sensitivity  and  ambiguity 
smoothly  varies  between  that  of  AMR  and  PSL.  In  this  way,  PMR  unifies  PSL  and 
AMR  within  a  common  theoretical  framework.  This  result  provides  insight  into  the 
fundamental  natures  of  active  and  passive  distributed  sensing. 
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PASSIVE  MIMO  RADAR  DETECTION 


I.  Introduction 

This  dissertation  advances  the  theory  of  passive  multiple-input  multiple-output 
(MIMO)  radar  detection.  It  also  unifies  passive  MIMO  radar  (PMR),  active  MIMO 
radar  (AMR),  and  passive  source  localization  (PSL)  sensor  networks  within  a  common 
theoretical  framework.  This  chapter  introduces  PMR  in  Sec.  1.1,  AMR  and  PSL  in 
Sec.  1.2,  and  discusses  the  relevance  of  passive  radar  to  the  United  States  Air  Force 
in  Sec.  1.3.  Current  approaches  to  passive  radar  detection  and  their  shortcomings 
are  discussed  in  Sec.  1.4,  and  the  methodology  used  to  address  these  shortcomings  in 
Sec.  1.5.  The  contributions  of  this  research  are  summarized  in  Sec.  1.6,  followed  by 
an  outline  of  the  remainder  of  this  dissertation  in  Sec.  1.7. 

1.1  Passive  Radar 

Passive  radar  is  a  type  of  radar  that  exploits  non-cooperative  radio  frequency 
(RF)  transmissions  to  detect,  localize,  track  and/or  image  objects  of  interest.  Non- 
cooperative  transmitters ,  also  known  as  illuminators  of  opportunity,  operate  sepa¬ 
rate  from  and  independent  of  the  passive  radar  system.  Example  non-cooperative 
transmitters  include  commercial  broadcast  transmitters,  cellular  phone  base  stations, 
navigation  satellites,  and  air  surveillance  radars  [88].  The  use  of  non-cooperative 
transmitters  distinguishes  passive  radar  from  active  radar.  Consequently,  a  passive 
radar  system  contains  only  receivers,  and  it  must  select  from  the  illuminators  within 
its  operating  environment.  This  enables  covert  operation,  yet  also  presents  challenges 
that  constrain  performance  in  terms  of  sensitivity,  resolution,  and  ambiguity  [31]. 
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Passive  radar  dates  to  the  earliest  days  of  radar.  In  1935,  Sir  Robert  Watson- Watt 
detected  an  aircraft  using  shortwave  broadcasts  from  the  Daventry  British  Broadcast¬ 
ing  Corporation  Empire  transmitter  [57].  During  World  War  II,  the  German  system 
Klein  Heidelberg  exploited  transmissions  from  the  British  Chain  Home  air  surveillance 
radar  to  passively  detect  Allied  bombers  [32],  In  the  1980s  and  1990s,  researchers 
in  the  United  States  and  the  United  Kingdom  revisited  using  commercial  broadcast 
transmitters  for  aircraft  detection  [34,  45,  46].  Several  demonstration  systems  re¬ 
sulted,  including  Lockheed  Martin’s  Silent  Sentry,  which  exploits  multiple  frequency 
modulated  (FM)  radio  transmitters  to  track  aircraft  [7].  Interest  in  passive  radar  grew 
rapidly  in  the  wake  of  these  investigations.  Fig.  1  depicts  the  number  of  passive  radar 
publications  that  appeared  in  the  English-language  technical  literature  from  1999  to 
2011.  Many  prototype  systems  were  reported  during  this  period.  Three  prominent 
examples  include:  Passive  Radar  Demonstrator  (PaRaDe),  developed  by  the  Insti¬ 
tute  of  Electronic  Systems  at  the  Warsaw  University  of  Technology  [59,  60];  COvert 
RAdar  (CORA),  developed  by  the  German  Fraunhofer  Institute  for  High  Frequency 
Physics  and  Radar  Techniques  (FHR)  [55];  and  Gruppenantenne  fur  militarische  Mo- 
bilfunk  Aufklarung 1  (GAMMA),  developed  by  the  German  Fraunhofer  Institute  for 
Communication,  Information  Processing  and  Ergonomics  (FKIE)  [65,  91]. 

While  most  reported  passive  radar  systems  use  only  a  single  illuminator  type,  a 
growing  number  of  multiband  systems  are  also  under  development  [8,  28,  56,  73]. 
The  term  multiband  refers  to  the  simultaneous  exploitation  of  different  illuminator 
types  across  multiple  frequency  bands.  For  instance,  the  Cassidian  Demonstrator  ex¬ 
ploits  FM  radio  and  digital  audio  broadcast  (DAB)  signals  in  the  very  high  frequency 
(VHF)  band  and  terrestrial  digital  video  broadcast  (DVB-T)  signals  in  the  ultra  high 
frequency  (UHF)  band  [73].  Multiband  operation  increases  frequency  diversity  and 
1  German  for  Array  Antenna  for  Mobile  Military  Reconnaissance 
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Year 

Figure  1.  Passive  Radar  Publications  by  Year  and  Type 


geometric  diversity  in  comparison  to  single-band  operation,  and  thereby  provides  both 
detection  and  localization  advantages.  Both  types  of  diversity  increase  the  number  of 
independent  samples  of  the  target’s  reflectivity  signature,  which  typically  decorrelates 
with  frequency  and  angle,  available  to  the  radar  system  for  processing.  The  result¬ 
ing  detection  improvement  is  termed  diversity  gain  [24],  Geometric  diversity  also 
increases  the  probability  the  target  Doppler  is  unmasked  from  clutter  interference 
in  at  least  one  bistatic  (transmitter-receiver)  channel;  this  is  termed  geometry  gain 
[26].  Finally,  the  additional  measurement  degrees-of-freedom  provided  by  geometric 
diversity  enables  direct  localization  of  targets  in  Cartesian  space  [39,  63]. 

In  addition  to  multiband  operation,  airborne  operation  has  also  received  signif¬ 
icant  interest  in  recent  years  [13,  22,  52],  Airborne  operation  extends  the  range  of 
viable  passive  radar  operating  modes  to  include  synthetic  aperture  radar  imaging  and 
ground  moving  target  indication,  while  at  the  same  time  complicating  the  detection 
problem  by  introducing  Doppler-spread  clutter.  To  assess  the  prospects  of  airborne 
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passive  radar,  the  North  Atlantic  Treaty  Organization  established  an  international 
task  group  in  2012  to  “comprehensively  research  the  current  state  of  the  art  of  air¬ 
borne  passive  radar,  identify  technological  challenges  including  hardware  and  signal 
processing  development,  identify  end  user  requirements,  and  address  operational  use¬ 
fulness”  [1].  While  only  FM  radio  appears  to  have  been  used  in  experimental  airborne 
passive  radar  trials  to  date  [13,  22],  DVB-T  has  also  been  identified  as  a  candidate 
illuminator  for  airborne  operation  [52], 

1.1.1  Network  Topology. 

Although  passive  radar  systems  have  only  receivers,  a  passive  radar  system  and 
the  non-cooperative  transmitters  exploited  by  that  system  are  sometimes  described 
as  a  passive  radar  network.  The  topology  of  such  networks  is  described  as  bistatic , 
multistatic ,  or  multiple-input  multiple- output,  depending  on  the  number  of  receiver 
sites  and  the  number  of  exploited  non-cooperative  transmitters.  These  topologies  are 
illustrated  in  Fig.  2.  Passive  bistatic  radar  (PBR)  consists  of  one  receiver  and  one 
non-cooperative  transmitter  (Fig.  2a).  Passive  multistatic  radar  consists  of  either:  (1) 
one  receiver  that  exploits  multiple  non-cooperative  transmitters  (Fig.  2b);  or  (2)  mul¬ 
tiple  geographically  separated  receivers  that  exploit  one  non-cooperative  transmitter 
(Fig.  2c).  Type  1  is  more  common  than  type  2  clue  to  the  availability  of  transmitters 
and  the  simplicity  of  using  one  receiver  site.  Finally,  passive  MIMO  radar  (PMR) 
consists  of  multiple  geographically  separated  receivers  and  multiple  non-cooperative 
transmitters  (Fig.  2d).  In  comparison,  PMR  is  the  most  general  passive  radar  topol¬ 
ogy,  as  both  bistatic  and  multistatic  systems  may  be  regarded  as  special  cases  of 
PMR.  Consequently,  PMR  is  the  primary  focus  of  this  work. 
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(a)  Passive  Bistatic  Radar 


(c)  Passive  Multistatic  Radar,  Single-Tx 


-  Target-Path 

-  Direct-Path 

- Target- Path 

-  Direct-Path 

(b)  Passive  Multistatic  Radar,  Single-Rx 


(d)  Passive  MIMO  Radar 


Figure  2.  Passive  Radar  Network  Topologies 


1.1.2  Signal  Environment. 

In  the  PMR  signal  environment,  direct-path  (i.e. ,  transmitter-to-receiver)  and 
target-path  (i.e.,  transmitter-to-target-to-receiver)  signals  from  each  transmitter  are 
incident  at  each  receiver.  This  is  illustrated  in  Fig.  2d.  The  transmitters  within  a 
PMR  network  typically  occupy  distinct  frequency  channels.  However,  some  commer¬ 
cial  illumination  types,  including  DVB-T  and  DAB,  may  operate  as  single  frequency 
networks,  in  which  multiple  geographically  separated  transmitters  emit  the  same  sig¬ 
nal  on  the  same  frequency  channel.  Exploitation  of  single  frequency  networks  compli¬ 
cates  the  detection  and  tracking  problem  because  it  introduces  measurement  origin 
uncertainty,  i.e.,  the  passive  radar  receiver  does  not  know  a  priori  which  detections 
go  with  which  transmitters  [79].  Sophisticated  algorithms  are  required  to  resolve  this 
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uncertainty  [9,  21].  This  work  considers  only  frequency-separable  transmissions. 

Many  commercial  transmitters  emit  continuously  or  with  high  duty  cycle,  e.g.,  ra¬ 
dio,  television,  and  cellular  phone  emitters.  Such  transmitters  are  termed  continuous 
wave  (CW)2.  Consequently,  direct-path  and  target-path  signals  are  often  received 
simultaneously.  This  causes  direct-path  self-jamming,  which  can  degrade  detection 
sensitivity  due  to  the  large  power  ratio  between  direct-path  and  target-path  signals 
typical  in  passive  radar  scenarios  [33].  Indeed,  power  ratios  exceeding  100  dB  have 
been  reported  [27].  To  counter  this  obstacle,  a  passive  radar  receiver  must  have 
large  dynamic  range  and  employ  one  or  multiple  direct-path  interference  suppression 
techniques  [33].  Two  common  interference  suppression  techniques  are:  (1)  adaptive 
beamforming  [36,  78,  91],  in  which  nulls  are  placed  in  the  interference  direction(s); 
and  (2)  adaptive  filtering  [15,  19,  35,  62,  69],  in  which  the  interfering  signals  are 
estimated  and  cancelled  via  subtraction.  In  addition,  clutter-path  (i.e.,  transmitter- 
to-clutter-to-receiver)  signals  are  also  typically  incident  at  each  receiver.  Clutter-path 
interference  presents  many  of  the  same  challenges  as  direct-path  interference  and  is 
countered  using  many  of  the  same  techniques. 

1.2  Related  RF  Sensor  Networks 

Although  PMR  is  the  primary  focus  of  this  dissertation,  AMR  and  PSL  sensor  net¬ 
works  are  also  discussed  for  comparison.  This  clarifies  the  unique  properties  of  PAIR, 
and  enable  AA1R,  PMR,  and  PSL  to  be  understood  within  a  common  theoretical 
framework.  AA1R  and  PSL  are  briefly  described  in  the  following  sections. 

2As  opposed  to  pulsed. 
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1.2.1  Active  MIMO  Radar. 


Active  MIMO  radar,  illustrated  in  Fig.  3,  is  defined  as  a  system  that  detects 
and  localizes  targets  using  multiple  distributed  transmit-receive  (Tx-Rx)  nodes  and 
multiple  separable  signals  [41].  Note  that  this  type  of  MIMO  is  also  described  in 
the  literature  by  the  terms  statistical  MIMO  [24],  noncoherent  MIMO  [66],  and  dis¬ 
tributed  MIMO  [85] 3 .  AMR  is  distinguished  from  PMR  by  the  use  of  cooperative 
transmitters  that  transmit  known  signals.  These  signals  are  specially  designed  to 
be  separable  in  one  of  several  possible  domains  (time,  frequency,  Doppler,  or  code) 
[89].  Due  to  this  separability,  the  scattered  returns  from  each  transmitted  signal  can 
be  isolated  at  each  receiver.  The  optimal  Neyman-Pearson  detector  for  AMR  detec¬ 
tion  performs  matched  filtering  per  transmitter-receiver  pair,  noncoherent  integration 
across  pairs,  and  thresholding  [14,  20,  24,  26,  43].  Noncoherent  integration  is  necessi¬ 
tated  by  the  non-isotropic  (non-coherent)  scattering  of  complex  targets  with  respect 
to  the  transmitter-receiver  pairs  in  MIMO  topologies. 

3This  work  does  not  refer  to  the  type  of  MIMO  radar  that  uses  closely-spaced  transmitters. 


Figure  3.  An  Active  MIMO  Radar  Sensor  Network 
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1.2.2  Passive  Source  Localization. 


Passive  source  localization  sensor  networks,  depicted  in  Fig.  4,  detect  and  localize 
targets  by  intercepting  and  processing  target-emitted  signals.  In  comparison  to  PMR, 
a  PSL  system  detects  target-emitted  rather  than  target-scattered  signals,  and  these 
signals  are  often  unknown  a  priori  to  the  PSL  system,  as  in  PMR.  The  detection 
and  localization  processing  in  PSL  sensor  networks  is  typically  based  on  pairwise 
processing  of  the  received  signals  across  receivers.  Conventional  approaches  to  target 
localization  follow  a  two-step  procedure  in  which  intermediate  parameters  are  (1) 
measured  and  then  (2)  combined  to  estimate  the  emitter  position  and  velocity.  Such 
approaches  have  utilized  one  or  multiple  measurement  types,  including:  angle-of- 
arrival  (AOA)  [25,  76];  time  difference  of  arrival  (TDOA)  between  receiver  pairs 
[16,  80];  and  frequency  difference  of  arrival  (FDOA)  between  receiver  pairs,  also 
termed  differential  Doppler  [16,  58].  More  recently,  centralized  approaches  that  detect 
and  localize  emitters  by  operating  directly  on  all  received  signals  have  also  been 
proposed  [5,  12,  83,  86,  87]. 


Figure  4.  A  Passive  Source  Localization  Sensor  Network 


1.3  Relevance  to  United  States  Air  Force 


The  mission  of  the  United  States  Air  Force  is  “to  fly,  fight,  and  win. .  .in  air,  space, 
and  cyberspace”  [2],  Airborne  passive  radar  contributes  to  this  mission  in  several 
ways.  First,  near-peer  adversaries  are  developing  anti-access/area  denial  (A2/AD) 
technologies  to  limit  the  projection  of  US  power  into  strategic  places  [4],  Passive 
radar  counters  these  technologies  by  providing  a  passive  means  to  attain  air  situ¬ 
ational  awareness  in  contested  environments  where  RF  emissions  are  precluded  by 
survivability  considerations.  Second,  airborne  passive  radar  is  less  susceptible  to 
electronic  attack  than  active  radar  in  contested  electromagnetic  environments.  The 
covertness  provided  by  passive  radar  makes  it  difficult  to  jam  a  passive  radar  system, 
and  thereby  provides  an  inherent  type  of  electronic  protection.  Third,  a  congested 
RF  spectrum  aids  rather  than  degrades  passive  radar  by  providing  the  passive  radar 
many  candidate  illuminators  to  exploit.  This  contrasts  with  active  radar,  which  must 
operate  within  spectral  gaps  to  prevent  destructive  interference.  Finally,  the  use  of 
transmitters  of  opportunity  enables  the  exploitation  of  restricted  frequency  bands  not 
otherwise  available  for  radar  use,  e.g.,  various  VHF  and  UHF  commercial  broadcast 
bands.  In  all  of  these  ways,  airborne  passive  radar  contributes  to  the  United  States 
Air  Force  core  functions  of  air  superiority  and  global  precision  strike  [23]. 

1.4  Current  Shortcomings 

Although  the  potential  of  passive  radar  has  been  demonstrated  experimentally,  a 
rigorous  theoretical  foundation  for  passive  radar  detection  does  not  yet  exist.  This  is 
particularly  true  for  PMR  systems,  as  the  majority  of  the  existing  research  has  focused 
on  bistatic  and  multistatic  topologies.  In  the  following,  two  current  approaches  to 
passive  radar  detection  are  described,  and  their  shortcomings  are  discussed.  The 
first,  referred  to  as  detection  with  references,  is  the  approach  used  by  all  reported 
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experimental  passive  radar  systems.  The  second,  referred  to  as  detection  without 
references,  is  a  recently  proposed  technique  for  detection  in  multiple-receiver  passive 
radar  topologies.  Due  to  its  prevalence,  detection  with  references  is  also  referred  to 
as  the  conventional  approach,  while  detection  without  references  is  referred  to  as  the 
alternative  approach.  Both  are  discussed  in  turn. 

1.4.1  PMR  Detection  with  References. 

The  conventional  approach  to  detection  in  passive  radar  is  based  on  matched  fil¬ 
tering.  For  each  bistatic  (transmitter-receiver)  pair,  the  passive  radar  receiver  isolates 
the  direct-path  and  target-path  signals  into  reference  and  surveillance  channels,  re¬ 
spectively.  This  is  accomplished  by  pointing  directional  antennas  at  the  transmitter 
and  the  anticipated  target  region,  respectively  [11,  67],  or  by  digital  beamforming  in 
multichannel  systems  [27,  60,  74,  78,  91],  as  shown  in  Fig.  5.  Adaptive  filtering  is  ap¬ 
plied  to  the  surveillance  channels  to  mitigate  direct-path  and  clutter-path  interference 
[15,  19,  35,  62,  69],  and  equalization  techniques  are  applied  to  the  reference  channels 
to  further  isolate  the  direct-path  signals  [17,  77,  92]4.  The  cross-ambiguity  function 
(CAF)  between  the  reference  and  surveillance  channels  is  then  computed,  and  tar¬ 
get  detection  is  declared  for  CAF  cells  that  exceed  a  threshold  value.  The  bistatic 
range,  bistatic  Doppler,  and  angle-of-arrival  measurements  from  all  such  detections 
are  then  fused  to  localize  and  track  targets  in  Cartesian  space  [21,  50,  61,  63,  73]. 
Note  that  this  approach  is  decentralized  in  that  detection  is  performed  separately  by 
each  bistatic  pair,  and  the  resulting  detections  are  fused  in  subsequent  processing. 

Although  this  approach  often  works  in  practice,  it  is  ad  hoc  because  it  simply 
mimics  the  matched  filtering  operation  used  in  active  radar.  Specifically,  calculation 
of  the  reference-surveillance  CAF  approximates  the  matched  filtering  operation  in  ac- 

4Demodulation-remodulation  of  the  reference  channels  is  also  often  performed  for  digital  trans¬ 
missions  [6,  18,  54,  68,  75,  92]. 


10 


Figure  5.  Reference  and  Surveillance  Channel  Formation 


tive  radar,  where  the  reference  channel  provides  an  estimate  of  the  a  priori  unknown 
transmit  signal,  ffowever,  matched  filtering  is  only  optimal  in  the  Neyman-Pearson 
sense  when  the  transmit  signal  is  known  exactly,  as  in  active  radar.  In  high  direct- 
path-to-noise  ratio  (DNR)  scenarios,  the  reference  channel  provides  a  high-quality 
estimate  of  the  transmitted  signal,  and  close  to  optimal  (matched)  detection  sensi¬ 
tivity  is  achieved.  This  is  the  case  in  many  scenarios  of  interest.  However,  if  the 
direct-path  signal  has  low  DNR,  then  the  reference  channel  is  mismatched  to  the 
originally  transmitted  signal,  and  detection  sensitivity  is  degraded  in  proportion  to 
this  mismatch.  Although  such  mismatch  is  mitigated  to  some  degree  by  reference 
channel  equalization,  it  cannot  be  removed  completely5.  Therefore,  the  conventional 
approach  is  best  suited  for  high-DNR  scenarios. 


1.4.2  PMR  Detection  without  References. 

An  alternative  detection  approach  has  recently  been  proposed  that  does  not  uti¬ 
lize  direct-path  signals  [84,  10].  In  this  approach,  target-path  signals  from  multiple 
widely-separated  receivers  are  isolated  and  cross-correlated  after  appropriate  delay- 
Doppler  compensation.  This  approach  applies  only  to  multiple-receiver  multistatic 
and  MIMO  topologies,  and  it  resembles  techniques  used  for  source  detection  in  PSL 

5Even  after  perfect  demodulation-remodulation  of  digital  signals,  mismatch  still  results  be¬ 
tween  the  transmitted  and  reconstructed  signals  due  to  transmitter-specific  effects  such  as  in¬ 
phase/quadrature  channel  mismatch  and  carrier  drift  [75]. 
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sensor  networks  [83,  40].  It  is  also  centralized  rather  than  decentralized  because  it 
jointly  processes  all  measured  target-path  signals  to  detect  and  localize  targets. 

Since  this  approach  does  not  process  direct-path  signals,  it  is  not  degraded  in  low- 
DNR  scenarios,  unlike  the  conventional  approach.  Rather,  low-DNR  scenarios  aid 
this  approach  by  minimizing  direct-path  interference,  which  simplifies  isolation  of  the 
target-path  signals.  Such  scenarios  are  rare  in  practice  [30],  but  could  result  if  there 
is  physical  blockage  of  the  line-of-sight  paths  between  transmitters  and  receivers  (e.g., 
the  Manastash  Ridge  Radar  [88,  p.  121]),  or  if  the  illuminators  utilize  highly  directive 
transmission  and/or  null  the  receiver  directions.  The  resulting  signal  environment  is 
illustrated  in  Fig.  6.  However,  by  ignoring  direct-path  signals,  this  approach  ignores 
a  potentially  useful  source  of  information  about  the  unknown  transmit  signals.  This 
is  particularly  true  in  high-DNR  scenarios,  in  which  the  direct-path  signal  provides 
a  high-quality  reference  that  can  be  used  for  (noisy)  matched  filtering,  as  in  the 
conventional  approach.  Thus,  this  alternative  reference-less  approach  is  better  suited 
for  low-DNR  scenarios  than  high-DNR  scenarios. 


Figure  6.  Passive  MIMO  Radar  without  Direct-Path  References 
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1.5  Research  Hypothesis  and  Methodology 


Both  of  the  detection  approaches  discussed  in  Sec.  1.4  represent  partial  rather  than 
general  solutions  to  the  passive  radar  detection  problem.  This  is  because  both  ap¬ 
proaches  are  based  on  incomplete  views  of  the  passive  radar  signal  environment.  The 
conventional  approach  assumes  reception  of  distortionless  direct-path  signals,  while 
the  alternative  approach  ignores  direct-path  signals  completely.  In  general,  direct- 
path  signals  are  present  but  distorted.  Consequently,  neither  approach  exploits  the 
passive  radar  signal  environment  to  full  advantage.  This  motivates  the  hypothesis  of 
this  research  program:  formulating  the  passive  radar  detection  problem  in  terms  of 
the  actual  received  signal,  which  minimally  includes  the  target-path  signal,  the  direct- 
path  signal,  and  receiver  noise,  will  result  in  a  detector  that  is  suitable  under  all  DNR 
conditions.  The  performance  of  such  a  detector  is  expected  to  approximate  the  perfor¬ 
mance  of  the  two  current  approaches  under  the  appropriate  DNR  conditions,  i.e. ,  its 
sensitivity  would  approximate  that  of  the  conventional  and  alternative  approaches  un¬ 
der  high-  and  low-DNR  conditions,  respectively.  Furthermore,  its  performance  would 
vary  smoothly  between  these  two  extremes  as  a  function  of  DNR.  Consequently,  it 
would  represent  a  general  solution  to  the  passive  radar  detection  problem. 

Accordingly,  the  methodology  of  this  research  program  is  to:  (a)  formulate  the 
PMR  detection  problem  to  account  for  target-path  signals,  direct-path  signals,  and 
receiver  noise  in  the  detection  hypotheses;  (b)  derive  the  generalized  likelihood  ratio 
test  (GLRT)  for  the  resulting  composite  hypothesis  testing  problem;  (c)  investigate 
the  probability  distributions  of  the  resulting  test  statistic  under  both  hypotheses; 
(d)  illustrate  the  detection  and  ambiguity  properties  of  the  resulting  detector  via 
numerical  simulation;  and  (e)  compare  the  detection  and  ambiguity  performance  of 
this  detector  to  that  of  analogous  GLRT  detectors  for  AMR  and  PSL  sensor  networks. 
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1.6  Research  Contributions 


The  main  contributions  of  this  research  are: 

1.  The  derivation  and  statistical  characterization  of  a  novel  detector  for  the  cen¬ 
tralized  PMR  detection  without  direct-path  references  problem  (Chapter  Iff) 

2.  The  derivation  and  statistical  characterization  of  a  novel  detector  for  the  cen¬ 
tralized  PMR  detection  with  direct-path  references  problem  (Chapter  IV) 

3.  The  introduction  of  a  unified  framework  for  detection  in  active  and  passive 
distributed  RF  sensor  networks  that  encompasses  active  MIMO  radar,  passive 
MIMO  radar,  and  passive  source  localization  sensor  networks  (Chapter  V) 

Each  contribution  is  described  in  detail  in  the  following  sections. 

1.6.1  PMR  Detection  without  References. 

A  novel  GLRT  detector  is  derived  for  centralized  PMR  detection  without  direct- 
path  references,  termed  the  surveillance-surveillance  GLRT  (SS-GLRT)  because  it 
performs  pairwise  processing  of  surveillance  channels  [38].  This  detector  extends 
prior  formulations  to  account  for  multiple  non-cooperative  transmitters  and  multi¬ 
channel  (array)  receivers.  These  extensions  significantly  enhance  probability  of  de¬ 
tection  and  reduce  detection  ambiguity.  Then,  using  recent  results  from  random 
matrix  theory,  the  exact  distribution  of  the  detection  test  statistic  is  identified  under 
both  hypotheses.  This  is  the  first  appearance  of  such  distributions  in  the  context  of 
radar  detection.  These  distributions  show  that  detection  sensitivity  is  only  a  function 
of  the  number  of  received  signal  samples,  the  number  of  transmitters  and  receivers 
in  the  PMR  network,  and  the  average  target-path  input  signal-to- noise  ratio  (SNR). 
Furthermore,  numerical  examples  demonstrate  important  properties  of  the  detector, 
namely,  that  (a)  receivers  and  transmitters  contribute  asymmetrically  to  detection 
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sensitivity  (which  is  uncommon  in  radar  detection),  (b)  integration  gain  grows  non- 
coherently  with  increasing  signal  length,  and  (c)  salient  features  of  its  ambiguity 
function  can  be  explained  in  terms  of  TDOA,  FDOA,  and  AOA. 

1.6.2  PMR  Detection  with  References. 

A  novel  GLRT  detector  is  derived  for  centralized  PMR  detection  with  direct-path 
references,  termed  the  reference-surveillance  GLRT  (RS-GLRT)  because  it  processes 
both  reference  and  surveillance  channels  [37].  This  detector  extends  the  formulation 
of  the  SS-GLRT  to  include  direct-path  signals  with  arbitrary  DNRs.  It  is  shown  that 
calculation  of  the  RS-GLRT  test  statistic,  which  is  expressed  in  terms  of  the  largest 
eigenvalues  of  complex  Wishart  matrices,  entails  operations  that  may  be  interpreted 
as  isolating  the  target-path  and  direct-path  signals  into  surveillance  and  reference 
channels,  respectively.  Although  the  exact  distribution  of  this  test  statistic  appears 
to  not  exist  under  either  hypothesis,  it  is  shown  that  this  statistic  is  approximately 
proportional  under  high-DNR  conditions  to  a  GLRT  statistic  for  AMR  detection, 
termed  the  matched  filter  GLRT  (MF-GLRT),  for  which  probability  distributions  are 
known  under  both  detection  hypotheses.  It  is  also  shown  that:  (a)  its  sensitivity  is 
only  a  function  of  the  number  of  received  signal  samples,  the  number  of  transmitters 
and  receivers  in  the  PMR  network,  the  average  target-path  input  SNR,  and  the 
average  input  DNR;  and  (b)  salient  features  of  its  ambiguity  function  can  be  explained 
in  terms  of  bistatic  range,  bistatic  Doppler,  and  AOA.  The  RS-GLRT  is  the  first 
PMR  detector  that  shows  how  detection  sensitivity  depends  on  target-path  SNR  and 
DNR.  Consequently,  the  RS-GLRT  represents  a  general  solution  to  the  passive  radar 
detection  problem  and  confirms  the  research  hypothesis. 
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1.6.3  Unified  Detection  Framework. 


Finally,  a  novel  theoretical  framework  for  detection  in  active  and  passive  dis¬ 
tributed  RF  sensor  networks  is  proposed.  This  framework,  illustrated  in  Fig.  7, 
encompasses  AMR,  PMR,  and  PSL  sensor  networks.  It  identifies  how  AMR,  PMR, 
and  PSL  are  related  by  simple  transformations  of  their  signal  environments.  These 
relationships  reveal  that  PMR  is  the  key  to  linking  AMR  and  PSL,  which  have  tra¬ 
ditionally  been  regarded  as  distinct,  and  suggest  that  PMR  detection  might  exhibit 
properties  of  detection  in  both  AMR  and  PSL  sensor  networks.  Analysis  of  the  RS- 
GLRT  and  SS-GLRT  detectors  shows  this  is  true,  i.e.,  PMR  detection  performance 
(sensitivity  and  ambiguity)  varies  between  that  of  AMR  and  PSL  as  a  function  of  the 
average  DNR.  In  high-DNR  scenarios,  RS-GLRT  performance  approaches  MF-GLRT 
performance,  which  is  the  GLRT  for  AMR  sensor  networks.  In  low-DNR  scenarios, 
RS-GLRT  performance  approaches  SS-GLRT  performance,  which  is  the  GLRT  detec¬ 
tor  for  both  PMR  without  references  and  PSL  sensor  networks.  This  is  explained  by 
observing  that  PSL  and  AMR  represent  two  extremes  in  terms  of  knowledge  about  the 
transmitted  signals,  i.e.,  the  signals  are  entirely  unknown  in  PSL  and  entirely  known 
in  AMR.  Direct-path  signals  provide  PMR  a  varying  degree  of  knowledge  about  the 
transmit  signals  that  is  quantified  by  the  DNR.  Thus,  PMR  unifies  PSL  and  AMR 
within  a  common  theoretical  framework.  This  result  provides  fundamental  insight 
into  the  natures  of  active  and  passive  distributed  RF  sensing. 


Unknown  Signals  Direct-Path  Unavailable  Transmitting  Target 

Passive  Ml  MO  (ent_  Passive  Source 

without  Ref.  *  Localization 


Known  Signals  Direct-Path  Available  Illuminated  Target 


Low-p 


Figure  7.  A  Unified  Framework  for  Active  and  Passive  Distributed  RF  Sensing 
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1.7  Chapter  Outline 


The  remainder  of  this  dissertation  is  organized  as  follows: 

Chapter  II:  presents  a  signal  model  for  PMR. 

Chapter  III:  considers  the  PMR  detection  without  direct-path  references  problem. 
The  SS-GLRT  detector  is  derived,  and  its  detection  and  ambiguity  performance 
is  analyzed  and  illustrated  via  numerical  simulation.  PMR  detection  without 
direct-path  references  is  shown  to  be  equivalent  to  PSL  detection. 

Chapter  IV:  considers  the  PMR  deteciton  with  direct-path  references  problem.  The 
RS-GLRT  detector  is  derived,  and  its  detection  and  ambiguity  performance  is 
analyzed  and  illustrated  via  numerical  simulation.  PMR  detection  performance 
is  shown  to  vary  between  that  of  detection  in  AMR  and  PSL  sensor  networks 
in  a  manner  that  depends  on  the  DNR. 

Chapter  V:  presents  the  unified  theoretical  framework  for  detection  in  active  and 
passive  distributed  RF  sensor  networks. 

Chapter  VI:  summarizes  the  main  contributions  of  this  research,  and  discusses  how 
it  can  be  extended  in  future  work. 
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II.  Signal  Model 


This  chapter  presents  a  signal  model  for  PMR  that  underlies  Chapters  III-V. 

2.1  Scenario 

Consider  a  scenario  with  Nt  transmitters,  Nr  receivers,  and  one  target.  Fig.  8  de¬ 
picts  the  geometry  of  the  if th  transmitter-receiver  pair,  which  consists  of  the  zth  trans¬ 
mitter  and  jth  receiver.  The  position  and  velocity  of  the  zth  transmitter  are  denoted 
by  d*  and  d*,  i  =  1 ...  Nt,  the  position  and  velocity  of  the  jth  receiver  by  C  and  r-7, 
j  =  1 ...  Nr,  and  the  target  position  and  velocity  by  t  and  t.  The  time  dependence 
of  d*,  dj  r-7,  r-7,  t,  and  t  is  implicit  to  simplify  notation;  in  general,  transmitters, 
receivers,  and  the  target  may  be  in  motion.  Define  Rq  (t)  as  the  direct-path  range 
between  the  ?’th  transmitter  and  jth  receiver,  Rq  (t)  =  ||rJ  —  d*||.  Similarly,  let  R\  (t) 
=  || t  —  d* ||  and  i?2(^)=  ||rJ  — 1||  denote  the  lengths  of  the  first  and  second  legs  of  the 
ijtli  target-path  channel  between  the  ith  transmitter,  target,  and  jth  receiver. 


Figure  8.  Geometry  for  the  ijth  Transmitter-Receiver  Pair 


The  jth  receiver  is  assumed  to  have  an  array  antenna  consisting  of  N3e  identical 
elements,  with  the  nth  element  at  location  r3n  =  r{+<5^,  where  =  rJ  is  the  location  of 
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the  reference  element,  and  83n  is  the  offset  vector  pointing  from  the  reference  element 
to  the  nth  element,  as  shown  in  Fig.  8.  Note  that  S{  =  0.  For  simplicity,  all  arrays 
are  assumed  to  have  an  equal  number  of  array  elements,  i.e.,  N3  =  Ne  for  all  j.  Let 
G3en(x.)  denote  the  element  pattern  of  the  nth  element  in  the  direction  of  x.  For 
simplicity,  it  is  assumed  that  element  patterns  are  equal  across  array  elements  within 
a  given  array,  i.e.,  G3e n(x)  =  G3e(x :)  for  all  n.  Such  a  condition  approximately  holds 
for  the  aligned  array  elements  within  a  planar  array  (e.g.,  GAMMA  [65]),  or  for  the 
azimuthally-isotropic  elements  within  a  circular  array  (e.g.,  PaRaDe  [59]).  Finally, 
let  k3n (x)  denote  the  unit  vector  pointing  from  the  nth  element  of  the  jth  receiver  to 
the  position  x,  i.e., 


ki(x)  = 


—  r-? 


X  -  Yrt 


(1) 


For  x  in  the  far  held,  k),(x)  «  k{(x)  =  kJ(x),  i.e.,  the  unit  vectors  pointing  from  the 
array  elements  to  a  location  in  the  far  held  are  approximately  equal. 


2.2  Signal  Environment 

The  ith  transmitter  generates  a  continuous  wave  (CW)1  narrowband  bandpass 
waveform  ul(t)  with  analytic  representation 

u^t)  =  ui(t)e3“c\  t  e  [0,  T]  (2) 


where  j  is  the  imaginary  unity,  cu*  is  the  carrier  frequency,  T  is  the  coherent  processing 
interval  (CPI),  and  ul(t)  is  the  complex  envelope  with  frequency  representation  U1(lo) 
and  bandwidth  Bl  in  Hertz,  U1(uj)  ~  0  for  all  |cu|  >  ttB1.  The  transmitted  signals 
{■ ul{t )  :  i  =  1 . . .  A)J  are  assumed  to  occupy  non-overlapping  frequency  channels  such 
that  they  are  separable  in  frequency.  The  complex  envelope  ul{t)  is  assumed  to  be 
lrThe  term  continuous  wave  is  used  to  denote  a  signal  with  continuous  transmission. 
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defined  such  that  the  energy  of  ul(t),  denoted  E\  is  given  by 

Ei=  f  \u\t)\2dt  =  [  \ui(t)\2dt  =  T  (3) 

Jo  Jo 

Note  that  this  condition  is  satisfied  for  unit-modulus  signals,  |ti?(t)|2  =  1. 

The  signal  ul(t)  is  presented  to  the  ith  transmitter’s  aperture  and  radiated  into 
the  environment.  This  signal  propagates  to  the  jth  receiver  along  the  direct-path 
(transmitter-to-receiver)  and  target-path  (transmitter-to-target-to-receiver)  channels, 
as  depicted  in  Fig.  8.  It  also  propagates  to  the  jth  receiver  along  multiple  clutter- 
path  (transmitter-to-clutter-to-receiver)  channels.  The  resulting  clutter-path  signals 
may  be  mitigated  by  adaptive  spatial  and  temporal  filtering  as  discussed  in  Sec.  1.1.2; 
consequently,  they  are  ignored  in  the  following.  The  signal  received  at  the  nth  element 
of  the  jth  receiver  is  the  superposition  of  the  direct-path  and  target-path  signals  from 
all  transmitters  within  the  receiver  bandwidth  and  receiver  noise.  It  is  given  by 

Nt  Nt 

ad  W  (*  ~  Td?n(*))  +  W  (*  ~  TtUf))  +  hn  (*)  (4) 

i=  1 

where  ajn(t)  and  <ifn(t)  denote  real  amplitude  scaling  coefficients  associated  with 
the  direct-path  and  target-path  channels,  respectively,  a 13  is  the  target’s  complex 
bistatic  reflectivity  associated  with  the  i jth  bistatic  pair,  r^n(t)  and  r^n(t)  denote 
the  propagation  time  delays  associated  with  the  direct-path  and  target-path  chan¬ 
nels,  respectively,  and  nJn(t)  is  wide-sense  stationary  (WSS)  bandpass  Gaussian  white 
noise  with  power  spectral  density  (PSD)  Pj(co)  =  N0  over  \uj  —  uj3c\  <  nB \  given  a 
receiver  bandwidth  B3  in  Hertz  and  center  frequency  c The  channel  coefficients 
adn(t)  an<4  a%n(J)  account  for  the  composite  amplitude  scaling  of  ul(t)  associated 
with  transmission,  propagation,  and  reception  along  the  direct-path  and  target-path 
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channels,  respectively.  They  are  defined  by 


^,p(rJ)A‘2G;(d-) 

P-rP(t)A-aGi(t) 

(4T)3(fl;(t)ff2(t))2 


(5) 

(6) 


where  Pc*rp  is  the  ith  transmitter’s  effective  radiated  power  toward  x,  A*  =  c//*  is 
the  wavelength  of  the  ith  transmitter,  c  is  the  speed  of  light,  and  /*  =  27T /ufc.  Note 
that  Pq  (t),  R\(t)  and  R32(t)  are  not  expected  to  change  significantly  over  the  interval 
[0,  T\.  Consequently,  let  alj  =  a^(t) |t_Q  and  a1/  =  a^'(t) |(_Q. 

The  signal  s£(t)  is  down-converted  and  channelized  in  frequency  to  extract  the 
complex  baseband  signal  for  each  transmit  channel.  The  resulting  complex  baseband 
signal  for  the  ith  channel,  denoted  s^(t),  is  given  using  (2)  and  (4)  by 


4(f)  =  LPF*  4(f)  e 


oj(0i-puit) 


(7) 


ij  j\  6°  —  cj'It  , 

=  a^e  v  c  d'n 


^-W)«*(t-4(t))+n«(t)  (8) 


where  Qi  in  (7)  denotes  an  unknown  phase  associated  with  the  jth  receiver’s  downcon- 
version  oscillator,  LPF*{-}  denotes  a  low  pass  filter  that  is  matched  to  the  bandwidth 
of  the  ith  transmit  channel,  and  ri%(t)  is  WSS  complex  baseband  Gaussian  white 
noise  with  PSD  P^(oj)  =  N0  over  |u;|  <  7 rBl.  Note  that  the  unknown  phase  63 
in  (7)  indicates  that  the  receivers  are  not  phase-synchronized  (coherent).  Achieving 
coherence  between  widely  distributed  receivers  is  difficult  in  practice.  However,  the 
receivers  are  considered  to  be  time-synchronized  in  the  sense  that  they  are  aligned  in 
time  to  within  a  small  fraction  of  the  resolution  time  of  the  each  transmit  channel, 
i.e.,  to  within  a  small  fraction  of  1  / Bl  for  all  i.  The  direct-path  and  target-path 
components  of  (8)  are  discussed  in  turn. 
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2.2.1  Direct-Path  Signal. 


Let  s^n(t)  denote  the  direct-path  component  of  (8), 

4, nit)  =  a«  u<(t  -  r -«,(<))  (9) 

where  r^n(t)  =  l?^n(t)/c,  and  R^n(t)  is  the  range  between  the  ith  transmitter  and 
the  nth  element  of  the  jth  receiver  array.  The  direct-path  range  can  be  expanded 
into  a  sum  of  two  terms, 


RZn(t)  =  Rij(t)  +  ARi>(t-,di)  (10) 

where  Rq  (t)  is  the  range  between  the  ith  transmitter  and  the  reference  element  of  the 
jth  receiver,  and  AR%(t;dt)  is  the  differential  range  of  the  nth  array  element  with 
respect  to  the  reference  element.  For  d'  in  the  far-held,  it  is  shown  in  Appendix  A 
that  AR%(t;dl)  may  be  approximated 

AR%(t;  d*)  «  -k^d*)  •  S{  (11) 

Note  that  for  d*  in  the  far-held,  kJ(d*)  is  approximately  constant  over  t  e  [0,  T\. 
Therefore,  AR%(dl)  =  A R%(t;  d*)|t=0  is  used  in  the  following.  From  (10),  the  propa¬ 
gation  delay  rjn(t )  can  then  be  expressed  as 

rUt)  =  Td(t)  +  ArS(  d4)  (12) 


where  rlJ  (t)  =  R,Z(t)/c  and  Ar^d*)  =  AR^(dl)/c. 


Inserting  (12)  into  (9), 


3 d,n 


(t)  =  aj  eX^-A^fdO)  ui  _  Tij _  Arij(d<))  e~^ 


'W 


(13) 
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Note  that  the  only  remaining  dependence  on  the  element  index  n  is  in  the  complex 
exponential  e_JW=Ar"  (d!).  Le^  denote  the  phase  of  this  complex  exponential, 

which  can  be  written  using  (11)  as 


C(cf)  A  -w*Ar«(d*) 


~  Ai^(d') 


^  ]  k'(d-)  ■  s; 


(14) 

(15) 

(16) 


Then,  applying  the  narrowband  approximation  to  (13),  in  which  the  complex  envelope 
of  the  signal  is  approximately  constant  across  the  elements  of  the  array  [81,  p.  34], 


4.0  *  u-  (t  -  sf  (i))  e-^ 


(17) 


Next,  consider  a  first-order  approximation  for  Rq  (t)  about  t  —  0, 


t=0+Rii(t) 


+  A  p ij  I  p ij  + 

L  -  -Ztf)  -Tin  L 


t= o 


Then,  using  Tdj  (t)  =  R^(t)/c  in  (17), 

4,0  =  a«eA+«<d'))ui(t„  H*t(t)/c)e-*in%m/e 

w  4eA+«ld'))u^(l  -  fl«/c)t  -  fl«/c)  e-*4W,+sSl<)/s 
=  eA+«,<d'))  „i(at  -  Tf)  e->"A  e»«i 


(18) 


(19) 

(20) 
(21) 


where  rd  =  Rq  /c,  ck  =  1  —  i?QJ/c  is  a  time  scaling  factor,  and  0Jd  is  the  Doppler 
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frequency  defined  as 


d?  =  -  (f )  t?i  (22) 

f  2i r\  (rJ  —  d*)  •  (rJ  —  d'') 

\  A*  /  ||rJ  —  d* || 

Note  that  a  =  1  —  Rq  /c  ~  1  for  the  range  of  expected  Rq  .  Consequently,  it  is  ignored 
in  the  following. 

To  summarize,  the  complex  baseband  direct-path  signal  received  by  the  nth  array 
element  of  the  jth  receiver  originated  by  the  ith  transmitter  is  given  by 


(23) 

t  = 


(a)  (b)  (c)  (d)  (e)  (f) 


(24) 


Equation  (24)  is  comprised  of  six  factors:  (a)  an  amplitude  scaling  factor,  (b)  an 
unknown  oscillator  phase  factor,  (c)  a  differential  carrier  phase  factor,  (d)  a  reference 
carrier  phase  factor,  (e)  the  time  delayed  complex  baseband  transmitted  signal,  and 
(f)  a  Doppler  modulation  factor.  To  simplify  notation,  (24)  may  be  expressed  as 


(25) 


where  7 ^  is  ij th  direct-path  channel  coefficient,  defined  using  (5)  and  (25)  as 


O  A  ij 


Id  =  ad e' 


=  e.^-1 


(4  tt)W)2 


(26) 

(27) 
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2.2.2  Target-Path  Signal. 


Let  s\°n{t)  denote  the  target-path  component  of  (8), 


3 t,n 


(t)  =  o'V 


u'(t  -  T%n(t)) 


(28) 


Following  a  similar  procedure  to  that  used  for  the  direct-path  signal  in  Sec.  2.2.1, 
s\Jn(t)  can  be  approximated  as 


(29) 


where  Tlt 3  is  the  bistatic  delay  with  respect  to  the  ith  transmitter,  target,  and  jth  re¬ 
ceiver, 


r«  4  (L)  (flf  +  ii|)  (30) 

=  ())  (l|t-d‘||  +  ||H'-t||)|(=0  (31) 

d()(t)  is  the  differential  carrier  phase  of  the  nth  element  with  respect  to  the  reference 
element,  defined  as 


<W)  =  -^Ar«(t) 
=  k«(t)-4 


(32) 

(33) 


and  uj\3  is  the  target-path  bistatic  Doppler  shift, 


,  A 

(dt  — 


^  {R\j  +  Rij) 


/2t r 


(t  -  d*)  •  (t  -  d*)  ( rj  -  t)  •  (F  -  t) 


It  —  dd 


F  -  t 


(34) 

(35) 


t= o 
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Similar  to  the  direct-path  signal  in  (25),  (29)  is  comprised  of  six  factors:  (a)  an  am¬ 
plitude  scaling  factor,  (b)  an  unknown  oscillator  phase  factor,  (c)  a  differential  carrier 
phase  factor,  (d)  a  reference  carrier  phase  factor,  (e)  the  time  delayed  complex  base¬ 
band  transmitted  signal,  and  (f)  a  Doppler  modulation  factor.  To  simplify  notation, 
(29)  may  be  expressed  as 


=  7?  ^ ’D  eI"?‘ 


(36) 


where  7^  is  the  ijth  target-path  channel  coefficient,  defined  using  (6)  and  (36)  as 


ij  A  ij  ij  lioi—  aj*rP) 

7/  =  a  JatJ  er\  c  ‘  / 


=  cdV 


®Wer«)  ^ /^p(t)A*2Gi(t) 

{4ir)3(R\Ri)2 


(37) 

(38) 


2.3  Discretization 


The  ith  baseband  signal  at  the  nth  element  of  the  jth  receiver  is  given  by 


sm  =  sln(t)  +  ^n(t)+nm  (39) 

Let  (39)  be  sampled  at  rate  /*  =  1  jT\  =  Bl  Hz  such  that  ,s^ [/]  =  (/T*),  resulting 

in  the  discrete  time  representation 

sn  M  =  S%n\f\  +  st?nR]  +  1  =  0,...,  V  —  1  (40) 
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where  Ll  =  \T /* J  is  the  total  number  of  samples,  and  the  direct-path  and  target-path 
signals  are  given  using  (25)  and  (36)  by 


yj  (di)  ul  [l  -  e°A\  e ivdl  (41) 

e?v?1  (42) 


where  uA  =  ujlj  j f\  and  ultJ  =  oj\j  / f l  are  normalized  Doppler  frequencies  measured  in 
radians  per  sample,  and  tA  =  ta  f *  and  =  Tty/®  are  normalized  delays  measured 
in  samples.  The  noise  sequence  n^[Z]  ~  CJ\f(0,a2),  where  <x2  =  N0Bl  is  the  average 
noise  power,  E  {n]{[l](nlY' [&])*}  =  cr28i-i>8j-ji8i-ki  and  8n  is  the  Kronecker  delta. 

Let  u*  G  Ci!xl  denote  the  sampled  transmit  waveform  vector  such  that  the  Ith 
element  is  defined  [u1]/  =  ul[l\  =  ul(lTls)  for  l  —  0, , . . ,  Ll  —  1.  It  can  be  shown  using 
(3)  that  1 1  u?  1 1 2  =  L\  Specifically, 

lluil2  =  E  l“‘(iU)|2  «  =  f  \u<(t)\2dt  —  f\T  —  V  (43) 

1=0  ^  Jo 

Let  s%>  sd,n >  and  st,n  be  defined  similarly,  i.e.,  [sj?']j  =  s%[l],  [s Jn]j  =  and 

[s((n]/  =  s&[t]  for  /  =  0, . . . ,  L1  —  1.  Furthermore,  define  D L{x)  G  CLxL  by 

Dl(x)  =  diag  ( [ej{0)x,  ejWx,  •  •  •  ,  e3{L~1)x ] )  (44) 


where  diag(x)  for  x  G  CL  is  the  L  x  L  square  matrix  with  diagonal  elements  x  such 
that  the  ( n,n)th  element  [diag(x)]n>n  =  [x]n.  Finally,  let  G  CLxL  denote  the 
unitary  discrete  Fourier  transform  (DFT)  matrix  such  that  the  ( m,n)th  element  is 


[W] 


^)mn 


(45) 
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for  m  =  0, . . . ,  L  —  1  and  n  =  0, . . . ,  L  —  1. 

Then,  the  direct-path  signal  (41)  can  be  represented  in  vector  form  by 

(Wf,Dii(-2</i‘)WI.)  u‘  (46) 

To  simplify  (46),  define  the  delay-Doppler  operator  T>(£,  v)  G  CLIxiI  as 

T>(£,  v)  =  DLi(z/)W^1Dii(— 27r£/ UyW Li  (47) 

which  is  the  matrix  representation  of  the  delay  by  £  and  Doppler  shift  by  v  operators 
applied  to  a  length- U  time  domain  sampled  signal.  Note  that  (47)  applies  a  circular 
delay  of  i  samples.  In  many  cases  the  wrapping  of  the  tail  end  of  the  signal  is 
insignificant  because  the  delay  is  only  a  small  fraction  of  the  total  signal  duration, 
i.e.,  f  Alternatively,  such  wrapping  may  be  avoided  by  zero-padding  the  received 
signals  to  account  for  the  longest  expected  delay.  Note  also  that  T>(£i,  z/1)T>(£2,  n2)  = 
V{£  1  +  £2,1^1  +  v2)i  and  VH(£,u)  =  T>(—£,—u).  Consequently,  the  delay-Doppler 
operator  is  unitary,  i.e.,  D-1(£,  v)  =  VH[£ ,  v)  such  that  VH{£ ,  v)V(£ ,  v)  =  I Li,  where 
1L  is  the  L  x  L  identity  matrix.  Using  (47)  in  p(46),  the  direct-path  signal  becomes 

=  u‘  (48) 

Similarly,  the  target  signal  (42)  can  be  represented  in  vector  form  as 

4n  =  ^e^V(e‘,V»)^  (49) 


Therefore,  the  ith  baseband  signal  at  the  nth  element  of  the  jth  receiver  array  in 


(40)  can  be  written  in  vector  form  using  (48)  and  (49)  as 


4  =  7?  e**<d‘l  ©(<?,  v»)  u‘  +  T«  e’«<«  ©(<« ,  tf )  u*  +  n« 


(50) 


where  ~  CAf(0Li,  a 2IL;)  and  0Li  is  the  zero  vector  of  length  Ll. 


2.4  Summary 

This  chapter  has  presented  a  detailed  derivation  of  the  PMR  signal  model  that 
underlies  the  remainder  of  this  dissertation.  In  particular,  Eq.  (50),  which  represents 
the  discretized  complex  baseband  signal  at  the  nth  array  element  of  the  jth  receiver 
associated  with  the  ith  transmitter  channel,  is  the  starting  point  for  the  mathematical 
developments  in  Chapters  III  and  IV.  It  also  underlies  signal  models  for  AMR  and 
PSL  sensor  networks  that  are  presented  in  Chapter  V  as  part  of  the  unified  detection 
framework.  The  form  of  Eq.  (50)  foreshadows  several  of  the  challenges  associated 
with  passive  radar  detection.  As  discussed  in  Sec.  1.1,  the  simultaneous  presence 
of  the  direct-path  and  target-path  signals  in  (50)  complicates  detection  due  to  the 
large  power  ratio  between  the  direct-path  and  target-path  signals,  |  1 2  / 1  At 7 1 2  V>  1, 

that  is  typical  in  passive  radar  scenarios.  In  addition,  the  presence  of  receiver  noise 

guarantees  that  the  direct-path  signal  will  not  be  able  to  be  isolated  without 
distortion.  Finally,  the  delays  and  Doppler  shifts  of  the  target-path  and  direct-path 
signals,  encoded  by  the  delay-Doppler  operators  and  T>(ttJ ,v\3)  require 

careful  consideration.  These  challenges  are  addressed  in  Chapters  III  and  IV. 
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III.  Passive  MIMO  Radar  Detection  without  References 


This  chapter1  addresses  the  problem  of  target  detection  in  PMR  networks  when 
direct-path  signals  are  unavailable,  shown  in  Fig.  9.  As  discussed  in  Sec.  1.1.2,  direct- 
path  signals  are  typically  the  dominant  signals  (and  interferers)  in  passive  radar  ap¬ 
plications  [30].  However,  it  is  possible  they  might  not  be  present  if  there  is  physical 
blockage  of  the  linc-of-sight  paths  between  transmitters  and  receivers,  or  if  illumina¬ 
tors  were  to  utilize  highly  directive  transmission  and/or  null  the  receiver  directions. 
In  such  cases,  the  absence  of  direct-path  signals  would  preclude  the  conventional  ap¬ 
proach  to  detection  described  in  Sec.  1.4.1,  and  would  necessitate  a  fundamentally 
different  type  of  processing.  This  has  given  rise  to  recent  theoretical  interest  in  target 
detection  in  PMR  networks  when  direct-path  signals  are  unavailable  [10,  84], 

Wang  and  Yazici  present  a  centralized  detector  for  this  problem  that  is  constrained 
to  be  linear  in  the  measurements  [84],  Their  formulation  assumes  that  targets  scatter 
isotropically,  which  is  unlikely  to  hold  for  complex  targets  in  multistatic  geometries 
[88].  Their  formulation  also  assumes  phase-synchronized  (coherent)  receivers,  which 

is  a  stringent  requirement  in  distributed  applications.  Bialkowski  et  al.  derive  a 
lrThe  material  in  this  chapter  appears  in  [38]. 


Figure  9.  Passive  MIMO  Radar  without  Direct-Path  Signals 
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generalized  likelihood  ratio  test  (GLRT)  without  constraining  the  detector  to  be  linear 
in  the  measurements  [10].  The  resulting  non-linear  detector  implicitly  accounts  for 
non-isotropic  target  reflectivity  and  non-coherent  receivers.  However,  neither  [84]  nor 
[10]  address  the  probability  distributions  of  their  detection  statistics. 

In  this  chapter,  the  results  of  [10]  are  extended  in  several  ways.  First,  the  for¬ 
mulation  of  [10]  is  extended  from  one  transmitter  to  multiple  transmitters,  and  from 
single-channel  receivers  to  multichannel  (array)  receivers.  The  GLRT  is  then  derived. 
These  extensions  significantly  enhance  probability  of  detection  and  reduce  detection 
ambiguity.  Second,  using  recent  results  from  random  matrix  theory,  the  exact  dis¬ 
tribution  of  the  test  statistic  is  identified  under  both  hypotheses.  To  the  author’s 
knowledge,  this  is  the  first  appearance  of  these  distributions  in  the  context  of  radar 
detection.  These  distributions  show  that  detection  sensitivity  is  only  a  function  of 
the  number  of  transmitters,  the  number  of  receivers,  the  number  of  signal  samples, 
and  the  average  input  SNR.  Numerical  examples  demonstrate  important  properties 
of  the  detector,  namely,  that  (a)  receivers  and  transmitters  contribute  asymmetrically 
to  detection  sensitivity,  which  is  uncommon  in  radar  applications,  and  (b)  integra¬ 
tion  gain  grows  non-coherently  with  increasing  signal  length.  Third,  an  equivalance 
is  demonstrated  between  PMR  detection  without  direct-path  signals  and  detection 
in  PSL  networks.  This  generalizes  PSL  to  the  scenario  in  which  multiple  transmit¬ 
ters  are  co-located  on  the  “source”  platform,  and  establishes  equivalence  between 
detection  ambiguity  and  “source”  localization.  Finally,  the  ambiguity  performance 
of  the  proposed  detector  is  investigated,  and  it  is  shown  that  salient  features  of  the 
detector  ambiguity  function  can  be  explained  in  terms  of  the  time-difference  of  arrival 
(TDOA),  frequency-difference  of  arrival  (FDOA),  and  angle-of-arrival  (AOA)  of  the 
target  signals. 
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3.1  Signal  Model 


Consider  a  PMR  network  like  the  one  shown  in  Fig.  9.  Let  there  be  Nt  transmit¬ 
ters  and  Nr  receivers.  Fig.  10  depicts  the  geometry  and  signal  environment  of  the 
ij th  bistatic  pair,  consisting  of  the  ith  transmitter  and  jth  receiver.  The  position 
and  velocity  of  the  ith  transmitter  are  denoted  by  d‘  and  dl,  and  the  position  and 
velocity  of  the  jth  receiver  by  rJ  and  rJ.  The  spatial  states  of  all  transmitters  and 
receivers  are  assumed  known.  Assume  the  transmitted  signals  are  narrowband  and 
separable  in  frequency,  and  that  each  receiver  is  equipped  with  an  Aj-element  array. 
The  incident  signal  at  the  nth  element  of  the  jth  receiver  is  channelized  in  frequency, 
demodulated  to  baseband  and  sampled  in  time.  Let  s%  G  Cilxl  denote  the  length- L* 
sampled  complex  baseband  signal  on  the  ith  frequency  channel  of  the  nth  element  of 
the  jth  receiver  array.  For  simplicity,  it  is  assumed  that  U  =  L  V  i,  i.e. ,  all  received 
signals  have  the  same  length  L.  An  expression  for  is  given  by  Eq.  (50)  in  Chap¬ 
ter  II.  Assuming  a  target  at  position  t  and  velocity  t,  and  ignoring  the  direct-path 
signal  in  (50),  s  j  is  given  by 

s«  =  7«  e«t>  v(e=,  v\i)  u<  +  n«  (51) 

where  u*  G  CLxl  is  the  complex  baseband  signal  emitted  by  the  ith  transmitter; 
T>(£]3 ,  ly7)  G  CLxL  is  the  unitary  linear  operator  that  accounts  for  the  delay  and 
Doppler  shift  imparted  to  the  ith  transmit  signal  as  it  propagates  to  the  jth  receiver 
along  the  target-path  channel,  defined  by  Eq.  (47)  in  Sec.  2.3;  £lt3  is  the  target-path 
propagation  delay  in  samples;  ry7  is  the  target-path  Doppler  shift  in  radians  per  sam¬ 
ple;  'd^( t)  is  the  differential  phase  at  the  nth  array  element  with  respect  to  a  reference 
element  due  to  plane  wave  propagation  from  the  target  direction,  defined  by  Eq.  (33) 
in  Sec.  2.2.2;  7^  is  a  complex  channel  coefficient  that  accounts  for  the  composite 
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Figure  10.  The  ijth  Bistatic  Pair  in  a  PMR  Network  without  Direct-Path  Signals 


scaling  of  iT  associated  with  the  ijth  target-path  channel,  defined  by  Eq.  (38)  in 
Sec.  2.2.2;  e  CLxl  is  circular  Gaussian  noise  distributed  as  CN(0l,  ct2Il)  with 
known  variance  cr2;  and  j  =  y/—l.  Here,  0 l  denotes  the  length- L  zero  vector,  and 
If  is  the  L  x  L  identity  matrix.  Noise  is  assumed  to  be  independent  across  transmit 
bands,  receivers,  and  array  elements,  i.e.,  Ejnjf =  a25n-mSi-kSj-ilL,  where 
(■)H  is  the  Hermitian  transpose  and  Sx  is  the  Kronecker  delta.  The  transmit  signal  u* 
is  defined  such  that  ||u*||2  =  L.  Additional  detail  on  the  terms  in  (51)  is  provided  in 
Chapter  II.  To  simplify  notation,  T>\3  is  used  in  place  of  T>(t3 ,  G3)  in  the  following. 

Fig.  10  shows  that  the  channel  coefficient  7^  may  be  expressed  as  the  product  of 
two  coefficients  that  correspond  to  the  two  legs  of  the  bistatic  target-path  channel, 
7 \3  =  7t  1 OY2 )  where  7)  [  and  7y,2  are  defined  using  (38)  as 


ryi  (t)  —  p-MR i(t)/c  /  -PerpW 

7w()  (t» 

7^2  (t)  =  cdv(0i-^(t)/c) 


\i2Gl(t) 


\  (47Ti?J2(t))' 


(52) 

(53) 


cu*  =  27 r/*,  /*  is  the  ith  transmitter’s  carrier  frequency  in  Hertz,  -R^(t)  =  ||t  —  d*||, 
R32(t)  =  Hr-7  —  1 1| ,  c  is  the  speed  of  light,  Pe*rp(t)  is  the  ith  transmitter’s  effective 
radiated  power  in  the  direction  of  t,  a*7  is  the  target’s  bistatic  reflectivity  as  observed 
by  the  ij th  Tx-Rx  channel,  07  is  a  random  phase  associated  with  the  jth  receiver, 
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X1  =  cj fl  is  the  wavelength  of  the  ith  transmitter,  and  GJe(t)  is  the  array  element  pat¬ 
tern  in  the  direction  of  t.  Similarly,  the  delay-Doppler  operator  T>\3  may  be  expressed 
as  the  composition  of  two  unitary  linear  operators  that  correspond  to  the  two  legs  of 
the  bistatic  target-path  channel,  i.e.,  T>\3  =  T>)32T>\  l ,  where  Vltl  and  T>\32  characterize 
the  delay  and  Doppler  shift  of  the  first  and  second  legs  of  the  target-path  channel, 
respectively. 

Let  a*J  (x)  £  CNeXl  denote  the  spatial  steering  vector  in  the  direction  of  x, 

a,:i(x)=  e^(x)  ...  T  (54) 

The  concatenation  of  the  time  series  vectors  from  all  Ne  elements  of  the  jtli  receiver 
in  the  ith  frequency  band,  s*-7  =  [(s*/)T  •  •  •  (s^JT]T  G  CNeLxl,  can  be  written  as 

s13  =  M ]3  u*  +  niJ  (55) 

where  iT-7  =  [(nf)T  •  •  •  (n %JT]T  £  CNeLxl,  and  the  matrix  G  CNeLxL  is  dehned 

Mt  =  It  (at  ®vt)  (56) 

where  a)J  =  ay  (t)  for  brevity,  and  Cg>  is  the  Kronecker  product. 

3.2  Surveillance-Surveillance  GLRT 

Consider  testing  for  the  presence  of  a  target  within  a  position- velocity  cell  (p,  p), 
termed  the  '‘cell  under  test ,’  where  p  and  p  are  the  hypothesized  emitter  position 
and  velocity,  respectively.  This  detection  problem  may  be  formulated  as  a  binary 


34 


hypothesis  test  between  alternative  ("Hi)  and  null  (Ho)  hypotheses: 


H\  :  s'-'  =  M;-/  u'  +  n  ' 

Ho  :  s'-'  =  !i;' 


(57) 


for  i  —  1 . . .  Nt  and  j  —  1 .. .  Nr.  In  (57),  the  matrix  M*-7  e  qn^lxl  jg  (]efpierj 


Mij 

p 


(58) 


where  7 is  the  target-path  channel  coefficient  associated  with  the  hypothesized  po¬ 
sition  p,  a)/  =  alJ(p)  for  brevity,  Vlpj  =  V(tJ ,  z/*J),  is  the  bistatic  delay  associated 
with  the  hypothesized  position  p  in  samples,  and  uft  is  the  bistatic  Doppler  shift 
associated  with  the  cell  under  test  (p,  p)  in  radians  per  sample.  Let  s*  denote  the 
concatenation  of  all  receiver  measurements  associated  with  the  ith  transmitter,  and 
let  s  denote  the  concatenation  of  all  s*  across  transmitters, 


s'  =  [(s'1)7’, . . . ,  (s'w,)r]r  e  cA'rtxl 
S  =  [(s ‘)T,.. .  ,(sNrf}T  e  cA',A''lxl 


(59) 

(60) 


Similarly,  let  7),  denote  the  vector  of  coefficients  associated  with  the  ith  transmitter, 
and  let  7p  denote  the  concatenation  of  all  Yp  across  transmitters, 


[7“  .  ■■l'vNr]T  £  C"’u 

(61) 

[(7jf  ...(7Dt]T  ec“” 

(62) 

Finally,  let  u  =  [(u1)11 . . .  (uArt)T]T  e  CNtLxl. 

Due  to  the  independence  of  the  receiver  noise  across  transmitter  channels,  the 
conditional  probability  density  function  (PDF)  of  s  under  Hi,  Pi(s|7p,  u),  factors 
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according  to 


(s  i  ipi u) = iiri(B*i7j,u* 


where 


1  f  l 

Pi(si  I  7p,  u*)  =  7- ^2WrLexp  -  —J2  ||s<J'  -  M;/u;  2 


Similarly,  the  PDF  of  s  under  "H0 ,  Po(s)>  is  given  by 


^0(S)  =  (7ra2)jVfiVrLeXP{-ipllSl|2} 


In  this  formulation,  the  transmit  signals  u  and  channel  coefficients  7p  are  considered 
deterministic  and  unknown.  Thus,  hypothesis  "Hi  is  composite  because  Pi(s|7p,u) 
is  parameterized  by  u  and  7p.  Therefore,  the  GLRT  is  derived,  which  replaces  each 
unknown  parameter  by  its  maximum  likelihood  estimate  (MLE)  in  the  likelihood 
ratio  test  [49].  This  derivation  follows  the  approach  of  Bialkowski  et  al.  in  [10],  yet 
extends  their  formulation  from  one  transmitter  to  multiple  transmitters,  and  from 
single-channel  receivers  to  multichannel  array  receivers. 

Let  Zi(7p,u|s)  =  log  pi  (s  |  7  ,  u)  denote  the  log- likelihood  function  under  'H\. 
Similarly,  let  /q(s)  =  logpo(s)-  Then,  the  GLRT  may  be  written  as 


max  Zi(7  ,  u  |  s)  —  l0(s)  ^  k  (66) 

{UpM  Ho 

where  k  is  chosen  to  maintain  a  desired  probability  of  false  alarm.  From  (63), 
hi'Jp,  u  |  s)  may  be  written  as 


hi'Yp,  u  |  s)  =  ^ri(7;,u*|sl 
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where,  from  (64)  and  ignoring  an  additive  constant, 


Nr 


i\(%,  ui s')  =  -4  E  IIs”  -  M”u‘ll2 

°  i= i 


(68) 


Substituting  (68)  into  (56), 


Nr 


i\(%,  u‘  |s*)  =  -4E  IIs”  -  ^”<a”  ®  ®?>‘f 

a  J=1 


(69) 


From  (69),  the  MLE  of  7*J,  7*J,  is  given  by 


/p 


((a^  ® 

|(ap  ®  VlJ )u*||2 


This  may  be  simplihed  to 


_ 

ip 


u 


where 


=  (Z>? 


>*s? 


(70) 


(71) 


(72) 


and 


1 

7^ 


(V e 

E 

n=l 


: 

Lap  J , 


(73) 


The  notation  [x]n  denotes  the  nth  element  of  x.  Note  that  s*J  may  be  interpreted  as  a 
surveillance  channel  for  the  7th  transmit  channel  that  is  formed  at  the  jtli  receiver  by 
beamforming  in  the  direction  of  p,  and  slJ  is  this  channel  after  removal  of  the  delay 
and  Doppler  shift  resulting  from  bistatic  propagation  with  respect  to  the  cell  under 
test  (p,  p).  Substituting  (71)  for  Yp  in  (69),  and  simplifying,  gives 


u1 


(74) 
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where  <E>*  =  [s*1, . . . ,  slsNr]  G  CLxNr.  Let  Ai(-)  denote  the  largest  eigenvalue  of  its 
matrix  argument,  and  let  vi(-)  denote  the  associated  eigenvector.  Then,  the  Rayleigh 
quotient  in  (74)  achieves  its  maximum  value,  Ai(<h*<f>(ff),  when  u  =  vi(<3>*<L*w)  [44,  p. 
176].  Therefore  u*  =  V!(<I>*$*H),  and 

|s‘)  =  (llsil2  -  A, (*■<))  (75) 

Noting  that  =  Ai(<h*w<h*),  and  that  typically  Nr  <C  L ,  it  is  more  compu¬ 

tationally  efficient  to  consider  the  Gram  matrix  G*s  =  (<3>*)H<3>*  G  CNrXNr,  giving 

iUi-y  W  |  S‘)  =  (Us4!!2  -  A,(GD)  (76) 

Therefore,  using  (67)  and  (76), 

1  Nt 

h(%,u  |s)  =  -^l|s||2+(^)EAi(GL)  (77) 

2—1 

It  can  similarly  be  shown  under  H0  that 

<o(s)  =  -U|s||2  (78) 

<7Z 


Using  (77)  and  (78)  in  (66),  the  resulting  GLRT  is  given  by 


Mss 


CT- 


Nt 

EAi<Gt) 


2=1 


Hi 

% 

Ho 


hi 


(79) 


The  statistic  £,ss  in  Eq.  (79)  is  termed  the  surveillance-surveillance  GLRT  (SS-GLRT) 
because  the  elements  of  the  matrices  {G*s  :  i  —  1 . . .  Nt}  consist  of  the  pairwise  inner 
products  of  delay-Doppler  compensated  surveillance  signals. 
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3.3  Distribution 


In  this  section,  PDFs  are  presented  for  the  test  statistic  in  (79)  under  'H\  and 
T-Lq.  This  development  considers  the  presence  or  absence  of  a  target  within  the  cell 
under  test  (p,  p),  where  p  and  p  are  the  hypothesized  target  position  and  velocity, 
respectively.  Under  this  condition,  the  surveillance  signal  in  (72)  reduces  to 

s\j  =  bx  /,('u;  +  (80) 

where  =  7 ^  y/We,  bi  =  1  under  "Hi  and  bi  =  0  under  "H0,  and  is  distributed  as 
CA/”(0l,  <t2Il).  Note  that  /i(J  is  the  channel  coefficient  y*-7  scaled  by  a  beamforming 
gain  of  y/Ng.  Then,  s*J  is  distributed  as 

slj  ~  CAf(h  fjj u\  (81) 


To  begin,  note  that  £ss  =  where 

e.  =  Ai(^G‘,)  (82) 

Because  the  receiver  noise  is  independent  across  transmit  channels,  the  PDF  of  t;ss, 
denoted  p6s(£),  is 


Pi.. (0  =  [pIs  * pL  *■■■* pZ]  (0  (83) 

where  *  denotes  convolution,  and  p\sX0  is  the  PDF  of  The  problem  is  then 
finding  p\aa(£)  under  both  hypotheses. 
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3.3.1  Alternative  Hypothesis. 


Using  (80),  the  Gram  matrix  G*s  under  the  alternative  hypothesis,  Hi,  becomes 


G'„  = 


(a 


i  1 


u*  +  n 


il\H 


tfu'+h's 


"■Vs 


rui+nisNr~ 


(84) 


(/i*Arpu*+n*Ar’-)J/J 

=  (u®(M®jT  +  N‘)"  +  N*)  (85) 

' - v - - V - " 

=  (*i)»  =  s* 

where  =  [/x*1, . . . ,  nlsNr]T  and  N*  =  [n*1, . . . ,  h\Nr].  Thus,  the  columns  of  ($*)'ff 
are  independent  Ar-variate  complex  Gaussian  vectors,  where  the  kth  column  is  dis¬ 
tributed  as  CA/"((/i*)*[u*]^,  cr2Ijvr)  •  Consequently,  G*s  is  a  non-central  uncorrelated 
complex  Wishart  matrix  [47,  90],  denoted  G*s  ~  Wjvr(L,  E*,12®),  where  X*  =  a2lNr 
is  the  covariance  matrix  of  the  columns  of  ($*)■**,  and 


12'  4  (S 


E{(41)h}E{4*}  = 


IM 


(Ai)*(  A 


i 


(86) 


is  a  rank-1  non-centrality  matrix  with  non- zero  eigenvalue  Qs  given  by 


L\\vT 

u2 


(87) 


and  eigenvector  (/z®)*  =  (/tx®)*/||/x®  ||. 

In  general,  the  PDF  of  the  largest  eigenvalue  0  of  a  matrix  X-1X,  where  X  ~ 
W s(t,  S,  12)  with  rank-1  noncentrality  matrix  12,  is  given  by  [48]  as 


/iW>;  M,C) 


e"c|4’(^)|tr(»-1(^)#(^))  U(<t>) 

r^-s  +  tnGrii-^ris-*) 


(88) 


where  Q  is  the  largest  eigenvalue  of  12,  T(-)  is  the  gamma  function  [3,  Eq.  (6.1.1)], 
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tr(-)  denotes  the  trace,  |  •  |  denotes  the  determinant,  U(-)  is  the  unit  step  function, 
<h(0)  is  the  s  x  s  matrix  with  entries 


{0t  me  4>  0p,1(j^  _  s  _|_  i-  £0)  n  =  l 

(89) 

01 -+s-m-n e-4>  n  >  1 

o-Fi(-;  •)  denotes  the  generalized  hypergeometric  function  pFq(ai  ■  ■  ■  ap;bi  ■  ■  ■  bq ;  z)  [29, 
Eq.  (9.14.1)]  with  p  =  0  and  q  =  1,  \I/(0)  is  the  s  x  s  matrix  with  entries 

I  Jo  y^e-y  oF^t-s+l;  ( y)  d y  n  =  1 
[*(0)kn  =  (90) 

7(i  +  s  —  m  —  n  +  1,  0)  n  >  1 

and  7(-,  •)  is  the  incomplete  gamma  function  [3,  Eq.  (6.5.2)].  Therefore,  the  PDF  of 
under  Hi  is 

p|ss(e;H1)  =  /1(^;lVr,L,C:j  (91) 

3.3.2  Null  Hypothesis. 

Under  the  null  hypothesis,  Ho,  G*s  in  (85)  reduces  to 

GL  =  (Nt)HNj  (92) 


which  is  a  central  uncorrelated  complex  Wishart  matrix,  denoted  Glss  ~  W^E*), 
with  S*  =  (y2l^r.  In  general,  the  PDF  of  the  largest  eigenvalue  0  of  a  matrix  E_1X) 
where  X  ~  Ws(i,  E),  is  given  by  [48]  as 


/o(0;s,t) 


nUi  T(t  -  *  +  i)r(s  -  *  + 1) 


(93) 
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where  ^c((f>)  and  3?c(</>)  are  s  x  s  matrices  with  elements 


[^,c(0)]m,n  =  'y(t  +  s-m-n  +  l,  4>)  (94) 

[<S>c(<t>)]m,n  =  0*+™^  (95) 

Therefore,  the  PDF  of  £lss  under  Hq  is 

pi.(Z;n0)  =  f0&Nr,L)  (96) 

3.3.3  Dependence  on  SNR. 

The  PDF  (£;  Hi)  in  (91)  depends  on  the  transmit  signal  u*  and  channel  coeffi¬ 
cients  7*  =  [i\/ \[We  through  only  the  received  signal  length  L  and  the  non- centrality 
parameter  Qs  in  (87),  which  may  alternatively  be  expressed  as 

C  =  iVSNR‘yg,  (97) 

where 

SNR;vg=||^||2/(iVra2)  (98) 

is  the  average  input  SNR  associated  with  the  ith  transmitter  after  surveillance  channel 
formation.  Consequently,  the  PDF  P£aa(£',Hi),  which  is  the  convolution  of  the  indi¬ 
vidual  plag(£’,  Hi)  across  transmitters,  is  a  function  of  only  Nt,  Nr,  L ,  and  {SNR,!lvg}. 
Thus,  detection  performance  does  not  depend  on  the  specific  structure  of  the  transmit 
waveforms  {u*},  but  only  their  energies  and  lengths.  In  addition,  it  does  not  depend 
on  the  specific  input  SNRs  at  each  receiver  channel;  rather,  it  depends  on  only  the 
average  SNR  for  each  transmitter  across  receivers. 
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3.4  Interpretation 


This  section  discusses  the  relationship  between  PMR  and  PSL  sensor  networks. 
PSL  sensor  networks  are  comprised  of  geographically  separated  receivers  that  detect 
and  localize  targets  by  intercepting  and  processing  target-emitted  signals.  An  exam¬ 
ple  PSL  sensor  network  is  illustrated  in  Fig.  11.  Comparison  with  Fig.  9  suggests  that 
PMR  sensor  networks  without  direct-path  signals  may  be  interpreted  as  PSL  sensor 
networks.  In  particular,  the  received  signal  in  (51)  may  be  expressed  in  terms  of 
the  channel  coefficients  and  delay-Doppler  operators  of  each  leg  of  the  target-path 
channel,  shown  in  Fig.  10,  as 


4 = 


u  +  n; 


(99) 


and  then  simplified  to 

+  (100) 

where  u*  =  l  T>lt  x  u*.  Eq.  (100)  states  that  the  received  signal  in  a  PMR  scenario 
is  equivalent  to  the  signal  that  would  have  been  received  if  the  signal  u*  had  been 
emitted  by  a  transmitter  at  the  target  state  (t,t).  This  hypothetical  situation  is 
illustrated  in  Fig.  12.  In  this  figure,  the  first  leg  of  the  target-path  channel,  which  is 
faded,  is  replaced  by  the  emission  of  u*  from  the  target.  Note  that  u*  is  simply  the 
incident  signal  at  the  target  after  scaling  and  delay-Doppler  compensation  associated 
with  one-way  propagation  from  the  transmitter  at  (d*,d*)  to  the  target  at  (t,t). 
Note  also  that  the  bistatic  reflectivity  ay",  which  typically  differs  between  transmitter- 
receiver  pairs  in  multistatic  geometries,  is  incorporated  into  the  channel  coefficient 
7^2 (t)  according  to  (53).  This  result  is  intuitive  because  a  target  may  be  interpreted 
as  “emitting”  each  incident  signal  via  the  scattering  process,  where  the  reflectivity 
alJ  represents  the  “antenna  gain”  in  the  direction  of  the  jth  receiver. 
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Figure  11.  A  PSL  Sensor  Network  with  Three  Receivers 


One  consequence  of  this  interpretation  is  that  knowledge  of  the  transmitter  states 
{(d*,d*)}  is  unnecessary  for  detection  in  PMR  networks  that  lack  direct-path  ref¬ 
erences.  This  can  be  seen  for  £ss  by  examining  the  jkth  element  of  G*s.  Let 
—  T>p2Vlp  l,  where  Vp  l  and  Vp2  denote  the  delay-Doppler  operators  of  the  first 
and  second  legs  of  the  target-path  channel  with  respect  to  the  ith  transmitter,  jth  re¬ 
ceiver,  and  hypothesized  state  (p,p).  Then 


[GtR=sf? 

(101) 

=  (vfsi>)H(vfs?) 

(102) 

= 

(103) 

=  {Tf£*,!)a(v%4) 

(104) 

Therefore,  may  equivalently  be  defined  to  consider  only  one-way  delay-Doppler 
compensation  with  respect  to  (p,p), 


ij  _  pO  aij 


p,  2 


(105) 


because  any  delay-Doppler  compensation  with  respect  to  the  first  leg  of  the  target- 
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< 


Figure  12.  Interpretation  of  PMR  without  References  as  a  PSL  Scenario 


path  channel  cancels  in  calculation  of  [G®s]  fe.  Note  that  Wang  and  Yazici  make  the 


same  observation  with  regard  to  calculation  of  their  statistic  [84].  Note  also  that 
the  GLRT  statistic  £ss  reduces  to  the  statistic  of  Vankayalapati  and  Kay  derived 
for  detection  in  PSL  sensor  networks  in  [83]  when:  single-channel  receivers  are  used 
(Ne  =  1);  the  target  emits  on  a  single  transmit  channel  (Nt  =  1);  and  (105)  is  used 
for  slJ .  This  confirms  the  equivalence  between  PSL  sensor  networks  and  PMR  sensor 
networks  that  lack  direct-path  references. 

Additional  insight  into  this  equivalence  is  provided  by  examining  the  ambiguity 
characteristics  of  the  proposed  detector,  i.e.,  examining  how  C,ss  varies  with  mismatch 
between  the  hypothesized  cell  under  test  (p,  p)  and  the  true  target  state  (t,t).  If 
there  is  mismatch  between  t  and  p,  then  the  target  angle-of-arrival  is  possibly  mis¬ 
matched  to  the  hypothesized  angle-of-arrival,  resulting  in  a  mismatch  loss.  In  this 
circumstance,  the  surveillance  signal  slJ  in  (73)  may  be  expressed  as 


(106) 


where  =  /3pt  nlJ  and  /3pt  =  /Ne  quantifies  the  angular  mismatch  loss.  Note 

that  \j3pt\  <  1,  and  f3pt  =  1  iff  =  a1/,  i.e.,  when  p  and  t  are  at  the  same  angle  with 
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respect  to  the  jth  receiver.  Inserting  (106)  into  (104),  and  ignoring  receiver  noise, 


a 


fj^U 


‘"vfv^w^v? 


u 


(107) 


Recall  that  Vf  =  T>(ttj ,  up*),  where  V(£,u)  is  the  delay-by-£  and  Doppler  shift- by- 
v  operator  defined  by  (47)  in  Sec.  2.3,  t/  is  the  target-path  propagation  delay  in 
samples,  and  is  the  target-path  Doppler  shift  in  radians  per  sample.  Similarly, 
Vp2  =  T>(P  2i  ,Jp,2 ) i  where  £Jp  2  is  the  direct-path  delay  between  the  hypothesized  state 
and  the  jth  receiver  in  samples,  and  vp2  is  the  corresponding  direct-path  Doppler 
shift  in  radians  per  sample.  Then,  (107)  may  be  expressed  after  simplification  as 


[G*„]jt  =  AffiV"©)  A^-Af^A^-A^V  (108) 

=  Af  (A4‘  -  A (?,  A,,*’!*  -  Ar!Jt)  (109) 

where  A£Jpk  and  Aftk  are  the  TDOAs  of  the  hypothesized  and  actual  target  signals 
with  respect  to  the  jth  and  kth  receivers,  respectively, 


A#  —  pi  pk 

LAtp  ~  lp,2  tp,2 

(110) 

«se  so 

1 

< 

(111) 

A Vp’jb  and  Av]'3k  are  the  FDOAs  of  the  hypothesized  and  actual  target  signals  with 
respect  to  the  ith  transmit  channel  and  jth  and  kt\i  receivers,  respectively, 


^  =  <2-^2  (112) 

A<4Ji  =  V?-V*  (113) 
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and  x*(A£,  Au)  is  the  ambiguity  function  (AF)  of  u*, 

x‘(Af,AO  =  ^[u-],[u-];+i,^  (114) 

1=0 

Noting  that  \1{A^  Au)  peaks  when  its  arguments  are  zero  leads  to  the  following 
iso-TDOA  and  iso-FDOA  conditions: 


< 

< 

(115) 

Az/’A  =  A  z/p 

P  L 

(116) 

In  words,  [G*J  -fe  peaks  when  (a)  the  hypothesized  TDOA  equals  the  actual  target 
TDOA  and  (b)  the  hypothesized  FDOA  equals  the  actual  target  FDOA.  This  is 
representative  of  distributed  PSL  networks,  which  localize  targets  along  contours  of 
constant  TDOA  and  FDOA  [16,  86]. 

The  preceding  analysis  considers  only  the  elements  of  G*s,  and  it  is  not  imme¬ 
diately  clear  how  these  entries  affect  the  maximum  eigenvalue  of  G*s.  However,  the 
simulation  results  in  Sec.  3.5.2  confirm  that  the  salient  ambiguity  properties  of  £ss, 
which  is  a  function  of  these  eigenvalues,  can  be  explained  in  terms  of  these  iso-TDOA 
and  iso-FDOA  conditions. 

3.5  Simulations 

This  section  illustrates  the  detection  and  ambiguity  performance  of  the  proposed 
detector  via  numerical  simulation.  Numerical  challenges  associated  with  calculating 
the  distributions  presented  in  Sec.  3.3  are  also  discussed. 
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3.5.1  Detection  Performance. 


This  section  illustrates  how  detection  performance  varies  with  SNR^  and  the 
system  parameters  Nt,  Nr,  and  L.  Recall  from  Sec.  3.3.3  that  the  distributions  of 
£ ss  are  functions  of  only  these  four  quantities;  accordingly,  these  quantities  determine 
detection  performance.  The  influence  of  Nt  and  Nr  on  detection  is  discussed  first, 
followed  by  L. 

3.5. 1.1  Number  of  Transmitters  and  Receivers. 

Fig.  13  depicts  probability  of  detection  (Pfi)  curves  as  a  function  of  SNR^  for  a 
scenario  with  L  =  100  samples  per  received  signal  and  a  varying  number  of  trans¬ 
mitter  and  receivers2.  Both  predicted  and  simulated  Pd  curves  are  shown  for  each 
transmitter-receiver  configuration.  The  predicted  curves  are  calculated  according  to 
the  distributions  presented  in  Sec.  3.3.  The  simulated  curves  are  calculated  using  105 
trials  under  "Ho  to  determine  a  detection  threshold  that  achieves  a  probability  of  false 
alarm  of  10~3,  and  then  using  5  x  104  trials  for  each  SNR^  to  determine  Pd-  For  conve¬ 
nience,  SNR*vg  are  assumed  to  be  equal  across  transmitters,  i.e.,  SNR!^  =  SNRavg  V  i. 
The  channel  coefficients  {/x*  :  i  —  1 . . .  Nt}  are  chosen  randomly  and  scaled  to  achieve 
the  desired  SNR!^  relative  to  the  fixed  noise  power  a2  =  1CT6.  The  transmit  sig¬ 
nals  {u*  :  i  =  1 ...  Nt}  are  also  chosen  randomly  according  to  u*  =  expjjfT},  where 
Gl  G  MLxl  is  a  random  phase  vector  with  i.i.d.  elements  uniformly  distributed  on 
[0,  27t] ,  such  that  ||iT||2  =  L.  Recall  from  the  discussion  in  Sec.  3.3.3  that  detection 
performance  is  independent  of  the  specific  values  of  /x*  and  u*,  and  depends  on  only 

2Note  that  SNR^vg  describes  the  average  SNR  following  surveillance  channel  formation.  This 
may  be  regarded  as  an  average  input  SNR  because  it  does  not  include  processing  gain.  In  contrast, 
detection  performance  is  often  quantified  in  active  radar  in  terms  of  the  output  SNR,  which  includes 
processing  gain.  Input  SNR  is  chosen  here  rather  than  output  SNR  because,  in  contrast  to  active 
radar  processing,  the  concept  of  processing  gain  is  not  well  defined  for  passive  processing  as  performed 
by  the  SS-GLRT.  In  particular,  the  SS-GLRT  statistic  cannot  be  divided  into  signal  and  noise  terms, 
thereby  preventing  calculation  of  their  power  ratio. 
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Figure  13.  Predicted  (Pred)  and  Simulated  (Sim)  Pd  Curves  for  a  PMR  Network  with 
a  Varying  Number  of  Transmitters  and  Receivers 


their  respective  energies. 

In  all  cases,  the  predicted  and  simulated  results  agree,  verifying  the  distributions 
presented  in  Sec.  3.3.  In  addition,  detection  sensitivity  appears  to  improve  mono- 
tonically  with  increasing  SNR4  ,  Nt,  and  Nr,  as  expected.  It  also  appears  that 
detection  sensitivity  is  improved  more  significantly  by  the  number  of  receivers  in  the 
network  than  by  the  number  of  transmitters.  Specifically,  detection  sensitivity  does 
not  depend  solely  on  the  number  of  bistatic  pairs;  rather,  for  a  given  total  num¬ 
ber  of  pairs,  detection  sensitivity  improves  with  an  increasing  number  of  receivers. 
For  example,  consider  the  results  for  (Nt,Nr)  =  (2,3)  and  (Nt,Nr)  =  (1,6).  Al¬ 
though  both  configurations  use  NtNr  =  6  bistatic  pairs,  (Nt,Nr)  =  (1,6)  exhibits 
superior  detection  sensitivity.  The  same  holds  when  comparing  (Nt,  Nr)  =  (4,  3)  and 
(Nt,Nr)  =  (2,6),  where  (Nt,Nr)  =  (2,6)  exhibits  superior  detection  sensitivity.  This 
result  contrasts  with  active  MIMO  radar  with  distributed  antennas,  where  detection 
sensitivity  depends  on  only  the  total  number  of  bistatic  pairs  in  the  system  [24] .  This 
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Figure  14.  Additional  Pd  Curves  for  PMR  Networks  with  a  Varying  Number  of  Trans¬ 
mitters  and  Receivers 

asymmetry  between  transmitters  and  receivers  is  also  evident  when  comparing  net¬ 
works  with  the  same  number  of  independent  TDOAs/FDOAs.  This  is  illustrated  in 
Fig.  14,  which  depicts  additional  Pd  curves  for  different  cases  than  were  previously 
considered  in  Fig.  13.  Consider  the  results  for  (Nt,  Nr)  =  (4,  2)  and  (Nt,  Nr)  =  (1,  5). 
Although  both  configurations  realize  Nt(Nr  —  1)  =  4  independent  TDOAs/FDOAs, 
(Nt,  Nr)  =  (1,  5)  exhibits  slightly  superior  detection  sensitivity. 

Although  the  specific  reason  for  this  asymmetry  is  unclear,  it  might  be  explained 
by  noting  that  increasing  the  number  of  transmitters  introduces  more  nuisance  pa¬ 
rameters  into  the  problem  in  the  form  of  unknown  transmitted  signals.  Conversely, 
increasing  the  number  of  receivers  provides  more  observations  of  each  unknown  trans¬ 
mitted  signal,  which  enables  each  signal  to  be  (implicitly)  estimated  more  precisely 
during  calculation  of  the  SS-GLRT  statistic. 
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3. 5. 1.2  Received  Signal  Length. 


The  previous  experiment  for  (Nt,  Nr )  =  (2,  3)  was  repeated  as  L  was  varied  over 
L  —  [1,3,10,30,100,300,1000,3000,10000].  The  resulting  Pd  curves  for  L  —  1  to 
L  =  1000  are  shown  in  Fig.  15.  As  expected,  detection  sensitivity  improves  with 
increasing  L.  This  improvement  may  be  quantified  in  terms  of  an  integration  gain 
Gmt(L),  defined  as  the  separation  between  the  Pd  curve  for  a  given  length  L  >  1  and 
the  Pd  curve  for  L  =  1  at  Pd  =  0.90.  For  instance,  the  difference  between  the  L  —  1 
and  L  —  30  curves  at  P,i  =  0.90  is  approximately  10  dB,  indicating  Gmt  ~  10  dB. 
Fig.  16  shows  G-'int  versus  L.  As  shown,  Gint  varies  between  L0-7  at  L  —  10  (10  dB) 
to  just  below  L 0,6  at  L  =  10000  (40  dB).  This  behavior  is  inconsistent  with  coherent 
integration  of  L  samples  in  active  radar  processing,  for  which  Gint(L)  =  L.  Rather, 
it  is  suggestive  of  non-coherent  integration  of  L  samples,  for  which  G;nt  generally 
approaches  L05  with  large  L  [71].  This  result  is  likely  due  to  the  fact  that  each 
element  of  G*s  is  the  inner  product  of  two  noisy  surveillance  signals.  Calculation  of 
this  inner  product  results  in  signal-noise  cross  terms  that  are  similar  to  those  produced 
in  square-law  non-coherent  integration.  In  contrast,  coherent  integration  is  realized 
in  active  radar  processing  by  matched  filtering,  in  which  each  noisy  surveillance  signal 
is  correlated  with  a  noiseless  reference  signal. 

3.5.2  Ambiguity  Performance. 

This  section  illustrates  the  ambiguity  performance  of  the  SS-GLRT  for  a  scenario 
with  two  stationary  transmitters  at  d1  =  [0.5,4]  km  and  d2  =  [—0.5, —4]  km,  three 
stationary  receivers  at  r1  =  [—4,2]  km,  r2  =  [—4,0.5]  km  and  r3  =  [—4, —2.5]  km, 
and  a  target  at  t  =  [4,  0]  km.  Both  single-element  and  6-element  uniform  linear  array 
(ULA)  receivers  are  considered.  The  transmitters  have  carrier  frequencies  of  8.0  and 
8.1  GHz,  respectively,  and  isotropic  Pe*rp  =  50  W.  Complex  baseband  signals  sjf  are 
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Figure  15.  Pd  Curves  for  a  PMR  Scenario  with  2  Transmitters,  3  Receivers,  and 
L  =  [1,  3, 10, 30, 100,  300, 1000]  Samples  per  Signal 


Figure  16.  Integration  Gain  Gjnt  as  a  Function  of  L  for  a  PMR  Scenario  with  2 
Transmitters  and  3  Receivers 


simulated  according  to  (51)  (ignoring  receiver  noise)  with  sampling  rate  fs  =  500  kHz, 
length  T  —  10  ms,  and  u*  =  exp {jO1},  where  9‘  G  MLxl  is  a  random  phase  vector  with 
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i.i.d.  elements  uniformly  distributed  on  [0,  27r]  and  length  L  =  fsT  =  5000.  The  target 
has  an  isotropic  10  dBsm  RCS,  i.e. ,  a =  VTO  for  all  i  and  j.  Note  that  the  receivers 
are  not  phase-synchronized  due  to  the  random  phase  6 J  in  the  channel  coefficient  7^ 
according  to  (38)  in  Sec.  2.2.2. 

In  contrast  to  usual  discussions  of  radar  ambiguity,  which  examine  the  effect  of 
delay  and  Doppler  mismatch  on  the  matched  filter  output  assuming  a  given  transmit 
waveform,  this  section  examines  the  effect  of  position  and  velocity  mismatch  between 
the  hypothesized  state  (p,  p)  and  the  actual  target  state  (t,t)  on  the  SS-GLRT  test 
statistic  £ ss .  This  perspective  provides  insight  into  how  system  properties,  such  as  the 
number  of  receivers  and  transmitters  and  their  relative  geometries,  affect  system-level 
ambiguity  in  the  detection  domain,  which  in  this  formulation  is  Cartesian  position- 
velocity  space.  Nonetheless,  as  discussed  in  Sec.  3.4,  system-level  ambiguity  in  Carte¬ 
sian  position-velocity  space  depends  on  waveform  ambiguity  in  delay-Doppler  space 
in  that  the  elements  of  the  Gram  matrix  G*s  sample  the  AF  of  the  ith  transmit  wave¬ 
form  u*,  y*(Af,  An).  Fig.  17  depicts  the  normalized  AF  for  one  of  the  signals  in  the 
present  scenario  as  a  function  of  delay  At  =  A ijfs  and  Doppler  A fd  =  fsAv/{2ii) 
mismatch.  As  shown,  delay  resolution  is  l//s  =  2  /rs,  Doppler  resolution  is  1/T  =  100 
Hz,  and  the  AF  floor  is  101og10(L)  «  37  dB  below  the  mainlobe  peak. 

3. 5. 2.1  Stationary  Target. 

Consider,  first,  a  stationary  target.  Fig.  18  depicts  £ss  as  a  function  of  hypothesized 
position  p  =  [ px,Py ]  when  the  hypothesized  velocity  is  matched  to  the  true  target 
velocity  p  =  t  =  0,  i.e.,  £<.,,(p,  p)|p=0-  Each  receiver  consists  of  a  single  element, 
Ne  =  1,  with  isotropic  element  pattern  G^(x)  =  1  V  j.  As  shown,  the  test  statistic 
peaks  at  the  true  target  position,  and  it  exhibits  a  ridge  along  the  target  iso-TDOA 
contours  associated  with  each  of  the  three  receiver-receiver  pairs.  Each  iso-TDOA 
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Figure  17.  Example  Transmit  Signal  Ambiguity  Function  y*(Ar,  A  fd)  (dB) 

contour  represents  the  positions  at  which  the  TDOA  of  the  hypothesized  state  equals 
that  of  the  actual  target  state,  i.e.,  Al3p  =  Ai{k  for  the  jkth  receiver-receiver  pair. 
This  is  illustrated  in  Fig.  19,  which  depicts  the  TDOA  A r3pk  =  At?k  /  fs  as  a  function 
of  p  for  j  —  2  and  k  —  3.  The  target  iso-TDOA  contour  A rpa  =  A rf3  is  denoted 
by  the  red  dash-dotted  hyperbola.  Note  that  Ai/p]k  =  Ai/t  jk  =  0  for  all  i,  j,  k  and 
p  in  this  example,  i.e.,  the  hypothesized  and  actual  FDOAs  equal  zero  everywhere, 
because  the  considered  scenario  is  stationary. 

This  result  is  consistent  with  the  analysis  presented  in  Sec.  3.4.  From  (109),  the 
jkth  element  of  G'ss  is  proportional  to  y*(A £Jpk  —  A^fc,0).  When  APpk  =  Ad{k,  i.e., 
along  the  target  iso-TDOA  contour,  this  element  samples  the  AF  peak.  When  there 
is  mismatch  between  At3p  and  Ad{k,  this  element  samples  the  AF  at  a  location  that 
is  offset  from  the  AF  peak  along  the  zero-Doppler  (A fd  =  0)  cut.  This  results  in  the 
appearance  of  ambiguity  “ridges”  along  the  iso-TDOA  lines  that  correspond  to  the 
target  signal  TDOA. 

The  effect  of  multichannel  receivers  is  seen  in  Fig.  20,  which  depicts  £ss(p,  p)|p=0 
when  each  receiver  is  a  6-element  ULA  facing  the  +px  direction  with  1.875  cm  cl¬ 
ement  spacing  and  unity  element  patterns.  In  comparison  to  Fig.  18,  beamforming 
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Figure  18.  £ss(p,  p)|p=0  (dB)  for  t  =  0  and  Single-Element  Receivers 
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Figure  19.  TDOA  A r^k  =  A iJpk / fs  (/j.s)  for  j  =  2  and  k  =  3 


during  surveillance  channel  formation  sharpens  the  response  of  £ss  around  t.  This  is 
expected  from  the  discussion  of  the  angular  mismatch  loss  plpl  in  Sec.  3.4.  Thus,  use 
of  multichannel  receivers  appears  to  reduce  detection  ambiguity  (and  improve  target 
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Figure  20.  £ss(P’P)lp=o  (dB)  for  t  =  0  and  Six-Element  ULA  Receivers 


localization  performance)  in  comparison  to  single-element  receivers. 


3. 5. 2. 2  Moving  Target. 

Next,  consider  a  moving  target  with  velocity  t  =  \px,py\  =  [—75,  75]  m/s.  Fig.  21 
depicts  £,ss  as  a  function  of  p  when  p  =  t  assuming  isotropic  single-element  receivers. 
As  expected,  still  peaks  at  p  =  t,  where  the  iso-TDOA  and  iso-FDOA  constraints 
are  satisfied  for  all  receiver-receiver  pairs.  However,  in  comparison  to  Fig.  18,  Fig.  21 
exhibits  the  additional  effects  of  FDOA  mismatch  between  receiver  pairs  when  p  ^  t. 
This  is  seen  by  the  manner  in  which  decreases  with  separation  from  the  target 
iso-FDOA  contours,  which  are  depicted  by  dashed  blue  lines,  in  addition  to  the  iso- 
TDOA  contours.  Such  FDOA  mismatch  results  even  though  p  =  t  at  all  p  because 
the  hypothesized  FDOA  that  is  induced  for  each  receiver-receiver  pair  still  varies  as  a 
function  of  p.  This  is  illustrated  in  Fig.  22,  which  illustrates  the  hypothesized  FDOA 
(kHz)  as  a  function  of  p  for  one  receiver-receiver  pair,  assuming  p  =  t.  The  target 
iso-FDOA  contour  is  shown  by  the  dashed  blue  line.  Unlike  iso-TDOA  contours, 
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Figure  21.  £ss(p,  p)|^=0  (dB)  for  t  =  0  and  Single-Element  Receivers 
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Figure  22.  FDOA  ^^A^*jfc(p,  p)|^  ^  (kHz)  for  i  =  1,  j  =  2,  and  k  =  3 


which  are  hyperbolas  in  two  dimensions,  iso-FDOA  contours  do  not  admit  a  simple 
analytic  description. 

Finally,  Fig.  23  depicts  ^ss(p,p)|.  when  each  receiver  is  the  6-element  ULA 
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Figure  23.  £ss(P’P)lp=o  (dB)  for  t  =  0  and  Six-Element  ULA  Receivers 


described  previously.  Similar  to  Fig.  20,  the  peak  around  the  true  target  position  is 
sharpened  due  to  the  attenuation  resulting  from  angular  mismatch  between  the  actual 
and  hypothesized  target  positions.  Thus,  the  statistic  £ss  is  influenced  by  an  angular 
constraint  with  multichannel  receivers  in  addition  to  the  iso-TDOA  and  iso-FDOA 
constraints. 


3.5.3  Numerical  Challenges. 

Calculating  the  distributions  required  for  Fig.  13  proved  challenging.  Direct  cal¬ 
culation  of  the  PDF  fi(£',Nr,L,Qs)  hi  (88)  proved  infeasible  for  large  ( and  order 
L  >  10  without  variable  precision  arithmetic  because  the  required  calculations  ex¬ 
ceeded  the  largest  double  precision  floating-point  number  defined  by  IEEE  Standard 
754,  which  is  approximately  1.79  7  69  x  10308.  Use  of  variable  precision  arithmetic 
in  Mathematica  alleviated  this  problem;  however,  the  resulting  calculations  were  ex¬ 
tremely  time  consuming.  It  was  found  that  estimating  this  PDF  by  numerically 
differentiating  the  associated  cumulative  distribution  function  (CDF),  which  is  also 
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given  by  [48],  provided  a  significant  speed  improvement.  This  approach  avoids  calcu¬ 
lation  of  the  matrix  inverse  \Er~1(0),  which  is  a  challenging  operation  because  \f,(0) 
is  often  poorly  conditioned  due  to  the  large  dynamic  range  of  its  entries.  The  PDF 
/i(£;  Nr,  L ,  Qs)  was  successfully  calculating  using  this  approach  over  the  range  of  con¬ 
sidered  SNRavg  for  L  <  100  when  Nr  <  6.  Calculations  for  larger  L  or  Nr  did  not 
result  in  valid  CDFs  and  remains  a  topic  for  future  investigation. 

3.6  Conclusion 

In  this  chapter,  a  GLRT  for  centralized  detection  in  PMR  networks  that  oper¬ 
ate  without  direct-path  reference  signals,  termed  the  surveillance-surveillance  GLRT, 
was  presented.  The  distributions  of  the  SS-GLRT  test  statistic  were  identified  un¬ 
der  both  hypotheses  using  recent  results  from  random  matrix  theory  that  pertain  to 
the  distributions  of  the  largest  eigenvalues  of  complex  Wishart  matrices.  This  detec¬ 
tion  problem  was  shown  to  be  equivalent  to  detection  in  PSL  sensor  networks  when 
the  source  simultaneously  emits  on  multiple  transmit  channels.  Through  numerical 
simulation,  it  was  shown  that  transmitters  and  receivers  contribute  asymmetrically 
to  detection  sensitivity  in  such  sensor  networks.  It  was  also  shown  that  increasing 
the  received  signal  length  improves  sensitivity  in  a  manner  that  is  consistent  with 
non-coherent  integration  gain  rather  than  coherent  integration  gain.  Analysis  and 
simulation  of  system-level  ambiguity  also  showed  that  the  SS-GLRT  localizes  targets 
in  angle,  TDOA,  and  FDOA,  in  a  manner  analogous  to  PSL  sensor  networks  that 
utilize  AOA,  TDOA,  and  FDOA. 
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IV.  Passive  MIMO  Radar  Detection  with  References 


This  chapter1  addresses  the  problem  of  target  detection  in  PMR  networks  when 
direct-path  signals  are  available,  shown  in  Fig.  24.  As  discussed  in  Sec.  1.4.1,  the 
conventional  approach  to  PMR  detection  approximates  the  matched  filter  processing 
used  in  AMR  networks.  In  particular,  for  each  bistatic  (transmitter-receiver)  pair, 
the  direct-path  and  target-path  signals  are  isolated  into  reference  and  surveillance 
channels,  respectively,  and  the  reference-surveillance  cross-ambiguity  function  (CAF) 
is  calculated.  Calculation  of  the  CAF  is  analogous  to  matched  filtering  in  active  radar, 
except  the  (noisy)  reference  channel  is  used  in  place  of  a  known  transmit  signal.  As 
discussed  in  Sec.  1.4.1,  this  approach  is  ad  hoc  for  PMR  because  matched  filtering 
is  only  optimal  for  known  transmit  signals.  In  low-DNR  scenarios,  this  approach 
is  severely  degraded  due  to  mismatch  between  reference  channels  and  the  originally 
transmitted  signals.  Consequently,  this  approach  does  not  represent  a  general  solution 
to  the  PMR  detection  with  references  problem,  but  is  only  suitable  for  high-DNR 
scenarios,  in  which  references  closely  approximate  the  unknown  transmit  signals. 

This  chapter  presents  a  novel  centralized  GLRT  for  the  PMR  with  references 
lrThe  material  in  this  chapter  appears  in  [37]. 


Figure  24.  A  Passive  MIMO  Radar  Sensor  Network 
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detection  problem.  This  detection  problem  is  formulated  in  terms  of  the  incident 
signal  at  each  element  of  a  passive  radar  array  receiver,  which  is  the  superposition 
of  a  target-path  signal,  a  direct-path  signal,  and  receiver  noise.  In  addition,  the  de¬ 
tection  problem  formulation  assumes  non-isotropic  target  scattering,  and  makes  no 
assumption  about  coherence  between  receivers.  The  resulting  detector,  termed  the 
reference-surveillance  GLRT  (RS-GLRT)  because  it  forms  and  processes  both  refer¬ 
ence  and  surveillance  channels,  represents  a  general  solution  to  the  PMR  detection 
problem  that  addresses  how  detection  sensitivity  depends  on  both  SNR  and  DNR.  In 
this  way,  the  RS-GLRT  confirms  the  research  hypothesis  presented  in  Sec.  1.5. 

The  detection  and  ambiguity  performance  of  the  RS-GLRT  is  compared  against 
that  of  two  closely  related  detectors,  termed  the  matched  filter  GLRT  (MF-GLRT) 
and  surveillance-surveillance  GLRT  (SS-GRLT),  which  are  GLRTs  for  centralized 
detection  in  AMR  and  PSL  sensor  networks,  respectively.  PSL  and  AMR  represent 
two  extremes  in  terms  of  knowledge  about  the  transmitted  signals,  i.e. ,  the  signals 
are  entirely  unknown  in  PSL  and  entirely  known  in  AMR.  It  is  shown  that  the  RS- 
GLRT  varies  between  both  extremes  as  a  function  of  the  average  DNR.  When  the 
DNR  is  low,  PMR  detection  performance  approximates  PSL  detection  performance. 
When  the  DNR  is  high,  PMR  detection  performance  approaches  AMR  detection 
performance.  It  is  also  shown  that  the  RS-GLRT  test  statistic,  which  is  a  function  of 
the  largest  eigenvalues  of  complex  Wishart  matrices,  is  approximately  proportional 
to  the  AMR  test  statistic  under  the  high-DNR  condition.  These  results  lay  the 
groundwork  for  the  unified  detection  framework  presented  in  Chapter  V. 

4.1  Signal  Model 

Consider  a  PMR  network  with  Nt  transmitters  and  Nr  receivers.  Fig.  25  de¬ 
picts  the  geometry  and  signal  environment  of  the  ?'jth  bistatic  pair,  consisting  of  the 
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Figure  25.  The  ij th  Bistatic  (Transmitter-Receiver)  Pair  in  a  PMR  Network 

ith  transmitter  and  jth  receiver.  The  position  and  velocity  of  the  ith  transmitter  are 
denoted  by  d*  and  d*,  and  the  position  and  velocity  of  the  jth  receiver  by  r-7  and 
r-7.  The  spatial  states  of  all  transmitters  and  receivers  are  assumed  known.  In  the 
passive  radar  signal  environment,  direct-path,  target-path,  and  clutter-path  signals 
are  incident  at  each  array.  It  will  be  shown  that  the  RS-GLRT  can  be  interpreted 
as  forming  surveillance  and  reference  channels,  and  clutter  can  be  mitigated  in  these 
channels  through  a  variety  of  techniques  discussed  in  Sec.  1.1.2.  Furthermore,  detec¬ 
tion  performance  in  a  clutter-free  noise-only  environment  represents  an  upper  bound 
against  which  the  performance  of  clutter-suppressing  detectors  can  be  compared.  As 
such,  only  direct-path  and  target-path  signals  are  considered  here. 

Assume  the  transmitted  signals  are  separable  in  frequency,  and  that  each  receiver 
is  equipped  with  an  A^-element  array.  The  incident  signal  at  the  nth  element  of 
the  jth  receiver  is  the  sum  of  the  direct-path  and  target-path  signals  from  each 
transmitter.  This  signal  is  channelized  in  frequency,  demodulated  to  baseband,  and 
sampled  in  time.  Chapter  II  describes  this  process  in  detail.  Let  slj  e  CL*xl  denote 
the  lcngth-L*  sampled  complex  baseband  signal  on  the  zth  frequency  channel  of  the 
nth  array  element  of  the  jth  receiver  array.  For  simplicity,  it  is  assumed  that  Ll  = 
L  V  i.  Assuming  a  target  at  position  t  and  velocity  t,  s b  is  given  by  Eq.  (50)  in 
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Sec.  2.3  as 


s*J  = 


V 


'J  ,.i d'  ■:  plj 


u 


+  7’VC(t) 


u'  +  n 


y 


(117) 


where  u*  G  CLxl  is  the  length- L  discrete-time  complex  baseband  signal  emitted  by 
the  ith  transmitter;  T>lJ  e  CLxL  and  e  CLxi  are  unitary  linear  operators  that 
account  for  the  delays  and  Doppler  shifts  imparted  to  the  ith  transmit  signal  as  it 
propagates  to  the  jth  receiver  along  the  direct-path  and  target-path  channels,  respec¬ 
tively;  'd^(d')  and  'd^(t)  are  differential  phases  at  the  nth  array  element  with  respect 
to  a  reference  element  due  to  plane  wave  propagation  from  the  ith  transmitter  and 
target  directions,  respectively;  7 ^  and  7 \3  account  for  the  complex  scaling  of  iT  as¬ 
sociated  with  the  ijth  direct-path  and  target-path  channels,  respectively;  njf  e  Cixl 
is  circular  Gaussian  noise  distributed  as  CN{Ol,  °‘2Il)  with  known  variance  cr2;  and 
j  =  \/—  1 .  Here,  0^  denotes  the  length- L  zero  vector,  and  In  is  the  L  x  L  identity 
matrix.  The  delay-Doppler  operators  T>lrj  and  T>\3  are  defined 


(ns) 

u«  =  ©(7,7)  (ns) 


where  T>(£,  v)  is  the  delay-by-£  and  Doppler  shift-by-u  operator  defined  by  Eq.  (47) 
in  Sec.  2.3,  fj  and  i]3  are  the  delays  of  the  direct-path  and  target-path  channels  in 
samples,  respectively,  and  z/)2  and  u2t3  are  the  Doppler  shifts  of  the  direct-path  and 
target-path  channels  in  radians  per  sample,  respectively.  Noise  is  assumed  to  be  in¬ 
dependent  across  transmit  bands,  receivers,  and  array  elements,  i.e. ,  E{n^(n^)^}  = 
<j26n-m6i-k8j-ilL,  where  (-)H  is  the  Hermitian  transpose  and  Sx  is  the  Kronecker  delta. 
The  transmit  signal  u*  is  defined  such  that  ||u*||2  =  L.  Additional  detail  on  the  terms 
in  (117)  is  found  in  Chapter  II. 
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Let  av(x)  E  CNeXl  denote  the  spatial  steering  vector  in  the  direction  of  x, 

aij(x)=  e^(x)  e^'(x)  •••  e^(x)  T  (120) 

For  convenience,  let  acj  =  alJ"(d*)  and  a(J  =  ay  (t).  Then,  the  concatenation  of  the 
time  series  vectors  from  all  Ne  elements  of  the  jth  receiver  in  the  ith  frequency  band, 
=  [(sf)T  •  •  •  (s^)r]r  E  CNeLxl ,  can  be  written  as 

s"  =  (M;/  +  M;r)u'  +  nij  (121) 

where  iTJ  =  [(iT/)T  •  •  •  (n^JT]T  E  CNeLxl,  the  matrices  Mj?  and  are  defined  as 

Mj‘  =  Ydj  (aj‘  ®  V*)  E  CN*LxL  (122) 

M\j  =  ^  (af  ®  e  C^Lxi  (123) 

and  <E)  is  the  Kronecker  product. 

4.2  Detectors 

This  section  presents  a  derivation  of  the  RS-GLRT  detector.  Following  this,  the 
SS-GLRT  and  MF-GLRT  detectors,  which  correspond  to  detection  in  PSL  and  AMR 
sensor  networks,  respectively,  are  presented. 
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4.2.1  Reference-Surveillance  GLRT. 

For  a  given  position- velocity  cell  under  test  (p,  p),  the  detection  problem  may  be 
formulated  as  a  binary  hypothesis  test  between  Hi  and  Ho  hypotheses  as 

Hi  :  =  (Mj  +  Mj')ui  +  nij  (124) 

H0  :  s V  =  Mju*  +  nlj  (125) 

for  i  =  1 . . .  Nt  and  j  —  1 ...  Nr.  In  (124),  is  defined 

M”  =  (a*)  <g>  Vf)  e  CJVeLxL  (126) 

where  y*-7  is  the  target-path  channel  coefficient  associated  with  p,  a^  =  au(p),  and 
is  the  delay-Doppler  operator  associated  with  the  hypothesized  state  (p,  p).  Let  s* 
denote  the  concatenation  of  all  receiver  measurements  associated  with  the  ith  trans¬ 
mitter,  and  let  s  denote  the  concatenation  of  all  s*  across  transmitters, 


. ,  (siNr)T]T  e  cNrLxl 

(127) 

[(sT,.. 

. ,  (s^)rf  e  C™xl 

(128) 

Similarly,  let  Y(dp)  denote  the  vector  of  channel  coefficients  associated  with  the 
ith  transmitter,  where  the  subscript  notation  (-)(d,p)  denotes  either  (-)d  or  (-)p,  and 
let  7(rfiP)  denote  the  concatenation  of  all  l\dp)  across  transmitters2, 

Hv)  =  Kh)  ■■■dX)]T  ecNrX1  <129> 

In P)  =  [(7«rl)T  ■  ■  ■  (7gp))T]T  6  CAW'>"  (130) 

2Note  that  vectors  are  in  bold  (e.g.,  7 ld  and  yd)  while  scalars  are  not  (e.g.,  7^).  This  distinction 
can  be  difficult  to  see  in  print;  in  such  cases,  note  that  symbols  are  also  distinguished  by  their 
sub/superscripts. 
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Finally,  let  u  =  [(u1)7^ . . .  (uiVt)T]T  e  CNtLxl. 

Due  to  the  independence  of  the  receiver  noise  across  transmitter  channels,  the 
conditional  PDF  of  s  under  Hi,  pi(s  |  7d,7p,u),  factors  according  to 


Nt 

pi  (s  i  id,  iP, «)  =  n  &  i  yd,YP, ui) 

i= 1 


(131) 


where 


P\Y  I  Yd,  YP,  u*)  =  ci  exp<^  -  — 


Mf'u 


(132) 


and  =  M*/  +  M^.  The  PDF  of  s  under  Ho,  p0(s  I  7d,  u),  is  defined  similarly. 

In  this  formulation,  the  transmit  signals  u  and  channel  coefficients  'yd  and  ~{p 
are  considered  deterministic  and  unknown.  Thus,  hypothesis  Hi  is  composite  be¬ 
cause  the  PDF  of  s  in  (131)  is  parameterized  by  yd,  7 p,  and  u.  Hypothesis  Ho 
is  also  composite,  parameterized  by  7d  and  u.  Therefore,  the  GLRT  is  derived, 
which  replaces  these  unknowns  with  their  MLEs  in  the  likelihood  ratio  test  [49]. 
Let  h(  7d,  7  p,  u  |  s)  =  log  pi  (s  |  7d,  7  P,  u)  and  /0(  jd,  u  |  s)  =  logp0(s  |  7d,  u)  denote  the 
log-likelihood  functions  under  Hi  and  Ho ■  Then,  the  GLRT  may  be  written  as 


Hi 

max  Z1(7d,  7  ,  u  |  s)  -  max  l0( 7d,  u  |  s)  ^  Krs  (133) 

{■Td-Tp-u}  {7d,u}  Ho 


Consider  the  log-likelihood  R( 7d,7p,u  |  s).  From  (131), 


Nt 

h  (7 d,  Ip,  u  |  s)  =  l\  (Yd,  Yp,  u*  |  s’)  (134) 


where,  from  (132)  and  ignoring  an  additive  constant, 


-1  r 

I  S')  =  E  lls‘J  -  M?«,f  (135) 

3  =  1 
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It  is  shown  in  Appendix  B  that 


M*V 


i  II 2 


t*>'\ 


+  §: 


*? 


Hr  iT 


I2  +  Ely 


(136) 


where  s*J  =  (Vlcl)HslJ  and  slJ  =  (V^)Hs1^  are  delay-Dopplcr  compensated  reference 
and  surveillance  signals  defined  by  (241)  and  (243),  Hr  and  Hls  are  complex  scalars 
defined  in  (246)  and  (247)  that  account  for  the  composite  scaling  of  the  reference 
and  surveillance  channels,  and  El^±  represents  the  energy  of  not  captured  by 
and  slsj,  defined  in  (230).  Note  that  Hr  an(l  Hs  are  defined  in  terms  of  y(/  and  y*-7, 
respectively.  Thus,  they  may  also  be  regarded  as  deterministic  unknowns,  and  they 
replace  and  y*J  in  the  following  development.  Accordingly,  let  H\rs)  and  H(r,s)  t>e 
defined  as 


M(r,s) 


[A*)---dfdTeC"'xl 


(137) 

(138) 


where  (-)(r,s)  denotes  either  (•),.  or  (-)s.  Using  (136)  in  (135), 


Nr 

E 

1  =  1 


t*>'\ 


+  s: 


V 


Hr**1 


I2  +  Ely 


(139) 


The  MLE  of  H(rs)  *s  given  from  (139)  by 


u 


)Hsij 

>  h(r,S) 


U' 


(140) 


Substituting  (140)  into  (139),  and  simplifying,  gives 


u^dh'phf'u'A 
IM|2  ) 


(141) 
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where  <3^  =  [<3>*  <3>*  ] ,  the  matrices  <3>*  and  <h*  are  defined  as 


®\r,s)  = 


il 

(r,s)  ? 


~iNr 

5h,s)J 


iLxNr 


(142) 


the  scalar  El  —  JT  1 1  s13  \ | 2  denotes  the  cumulative  energy  of  the  measurements  associ¬ 
ated  with  the  ith  transmitter,  and 


list'll2  =  Ill'll2  +  Ill'll2  +  E^X  (143) 

Let  Ai(-)  denote  the  largest  eigenvalue  of  its  matrix  argument,  and  let  v^-)  denote  the 
associated  eigenvector.  Then,  the  Rayleigh  quotient  in  (141)  achieves  its  maximum 
value,  Xi  (<h,1<h,1H) ,  when  iT  =  Vi  (T) <!>',")  [44,  p.  176].  Therefore,  u*  = 
and  (141)  becomes 

i\(K,  Ai,a‘|s‘)  =  ~li  (E‘  -  (144) 

Noting  that  Ai  (<3>,1<3>ff)  =  Ai($*1H$*1),  and  that  typically  2 Nr  <C  L,  it  is  more  efficient 
to  consider  the  Gram  matrix  G]  =  (<h*1)H<h^  e  ^2Nrx2Nr ,  giving 

l\(K, K,  U*  |  s')  =  -T  (e'  -  A^Gi))  (145) 

Therefore,  using  (134)  and  (145), 

1  Nt 

i,  (£,,£>  I  s)  =  --2J2  (E‘  -  MG’,))  (146) 

i= 1 

By  a  similar  procedure,  it  can  be  shown  under  "Ho  that 

1  Nt 

WAr,  U  |  S)  =  -  —  (fi‘  -  Aj  (G*r)) 

i= 1 
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(147) 


where  G*r  =  (<h*)i/<h*  G  CNrXNr.  Using  (146)  and  (147),  the  RS-GLRT  results  from 
(133)  as 

=  -^yRA,(Gi)  -  A, (Gy)  i  Kr,  (148) 

a  7  no 

4.2.2  Surveillance-Surveillance  GLRT. 

Alternatively,  a  PMR  system  may  use  only  the  target-path  signals  for  target 
detection.  The  resulting  detector  is  a  special  case  of  the  RS-GLRT  in  (148)  when 
only  surveillance  channels  are  formed  at  each  receiver,  i.e. ,  s*J  =  0  for  all  %  and  j, 

1  ^  Hi 

&»  =  -EAi(G«)  £  «**  <149) 

where  Glss  =  (4>*)if<h* .  This  is  the  SS-GLRT  detector  discussed  previously  in  Chapter 
III,  where  it  was  shown  that  PMR  detection  in  the  absence  of  direct-path  signals  is 
equivalent  to  PSL  detection,  as  the  target  in  PMR  can  be  interpreted  as  “emitting” 
the  scattered  target-path  signals  [40]. 

4.2.3  Matched  Filter  GLRT. 

For  comparison,  the  detector  that  results  when  the  transmit  signals  u  are  assumed 
known,  termed  the  matched  filter  GLRT  (MF-GLRT),  is  also  considered.  It  can  be 
derived  by  following  the  RS-GLRT  derivation  through  (141)  under  "Hi,  and  performing 
the  analogous  steps  under  "Ho,  resulting  in 

1  Nt  Nr  Hi 

tmf  =  |(U*)H§?|2  ^  Kmf  (150) 

Assuming  u  is  known  represents  an  AMR  scenario.  Note  that  Cmf  is  formed  by  delay- 
Doppler  compensating  the  surveillance  signals,  matched  filtering  each  bistatic  pair, 
and  non-coherently  integrating.  This  structure  is  common  in  AMR  detection  [24,  43]. 
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4.3  Distributions 


This  section  presents  PDFs  for  and  £,rs  under  Hi  and  "Ho3-  The  development 
considers  the  presence  or  absence  of  a  target  within  the  cell  under  test  (p,  p).  Under 
this  condition,  the  surveillance  signal  sz/  in  (243)  reduces  to 

s?  =  &i/#u*  +  n*J  (151) 

where  is  distributed  as  CM(Ol,  u2Il),  b\  —  1  under  Hi,  and  b\  —  0  under  Ho- 
Consequently,  slJ  is  distributed  as 

slj  ~  CM{bi  flu*,  a2IL)  (152) 

Similarly,  slJ  ~  CM  ul ,  (J2Il)  under  both  hypotheses. 

4.3.1  Matched  Filter  GLRT  Distributions. 

First,  consider  the  MF-GLRT  statistic  £mf  in  (150).  From  (152),  the  inner  product 
(u*)Hs*J  is  distributed  as  CM(biHljL,a2L),  and  the  statistic  can  be  written  as 

W  =  ^2^  ^(2ATtArr),Cm/  (153) 

where  xfk)  (  *s  a  non-central  chi-squared  random  variable  with  k  degrees  of  freedom 
and  non-centrality  parameter  (,  and 

o  Nt  Nr 

Cmf=ML  EE  |6iM?i|2  (154) 

i= 1  3= 1 
Nt 

=  2  biLNr  E  SNRlg  (155) 

i= 1 

3Note  that  PDFs  for  £ss  have  already  been  presented  under  Hi  and  Ho  in  Sec.  3.3. 
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where  SNR^vg  =  || /x*  \\2/(Nra2)  is  the  average  input  target-path  SNR  associated  with 
the  ith  transmitter  after  surveillance  channel  formation.  Consequently,  the  PDF  of 
Cm.f  i  denoted  pmf(£),  is  given  by 

Pmf(0  =  (£)  fX*  (f ;  2 NtNr,  Cm/)  (156) 

where  fx 2(x;  k,()  is  the  PDF  of  xfk)x,  he., 

,  \  (fc— 2)/4  _ 

fx2(x;k,C)  = -e~{x+^/2  I^J  h/2-i(VCx)  (157) 

and  Iv(z)  is  the  modified  Bessel  function  of  the  first  kind  of  order  v.  Note  that 
Cm/  =  0  under  "Ho,  and  (157)  reduces  to  the  central  chi-squared  distribution  with  k 
degrees  of  freedom. 

4.3.2  Reference-Surveillance  GLRT  Distributions. 

Next,  consider  the  RS-GLRT  statistic  £rs  in  (148).  The  Gram  matrix  G*r  follows 
the  non-central  uncorrelated  complex  Wishart  distribution  under  both  "Hi  and  "H0- 
In  particular,  Glrr  ~  WNr(L,  £*r,  f2*r),  where  cr2lNr  and  f Vrr  is  a  rank-1  non¬ 

centrality  matrix  with  non-zero  eigenvalue  Q.r  given  by 

Crr  =  =  LNr  DNR;vg  (158) 

and  DNRCvg  =  ||/pr||2/(fVrcr2)  is  the  average  input  direct-path  SNR  associated  with 
the  ith  transmitter  after  reference  channel  formation.  This  distribution  is  unchanged 
between  hypotheses  because  the  Gram  matrix  G*r  depends  only  on  the  direct-path 
signals.  The  Gram  matrix  G)  also  follows  the  non-central  uncorrelated  complex 
Wishart  distribution  under  both  hypotheses,  i.e. ,  G*  ~  w2Nr(L,x\,n[),  where  = 
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a2l2Nr  and  fij  is  a  rank-1  non-centrality  matrix  with  a  non-zero  eigenvalue  C{  that 
varies  between  hypotheses  as 

{LNr  (DNR’  +  SNR*  )  under  U  i 

(159) 

LNr  DNR4vg  under  H0 

Therefore,  the  individual  PDFs  of  Ai(G*rr)  and  Ai(Gj)  may  be  found  using  the  result 
in  [48].  However,  Ai(Gj)  and  Ai(G*rr)  are  not  independent  because  G*r  is  a  diag¬ 
onal  block  of  Gj .  The  joint  distribution  of  the  eigenvalues  of  a  complex  Wishart 
matrix  and  the  eigenvalues  of  its  diagonal  blocks  does  not  appear  to  exist,  precluding 
determination  of  the  exact  distribution  of  £rs  for  now. 

Nonetheless,  it  is  now  shown  that  £,rs  approximates  £m/  within  a  scalar  constant 
assuming  high-DNR,  i.e.,  DNR^g  >  1  Vi,  and  assuming  a  high  average  power  ratio 
Jf  between  the  direct-path  and  target-path  signals,  i.e.,  p*  =  ||^,r||2/||/i.*  ||2  >>  1  V  i. 
Let  An(-)  denote  the  nth  largest  eigenvalue  of  its  matrix  argument.  Linder  high-DNR, 
Gj  is  approximately  rank-1,  i.e.,  Ai(Gj)/Aj(Gj)  3>  1  for  2  <  j  <  2Nr.  Consequently, 

Ai(Gj)  «  \J A2(Gj)  H - 1-  A^^Gj)  =  ||Gj||F  (160) 


where  ||-||F  denotes  the  Frobenius  norm  of  its  matrix  argument.  Note  that  Gj  can  be 
partitioned  into  four  blocks, 


S  S 

s  r 

_A 

I 

0 

Co  «*• 

Co 

Gjr 

<biH  <f)i 

r  s 

i 

•<s>  S- 

o 

Cii 

VJTrs 

Qi 

(161) 
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where  (G*rJ  =  Gjr.  Accordingly, 
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Let  Fis=  1 1 G*  J 


Fi 


I|g;.i 

|F,  and  Flrr  =  |  G/||F.  Then,  letting  Flrs 

=  pi  /  pi 

rs/  rr’ 

l|G)| 

/  pi2  _L_  2  Pi2 

1  _  Tpi  /  i  |  ss  '  rs 

If  “  *rr\  1  +  7^2 

V  1  rr 

(163) 

~F^l  +  2Ff1 

(164) 

=  (l  +  F'r,  +  o(Kl)) 

(165) 

«  p(i + dt) 

(166) 

The  approximation  in  (164)  results  because  2 Ffs  3>  Ffs  under  the  high-p*  condition. 
Eq.  165  results  from  applying  the  Taylor  series  expansion  for  y/l  +  x  around  x  —  0, 
a/1  +  x  =  1  Px/2  +  0(x2),  which  converges  for  |x|  <  1.  Note  that  2 Ffs  <C  Ffr  under 
the  high-p*  condition,  which  justifies  the  approximation  in  (166). 


Similarly,  the  matrix  G*r  is  also  approximately  rank-1  under  the  high-DNR  con¬ 


dition.  Consequently, 


MG/.) 


~  pi 
rr 


(167) 


Therefore,  using  (166)  and  (167), 
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Furthermore,  ignoring  the  noise  in  s/  due  to  the  high-DNR  condition,  i.e. ,  s/  «  p/u*, 

•  2 

Fls  may  be  expanded  as 
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(169) 


Similarly, 
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(171) 


Using  (170)  and  (171)  in  (168), 
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Consequently,  £rs  can  be  approximated  as 

Z™  ~  (2)  ^(2JVtJVr),Cm/  (174) 

and  its  distribution  is  prs(£)  =  2/x2(2£;  2NtNr,  Cm/)-  Since  detection  performance  is 
unaffected  by  scaling  of  the  test  statistic  (provided  the  detection  threshold  is  similarly 
scaled),  the  RS-GLRT  detector  is  equivalent  to  the  MF-GLRT  detector  under  the 
high-DNR  and  high-jo*  conditions. 

The  validity  of  (174)  is  verified  in  Fig.  26,  which  compares  normalized  histograms 
generated  from  105  simulated  realizations  of  £rs  and  the  associated  predicted  PDFs 
according  to  (174)  under  both  hypotheses.  The  simulation  involves  a  scenario  with 
2  transmitters,  3  receivers,  and  L  =  1000  samples  per  transmit  signal.  The  channel 
scale  factors  and  were  selected  randomly  and  scaled  to  achieve  SNR^vg  =  —20 
dB  and  DNR,?lvg  =  +20  dB  for  all  i,  assuming  0 2  =  10~6,  such  that  jo*  =  +40  dB.  As 
shown,  the  predicted  PDFs  closely  match  the  empirical  results. 


4.3.3  Dependence  on  SNR  and  DNR. 

It  is  worth  noting  that  all  the  presented  PDFs  depend  only  on  SNR*ivg,  DNR*avg, 
Nt,  and  Nr.  This  is  seen  by  examining  the  non-centrality  parameter  Qmf  in  (155)  and 
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(a) 


(b) 


Figure  26.  Empirical  and  Predicted  PDFs  for  £rs  Under  (a)  "Ho  and  (b)  Tii 

the  complex  Wishart  non-centrality  matrix  eigenvalues  Qs,  Qr,  and  in  (87),  (158), 
and  (159),  respectively,  all  of  which  are  functions  of  LNrSNRlavg  and/or  LNrDNRavg. 
The  dependence  on  Nt  is  seen  in  the  non-centrality  parameter  (mf  in  (155)  and  by 
the  convolution  over  Nt  PDFs  in  (83).  Consequently,  detection  performance  does 
not  depend  on  the  specific  structure  of  the  transmit  waveforms  or  on  the  specific 
values  of  the  target-path  and  direct-path  SNRs  across  receivers.  Rather,  it  depends 
only  on  the  transmit  waveform  energy  ||u*||2  =  L  and  the  average  target-path  and 
direct-path  SNRs  across  receivers.  Note  that  the  LNr  factor  may  be  interpreted  as 
signal  processing  gain  resulting  from  coherent  integration  in  time  and  across  receiver 
channels.  In  addition,  detection  performance  is  expected  to  increase  monotonically 
in  LNr  SNR*vg.  This  is  demonstrated  by  the  empirical  results  presented  in  Sec.  4.5. 
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4.4  Interpretations 


This  section  discusses  the  relationship  between  the  statistics  £rs,  £ss,  and  £m/.  In 
particular,  the  Gram  matrix  in  £rs  is  examined,  and  it  is  shown  that  calculation 
of  its  block  matrices  G\s  and  G*r  entails  operations  that  are  characteristic  of  PSL 
and  AMR,  respectively.  The  following  discussion  distinguishes  between  the  true  target 
state  (t,  t)  and  the  hypothesized  cell  under  test  (p,  p).  If  there  is  a  mismatch  between 
t  and  p,  then  the  target  angle-of- arrival  is  possibly  mismatched  to  the  hypothesized 
angle-of-arrival,  resulting  in  a  mismatch  loss.  Using  (231),  it  can  be  shown  that  the 
surveillance  signal  following  beamforming,  slJ ,  is  given  by 

s?  =  Wuf  +  n?  (175) 


where  =  Cpf^J't  and  £*■?  is  the  angular  mismatch  loss, 
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and  the  scalars  (3pt ,  [3(jt  and  fijp  are  defined  by 
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(176) 


(177) 


Note  that  <  1,  and  Q}  =  1  when  =  a(J,  i.e.,  when  p  and  t  are  at  the  same 
angle  with  respect  to  the  jth  receiver. 


4.4.1  Surveillance-Surveillance  Processing. 

Consider  G'ss,  the  upper-left  block  of  G)  in  (161).  The  elements  of  this  matrix 
consist  of  the  pairwise  inner  products  of  the  surveillance  signals  {s*-7  :  j  —  1 . . .  Nr}. 
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Specifically,  the  jkth  element  of  G*s  may  be  expressed  as 


[GLh  =  sf  sf  (178) 

=  afVj’Vfa?  (179) 

=  s  fv(e>,  f f )  sf  (180) 

-f)sf  (181) 

^Aeyjk  ±A4'jk 

=  y?(A<y‘,Ai4'!‘)  (182) 


where  x(jfe(A7,  Az/)  is  the  CAF  between  the  ith  surveillance  channels  of  the  jth  and 
fcth  receivers,  defined  as 

L—l 

At  (AC,  AO  =  V  (183) 

1=0 

Ad^k  is  the  hypothesized  TDOA  at  the  jth  and  hth  receivers,  and  Az/*jfc  is  the 
corresponding  hypothesized  FDOA.  Thus,  the  elements  of  G*s  may  be  interpreted  as 
samples  of  the  pairwise  surveillance-surveillance  CAFs.  This  is  termed  surveillance- 
surveillance  processing. 

The  elements  of  G*s  may  also  be  expressed  in  terms  of  the  ambiguity  function 
(AF)  of  the  ith  transmit  signal  u®.  Substituting  (175)  for  slJ  and  s*fc  in  (181),  and 
ignoring  receiver  noise,  [G*Jjfc  may  be  expressed  as 


=  sf  Z>(A/«\  Az/j*)  sf  (184) 

=  (f  Cf  u‘)*®(Af  Af  l“)  (/if  Bf  u')  (185) 

=  iiffifuiHVH(A4J\Av;Jk)V(Af/^AC/k)u‘  (186) 

=  fifi if  x‘(A/f -Af‘,  A^-Ai**)  (187) 
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where  x*(A£,  Au)  is  the  AF  of  u*,  defined  as 

(188) 

1=0 

Noting  that  x*(Af,  Au)  peaks  when  its  arguments  are  zero  leads  to  the  following 
iso-range  and  iso-Doppler  conditions: 


A&k  =  A( 

P  L 

A//';'7''  =  A  uj’jk 

P  L 


(189) 

(190) 


In  words,  [G*J  -fe  peaks  when  (a)  the  hypothesized  TDOA  equals  the  actual  target 
TDOA  and  (b)  the  hypothesized  FDOA  equals  the  actual  target  FDOA.  This  is  rep¬ 
resentative  of  PSL  sensor  networks,  which  localize  targets  along  contours  of  constant 
TDOA  and  FDOA  [16,  86], 

4.4.2  Reference-Surveillance  Processing. 

Next,  consider  G*r,  the  upper-right  block  of  G)  in  (161).  The  elements  of  this 
matrix  consist  of  the  pairwise  inner  products  of  the  surveillance  signals  {s*J  :  j  = 
1 . . .  Nr }  with  the  reference  signals  {s*-7  :  j  =  1 . . .  Nr}.  It  can  similarly  be  shown  that 
the  elements  of  G*r  can  be  expressed  as 


(191) 

sfV'jVfs? 

(192) 

,^)VH  (<jf,  v*)4k 

(193) 

sfvie1’ 

(194) 

(195) 
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where  A =  if  —  i1^  is  relative  delay,  defined  as  the  difference  between  the  hypoth¬ 
esized  bistatic  delay  if  and  the  direct-path  delay  i'f,  Az/’jfe  =  uf  —  z/f  is  relative 
Doppler,  defined  as  the  difference  between  the  hypothesized  bistatic  Doppler  vfj  and 
the  direct-path  Doppler  z/f,  and  x*’/fc(A£,  Az/)  is  the  CAF  between  the  ith  surveil¬ 
lance  channel  of  the  jth  receiver  and  the  ith  reference  channel  of  the  A;th  receiver, 
defined  as 

x*,f  (Af,  A„)  =  £  (196) 

1=0 

Thus,  the  elements  of  G*r  may  be  interpreted  as  samples  of  the  pairwise  reference- 
surveillance  CAFs.  This  is  termed  reference-surveillance  processing. 

Similarly,  [G*r]  ■k  may  also  be  expressed  (ignoring  noise)  as 


[Gy*  =  sf  X>(A$‘,  A,/?)  s?  (197) 

=  (x«8) 

=  -v?)  (199) 


Noting  that  x*(Af,  Az/)  peaks  when  its  arguments  are  zero  leads  to  the  following 
iso-range  and  iso-Doppler  conditions: 


if  (200) 

v?  (201) 

In  words,  [G*r]  -fc  peaks  when  (a)  the  hypothesized  bistatic  delay  equals  the  actual 
target  bistatic  delay  and  (b)  the  hypothesized  bistatic  Doppler  equals  the  actual  target 
bistatic  Doppler.  This  is  representative  of  the  processing  associated  with  each  bistatic 
pair  in  active  multistatic  and  MIMO  radar,  which  localize  targets  along  contours  of 
constant  bistatic  range  and  bistatic  Doppler  [64,  42], 
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4.5  Simulations 


This  section  compares  the  detection  and  ambiguity  performance  of  the  RS-GLRT, 
SS-GLRT,  and  MF-GLRT  detectors  via  numerical  simulation.  Discussion  of  the  sig¬ 
nificance  of  these  simulation  results  follows  in  Sec.  4.6. 

4.5.1  Detection  Performance. 

This  section  illustrates  how  the  detection  performance  of  each  GLRT  varies  with 
SNRavg,  DNRavg,  and  L.  Fig.  27  depicts  probability  of  detection  (Pd)  curves  as 
a  function  of  SNR^g  for  a  scenario  with  Nt  —  2  transmitters,  Nr  =  3  receivers,  and 
L  =  1000  samples  per  signal.  Each  curve  is  generated  empirically  using  105  trials  under 
"Ho  determine  the  detection  threshold  that  achieves  a  probability  of  false  alarm 
(Pfa)  of  10-3.  Pd  is  calculated  under  "Hi  using  104  trials  for  each  value  of  SNRavg.  A 
separate  RS-GLRT  curve  is  calculated  for  each  I)NR*  vg  as  it  is  varied  from  -40  to  +20 
dB  in  5  dB  increments.  For  convenience,  SNRavg  and  DNRavg  are  assumed  to  be  equal 
across  transmitters,  i.e.,  SNRavg  =  SNRavg  V  i  and  DNR*vg  =  DNRavg  V  i.  As  shown 
in  Fig.  27,  the  detection  sensitivity  of  the  RS-GLRT  improves  with  increasing  DNRavg. 
Specifically,  sensitivity  asymptotes  at  low  and  high  values  of  DNRavg  and  improves 
monotonically  with  increasing  DNRavg  at  intermediate  DNRavg.  This  asymptotic 
sensitivity  is  slightly  inferior  to  that  of  the  SS-GLRT  at  low  DNRavg  and  equal  to 
that  of  the  MF-GLRT  at  high  DNRavg. 

Another  view  of  these  results  is  shown  in  Fig.  28,  which  plots  the  RS-GLRT  Pd  as 
a  function  of  SNRavg  and  DNRavg.  Pd  =  0.5  and  Pd  =  0.9  iso-contours  are  depicted 
as  solid  black  lines.  SNRavg  values  at  which  the  SS-GLRT  and  MF-GLRT  achieve 
Pd  =  0.5  and  Pd  =  0.9  are  shown  as  vertical  dash-dotted  red  (SS-GLRT)  and  dashed 
blue  (RS-GLRT)  lines.  Lines  of  constant  p  are  also  depicted.  As  shown,  detection 
performance  may  be  divided  into  three  regions:  a  low-p  region;  a  high-DNR  region; 


80 


Figure  27.  Probability  of  Detection  Curves  for  Nt  =  2,  Nr  =  3,  and  L  =  1000 

and  a  transition  region.  The  low-p  border  is  defined  as  the  power  ratio  at  which  the 
SNRavg  of  the  RS-GLRT  Pd  =  0.9  iso-contour  decreases  one  1  dB  below  its  asymptotic 
value.  Similarly,  the  High-DNR  border  is  defined  as  the  value  of  DNRavg  at  which 
the  SNRavg  of  the  Pd  =  0.9  iso-contour  increases  1  dB  above  its  asymptotic  value. 

The  rationale  for  defining  the  region  boundaries  in  this  manner  is  seen  by  exam¬ 
ining  the  effect  of  L  on  detection  performance.  Fig.  29  depicts  Pd  =  0.9  iso-contours 
for  the  three  GLRTs  over  L  =  [1,  3, 10,  30, 100, 300, 1000,  3000, 10000].  As  shown,  the 
upper  boundary  occurs  at  approximately  the  same  DNRavg  for  all  L;  the  average 
for  L  >  10  at  DNRavg  =  1.21  dB  is  depicted.  Similarly,  the  lower  boundary  occurs 
at  approximately  the  same  p  for  all  L;  the  average  for  L  >  10  at  p  =  —4.51  dB 
is  depicted.  Thus,  the  upper  and  lower  region  boundaries  are  defined  in  terms  of  a 
constant  DNRavg  and  constant  p,  respectively.  Note  that  the  single-sample  ( L  =  1) 
performance  of  all  three  GLRTs  converges. 

Fig.  29  shows  that  detection  sensitivity  improves  with  increasing  L  for  every  detec- 
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Figure  28.  2D  View  of  RS-GLRT  P ^  for  Nt  =  2,  Nr  =  3,  and  L  =  1000 

tor,  however  the  rate  of  sensitivity  improvement  with  increasing  L  differs  between  de¬ 
tectors.  This  improvement  may  be  quantified  in  terms  of  an  integration  gain  Gint(L), 
defined  as  the  separation  between  the  Pd  curve  for  a  given  detector  and  length  L  >  1 
and  the  Pd  curve  for  that  same  detector  with  L  =  1  at  Pd  =  0.90.  For  instance,  the 
difference  between  the  L  =  1  and  L  =  100  curves  at  Pd  =  0.90  for  the  MF-GLRT 
detector  is  20  dB,  indicating  an  integration  gain  of  20  dB.  Fig.  30  shows  G-mt  versus 
L  for  the  MF-GLRT,  SS-GLRT,  and  the  RS-GLRT  in  both  the  high-DNR  and  low-p 
regions.  As  shown,  the  MF-GLRT  exhibits  an  integration  gain  of  G-mt(L)  =  L ,  which 
reflects  coherent  integration.  The  SS-GLRT  exhibits  Gmt  that  varies  between  L 0  7  at 
L  =  10  (10  dB)  to  just  below  L 0,6  at  L  =  10000  (40  dB).  This  reflects  non-coherent 
integration,  and  was  discussed  previously  in  Sec.  3. 5. 1.2.  Finally,  the  RS-GLRT  re¬ 
alizes  an  integration  gain  that  is  coherent  in  the  high-DNR  region  and  non-coherent 
in  the  low-p  region.  The  significance  of  this  result  is  discussed  in  Sec.  4.6. 
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Figure  29.  Depiction  of  RS-GLRT  Pd  =  0.9  iso-contours  for  Nt  =  2,  Nr  =  3,  and  L  =  [1,3,10,30,100,300,1000,3000,10000].  The 
SNRavg  values  at  which  the  SS-GLRT  and  MF-GLRT  achieve  Pd  =  0.9  are  depicted  by  dash-dotted  red  lines  and  dashed  blue 
lines,  respectively,  for  each  L. 


Figure  30.  G int  as  a  Function  of  L  for  a  PMR  Scenario  with  Nt  =  2  and  Nr  =  3 
4.5.2  Ambiguity  Performance. 

This  section  illustrates  the  ambiguity  performance  of  the  considered  GLRTs  by 
calculating  their  test  statistics  as  a  function  of  the  hypothesized  target  position 
P  =  [Px,Py]-  The  scenario  consists  of  two  stationary  transmitters  at  d1  =  [0.5,4]  km 
and  d2  =  [—0.5,  —4]  km,  three  stationary  receivers  at  r1  =  [—4,  2]  km,  r2  =  [—4,  0.5] 
km  and  r3  =  [—4,  —2.5]  km,  and  one  target  at  t  =  [4,  0]  km.  Stationary  and  mov¬ 
ing  target  scenarios  are  considered.  The  transmitters  have  carrier  frequencies  of  8.0 
and  8.1  GHz,  respectively,  and  isotropic  Pe*rp  =  50  W.  All  receivers  are  6-element 
uniform  linear  arrays  facing  the  +px  direction  with  1.875  cm  element  spacing  and 
unity  element  patterns,  i.e.,  G3e(-)  =  1.  Complex  baseband  signals  s]3  are  simulated 
according  to  Eq.  (117)  with  fs  =  500  kHz,  T  =  2  ms,  a2  =  2.0019  x  10~14  (-106.99 
dBm),  and  ul  =  exp{j0*},  where  9L  e  MLxl  is  a  random  phase  vector  with  i.i.d.  el¬ 
ements  uniformly  distributed  on  [0,  2i r],  and  L  =  fsT  =  1000.  The  target  has  an 


84 


isotropic  10  dBsm  RCS,  i.e.,  a y'  =  VTO  for  all  i  and  j.  Note  that  the  receivers  are  not 
phase-synchronized  clue  to  the  random  phase  6 i  in  the  direct-path  and  target-path 
channel  coefficients,  7^  and  7^,  according  to  Eqs.  (27)  and  (38)  in  Secs.  2.2.1  and 
2.2.2,  respectively.  Note  that  this  scenario  is  the  same  as  the  scenario  considered  in 
Sec.  3.5.2,  except  that  here  the  CPI  duration  T  is  2  ms  ( L  =  1000)  rather  than  10 
ms  (L  =  5000). 

Three  SNR-DNR  scenarios  are  considered,  given  by  (SNRavg,  DNRavg)  =  (—5,  —35) 
dB,  (—10,  —10)  dB,  and  (—15, 15)  dB.  These  SNR-DNR  scenarios  fall  within  the  low- 
71,  transition,  and  high-DNR  regions,  respectively,  and  they  are  identified  in  Fig.  28 
by  star  symbols.  Note  that  each  SNRavg  is  between  8  to  10  dB  above  the  SNRavg  that 
is  required  at  that  DNRavg  to  achieve  Pd  =  0.90.  High  SNRavg  values  are  chosen  so 
that  the  ambiguity  response  may  be  clearly  seen  above  noise.  In  order  to  achieve  each 
SNR-DNR  scenario,  it  is  necessary  to  scale  the  direct-path  and  target-path  signals 
that  result  from  simulation  of  the  scenario  described  in  the  previous  paragraph.  In 
particular,  simulation  of  this  scenario  according  to  the  signal  model  of  Sec.  4.1  results 
in  SNRavg  =  —43.14  dB  and  DNRavg  =  36.60  dB,  where  these  averages  are  taken 
across  all  surveillance  and  reference  channels,  respectively.  Therefore,  to  achieve 
(SNRavg,  DNRavg)  =  (—15,15)  dB,  every  target-path  signal  is  scaled  by  +28.14  dB 
(—43.14  +  28.14  =  —15  dB),  and  every  direct-path  signal  is  scaled  by  —21.60  dB 
(36.60  —  21.60  =  15  dB).  The  other  two  SNR-DNR  scenarios  are  achieved  similarly. 

4. 5. 2.1  Low-p  Region. 

Let  (SNRavg,  DNRavg)  =  (—5,  —35)  dB,  which  is  within  the  low-p  region.  Consider, 
first,  a  stationary  target  scenario,  i.e.,  t  =  [ix,iy]  =  02.  Fig.  31  depicts  the  RS- 
GLRT  and  SS-GLRT  statistics,  £,rs  and  £ss,  respectively,  as  a  function  of  hypothesized 
position  p  when  the  hypothesized  velocity  is  matched  to  the  true  target  velocity, 
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Figure  31.  Low-p  Ambiguity  of  the  (a)  RS-GLRT  and  (b)  SS-GLRT  for  a  Stationary 
Target 

p  =  t  =  o2.  As  shown,  both  statistics  exhibit  a  similar  response.  In  particular,  both 
test  statistics  peak  at  the  true  target  position,  and  the  effect  of  beamforming  during 
surveillance  channel  formation  is  evident  in  the  angular  masking  of  their  responses 
with  respect  to  the  receivers.  In  addition,  the  orientation  of  their  main  ambiguity 
responses  in  the  target  vicinity  are  aligned  with  the  target  iso-TDOA  hyperbolas, 
which  are  depicted  as  red  dashed-dotted  lines.  Each  iso-TDOA  contour  represents 
the  positions  at  which  the  TDOA  of  the  hypothesized  state  equals  the  TDOA  of  the 
actual  target  state  with  respect  to  a  given  receiver-receiver  pair,  i.e. ,  A£^jk  =  A 
for  the  jkth.  receiver-receiver  pair.  This  is  illustrated  in  Fig.  32a,  which  depicts  the 
TDOA  A Tp’jk  =  A t^k / fs  as  a  function  of  p  for  i  =  1,  j  =  2,  and  k  =  3.  The  target 
iso-TDOA  contour  At/’22  =  At/’23  is  denoted  by  the  red  dash-dotted  hyperbola. 

Next,  consider  a  moving  target  with  velocity  t  =  [tx,iy]  =  [—375,  375]  m/s. 
Fig.  33  depicts  the  RS-GLRT  and  SS-GLRT  statistics,  £rs  and  £,ss,  respectively,  as 
a  function  of  hypothesized  position  p  when  the  hypothesized  velocity  is  matched  to 
the  true  target  velocity,  p  =  t.  As  in  the  stationary  target  scenario,  both  statistics 
exhibit  similar  ambiguities.  However,  in  contrast  to  the  stationary  target  scenario, 


(a)  TDOA  A r^k  =  M^k/fs  (ns) 


(b)  FDOA  (|;)  A^A:(p,p) 


.  (kHz) 

P=t 


Figure  32.  Example  TDOA  and  FDOA  Plots  for  i  =  1,  j  =  2,  and  k  =  3 


the  ambiguity  peaks  are  sharpened  around  the  true  target  location  due  to  FDOA 
mismatch  in  addition  to  TDOA  mismatch.  This  is  seen  by  the  manner  in  which  the 
peak  responses  are  aligned  with  the  target  iso-FDOA  contours,  which  are  depicted 
in  Fig.  33  by  blue  dashed  lines.  Each  iso-FDOA  contour  represents  the  positions  at 
which  the  FDOA  of  the  hypothesized  state  equals  the  FDOA  of  the  actual  target  state 
with  respect  to  a  given  receiver-receiver  pair,  i.e. ,  =  A/yG7  for  the  jkth  receiver- 

receiver  pair.  This  is  illustrated  in  Fig.  32b,  which  depicts  the  FDOA  (A)  Az/*jfe  in 
kHz  as  a  function  of  p  for  i  —  1,  j  =  2,  and  k  —  3,  and  assuming  p  =  t.  The  target 
iso-FDOA  contour  AzA23  =  A v]'2i  is  denoted  in  this  figure  by  the  blue  dashed  line. 

The  similarity  of  the  RS-GLRT  and  SS-GRLT  responses  for  both  stationary  and 
moving  targets  indicates  that  the  ambiguity  characteristics  of  the  RS-GLRT  are  dom¬ 
inated  by  surveillance-surveillance  processing  in  the  low-p  region.  This  is  consistent 
with  the  detection  sensitivity  results  in  Sec.  4.5.1.  It  is  interesting  to  note,  however, 
that  the  RS-GLRT  ambiguity  response  rolls-off  more  rapidly  away  from  the  peak  re¬ 
sponse  than  the  SS-GLRT.  This  is  seen  by  noting  that  14  dB  and  4  dB  of  dynamic 
range  are  shown  for  £rs  and  £,ss,  respectively,  in  Figs.  31  and  33. 
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Figure  33.  Low-p  Ambiguity  of  the  (a)  RS-GLRT  and  (b)  SS-GLRT  for  a  Moving 
Target 


4. 5. 2. 2  High-DNR  Region. 

Next,  let  (SNRavg,  DNRavg)  =  (—15,15)  dB,  which  is  within  the  high-DNR  re¬ 
gion  as  illustrated  in  Fig.  28.  Consider,  first,  a  stationary  target  scenario,  i.e., 
t  =  [tx,iy]  =  02.  Fig.  34  depicts  the  ambiguity  responses  of  the  RS-GLRT  and  MF- 
GLRT  statistics,  and  £m/,  respectively,  as  a  function  of  hypothesized  position  p 
when  the  hypothesized  velocity  is  matched  to  the  true  target  velocity,  p  =  t  =  02.  As 
shown,  the  ambiguity  responses  of  both  statistics  are  nearly  identical.  In  addition, 
the  peak  responses  are  aligned  along  the  target  iso-range  ellipses,  which  are  depicted 
as  green  dash-dotted  lines.  Each  iso-range  contour  represents  the  positions  at  which 
the  bistatic  range  of  the  hypothesized  state  equals  the  bistatic  range  of  the  actual 
target  state  with  respect  to  a  given  transmitter- receiver  pair,  i.e.,  for  the 

ij'tli  transmitter-receiver  pair.  This  is  illustrated  in  Fig.  35a  for  one  transmitter- 
receiver  pair,  which  depicts  the  bistatic  range  ctj / fs  in  km  as  a  function  of  p  for 
i  —  1  and  j  =  2.  The  target  iso-range  contour  ft2  =  i]2  is  denoted  in  this  figure  by 
the  green  dash-dotted  ellipse. 
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Figure  34.  High-DNR  Ambiguity  of  the  (a)  RS-GLRT  and  (b)  MF-GLRT  for  a  Sta¬ 
tionary  Target 
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Figure  35.  Bistatic  Range  and  Bistatic  Doppler  for  i  =  1  and  j  =  2 


Next,  consider  a  moving  target  with  velocity  t  =  [tx,iy]  =  [—375,  375]  m/s. 
Fig.  36  depicts  the  RS-GLRT  and  MF-GLRT  statistics,  £rs  and  t;ss,  respectively,  as 
a  function  of  hypothesized  position  p  when  the  hypothesized  velocity  is  matched  to 
the  true  target  velocity,  p  =  t.  Zoomed  versions  of  both  statistics  in  the  immediate 
target  vicinity  are  also  depicted  to  give  a  close  view  of  the  peak  ambiguity  responses. 


(a)  £rS(p,p)|p=t  (dB) 


(b)  £m/(P’P)lp=t  (dB) 


(c)  Crs(p,p)|p=t  (dB)  -  Target  Vicinity  (d)  £m/(P- P)lp=t  (dB)  ‘  Target  Vicinity 


Figure  36.  High-DNR  Ambiguity  of  the  (a,c)  RS-GLRT  and  (b,d)  MF-GLRT  for  a 
Moving  Target 


As  in  the  stationary  target  scenario,  both  statistics  exhibit  nearly  identical  ambigu¬ 
ities.  In  addition,  the  peak  responses  are  sharpened  around  the  true  target  location 
in  comparison  to  Fig.  34  due  to  the  effect  of  bistatic  Doppler  mismatch  in  addition 
to  bistatic  range  mismatch.  This  is  seen  by  the  manner  in  which  the  peak  responses 
only  exist  in  the  immediate  vicinity  of  both  the  target  iso-range  and  iso-Doppler  con¬ 
tours  (the  iso-Doppler  contours  are  depicted  by  cyan  dashed  lines).  Each  iso-Doppler 
contour  represents  the  positions  at  which  the  bistatic  Doppler  of  the  hypothesized 
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state  equals  the  bistatic  Doppler  of  the  actual  target  state  with  respect  to  a  given 
transmitter-receiver  pair,  i.e.,  =  vlt3  for  the  i/th  transmitter-receiver  pair.  This  is 

illustrated  in  Fig.  35b,  which  depicts  the  bistatic  Doppler  (^)  in  kHz  as  a  function 
of  p  for  i  —  1  and  j  =  2,  assuming  p  =  t.  The  target  iso-Doppler  contour  z//2  =  z/t12 
is  denoted  in  this  figure  by  the  cyan  dashed  line. 

The  near  equality  of  the  RS-GLRT  and  MF-GLRT  responses  for  both  stationary 
and  moving  targets  indicates  that  the  ambiguity  characteristics  of  the  RS-GLRT  are 
dominated  by  reference-surveillance  processing  in  the  high-DNR  region.  Again,  this  is 
consistent  with  the  detection  sensitivity  results  in  Sec.  4.5.1,  in  which  the  sensitivity 
of  the  RS-GLRT  approached  that  of  the  MF-GLRT  in  the  high-DNR  region. 

4. 5. 2. 3  Transition  Region. 

Finally,  let  (SNRavg,  DNRavg)  =  (-10,-10)  dB,  which  is  within  the  transition 
region  as  shown  in  Fig.  28.  Fig.  37  depicts  the  ambiguity  response  of  the  RS-GLRT 
statistic,  £rs,  for  the  stationary  target  scenario,  in  which  p  =  t  =  O2.  As  shown, 
the  ambiguity  response  exhibits  significant  sidelobes  along  both  the  iso-TDOA  and 
the  iso-bistatic  range  contours.  Fig.  38  similarly  depicts  the  RS-GLRT  ambiguity 
response  for  the  moving  target  scenario,  in  which  p  =  t  =  [—375,  375]  m/s.  In  com¬ 
parison  to  the  stationary  target  scenario  in  Fig.  37,  the  moving  target  peak  ambiguity 
response  in  Fig.  38  is  additionally  constrained  to  exist  within  the  immediate  vicinity 
of  the  iso-FDOA  and  iso-bistatic  Doppler  contours.  Thus,  the  RS-GLRT  ambiguity 
response  in  the  transition  region  exhibits  properties  of  the  RS-GLRT  ambiguity  re¬ 
sponses  in  both  the  low-p  and  high-DNR  regions,  in  that  it  appears  to  be  constrained 
in  TDOA,  FDOA,  bistatic  range,  and  bistatic  Doppler4.  This  indicates  that  both 
surveillance-surveillance  and  reference-surveillance  processing  contribute  significantly 

4RS-GLRT  ambiguity  is  also  implicitly  constrained  in  angle  due  to  the  beamforming  operation 
during  surveillance  channel  formation. 
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Figure  37.  Transition  Region  Ambiguity  of  the  RS-GLRT  for  a  Stationary  Target 
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Figure  38.  Transition  Region  Ambiguity  of  the  RS-GLRT  for  a  Moving  Target 


to  the  ambiguity  properties  of  the  RS-GLRT  in  the  transition  region.  Again,  this 
result  is  consistent  with  the  detection  sensitivity  results  in  Sec.  4.5.1,  where  the  tran¬ 
sition  region  RS-GLRT  sensitivity  was  shown  to  be  intermediate  to  its  sensitivity  in 
the  low-p  and  high-DNR  regions. 
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4.6  Discussion 


The  results  presented  in  Sec.  4.5  demonstrate  that  the  RS-GLRT  bridges  the  gap 
between  the  detection  and  ambiguity  performance  of  the  SS-GLRT  and  MF-GLRT 
detectors  in  a  manner  that  depends  on  the  quality  of  the  direct-path  reference  signals. 
This  result  is  expected  from  the  analysis  presented  in  Sec.  4.4,  where  it  was  shown 
that  the  RS-GLRT  statistic  encompasses  the  signal  processing  operations  that  are 
characteristic  of  both  PSL  and  AMR,  termed  reference-surveillance  and  surveillance- 
surveillance  processing,  respectively.  Both  are  discussed  in  turn. 

First,  reference-surveillance  processing  entails  calculation  of  the  pairwise  reference- 
surveillance  CAFs.  This  is  the  characteristic  operation  of  matched  filtering  in  AMR, 
and  in  the  conventional  approach  to  PMR  detection  using  reference  signals.  Geomet¬ 
rically,  each  Tx-Rx  pair  constrains  the  target  location  in  bistatic  range  and  bistatic 
Doppler,  as  discussed  in  Sec.  4.4.2  and  illustrated  in  Sec.  4.5.2.  The  detection  results 
in  Sec.  4.5.1  show  that  this  type  of  processing  is  dominant  for  RS-GLRT  detection 
when  the  DNR  is  positive,  and  that  coherent  integration  gain  is  realized  with  increas¬ 
ing  signal  length  in  this  region.  Furthermore,  RS-GLRT  performance  asymptotically 
approaches  MF-GLRT  performance  with  increasing  DNRavg.  This  is  expected  from 
the  analysis  in  Sec.  4.3.2,  where  it  was  shown  that  C,rs  approximates  E,mf  within  a  scalar 
constant  under  the  high-DNR  and  high-p  conditions.  This  justifies  the  conventional 
approach  to  PMR  detection,  in  which  the  direct-path  signal  in  the  reference  chan¬ 
nel  is  used  as  the  reference  in  an  approximate  matched  filtering  operation,  provided 
DNRavg  is  positive. 

Second,  surveillance-surveillance  processing  entails  calculation  of  the  pairwise 
surveillance-surveillance  CAFs.  This  is  the  characteristic  operation  of  PSL  and  PMR 
without  a  reference  signal.  Geometrically,  each  Rx-Rx  pair  constrains  the  target  lo¬ 
cation  in  TDOA  and  FDOA,  as  discussed  in  Sec.  4.4.1  and  illustrated  in  Sec.  4.5.2. 
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The  detection  results  in  Sec.  4.5.1  show  that  this  type  of  processing  is  dominant  for 
RS-GLRT  detection  when  SNRavg  exceeds  DNRavg  by  approximately  5  dB,  i.e.,  when 
p  <  5  dB.  In  this  region,  the  RS-GLRT  realizes  a  non-coherent  integration  gain  with 
increasing  signal  length,  similar  to  the  SS-GLRT.  Note  that  low-p  scenarios  are  ex¬ 
pected  to  be  rare  in  practice,  as  both  theoretical  analysis  and  experimental  results 
indicate  p  is  often  exceedingly  (and  problematically)  large  [27,  33]. 

It  is  interesting  to  note  that  the  asymptotic  detection  performance  of  the  RS- 
GLRT  in  the  low-p  region  is  slightly  inferior  to  that  of  the  SS-GLRT.  This  might 
be  explained  by  the  fact  that  the  RS-GLRT  is  adversely  affected  at  low  DNR  by 
extremely  noisy  reference  signals  that  are  completely  ignored  by  the  SS-GLRT.  In¬ 
corporating  these  additional  reference  signals  degrades  the  implicit  estimation  of  the 
unknown  signals  that  is  performed  in  calculation  of  the  RS-GLRT  statistic.  It  also 
necessitates  the  implicit  estimation  of  the  unknown  reference  channel  coefficients, 
which  are  additional  nuisance  parameters  that  do  appear  in  the  SS-GLRT  formu¬ 
lation.  Therefore,  this  result  suggests  that  there  are  cases  in  which  it  is  better  to 
ignore  additional  information  that  pertains  to  nuisance  parameters  within  the  prob¬ 
lem  formulation.  This  is  consistent  with  a  result  presented  by  Ramirez  et  al.  in  [70] , 
where  it  was  shown  that  ignoring  a  priori  knowledge  about  the  covariance  structure 
of  random  observations  in  a  particular  detection  problem  is  optimal  under  low-SNR 
conditions.  In  their  application,  incorporating  this  a  priori  knowledge  via  the  deriva¬ 
tion  of  a  GLRT  actually  degrades  performance  because  estimation  of  the  resulting 
nuisance  parameters  cannot  be  performed  accurately  under  low-SNR  conditions. 

These  results  clarify  the  fundamental  role  of  direct-path  signals  in  PMR  detec¬ 
tion.  On  one  hand,  these  signals  might  not  be  expected  to  aid  detection  perfor¬ 
mance  because  they  do  not  differ  between  the  hypotheses  in  the  detection  problem 
of  Eqs.  (124)-(125).  On  the  other  hand,  they  provide  valuable  information  about 
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the  unknown  transmitted  signals,  which  are  considered  nuisance  parameters  in  our 
formulation.  In  other  words,  they  shed  light  on  what  the  target  return  should  look 
like  by  providing  (potentially)  high-quality  estimates  of  the  unknown  transmitted 
signals  to  the  PMR  system.  Fig.  29  demonstrates  that  the  benefit  provided  by  the 
direct-path  signals  in  PMR  detection  depends  on  the  length  of  the  received  signals. 
This  benefit  may  be  measured  in  terms  of  a  direct-path  gain,  Gdp,  defined  as  the 
difference  between  the  MF-GLRT  and  SS-GLRT  integration  gain  (Gint)  curves  in 
Fig.  30.  The  SS-GLRT  is  chosen  rather  than  the  RS-GLRT  because  the  SS-GLRT 
represents  the  best  possible  PMR  detection  performance  that  can  be  achieved  when 
direct-path  signals  are  ignored  completely.  As  shown  in  Fig.  30,  Gdp  ranges  from  0 
dB  for  L  —  1  to  16.47  dB  for  L  =  10000.  Thus,  direct-path  signals  may  significantly 
improve  detection  performance,  particularly  for  long  signal  lengths. 

4.7  Conclusion 

This  chapter  has  presented  a  GLRT  for  centralized  passive  MIMO  radar  detection. 
The  performance  of  this  detector  has  been  compared  to  similar  GLRTs  for  centralized 
detection  in  active  MIMO  radar  and  passive  source  localization  sensor  networks.  It 
has  been  shown  that  the  detection  and  ambiguity  performance  of  the  passive  MIMO 
radar  GLRT  varies  between  that  of  GLRTs  for  active  MIMO  radar  and  passive  source 
localization  in  a  manner  that  depends  on  the  average  direct-path  SNR,  which  is  a 
measure  of  the  degree  of  knowledge  about  the  (a  priori  unknown)  transmitted  signals. 
With  high  direct-path  SNR,  passive  MIMO  radar  detection  approaches  that  of  active 
MIMO  radar.  With  low  direct-path  SNR,  passive  MIMO  radar  detection  is  similar  to 
that  of  passive  source  localization.  In  this  way,  passive  MIMO  radar  generalizes  both 
active  MIMO  radar  and  passive  source  localization,  unifying  them  within  a  common 
theoretical  framework.  Such  a  framework  is  discussed  further  in  Chapter  V. 
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V.  Unified  Detection  Framework 


This  chapter  introduces  a  unified  theoretical  framework  for  detection  in  active  and 
passive  distributed  RF  sensor  networks.  This  framework,  summarized  in  Fig.  39,  en¬ 
compasses  active  and  passive  sensing  networks  including  AMR,  PMR,  and  PSL.  The 
main  features  of  this  framework  have  already  been  introduced  through  the  derivation, 
interpretation,  and  numerical  analysis  of  the  SS-GLRT  and  RS-GLRT  detectors  in 
Chapters  Ill  and  IV,  respectively.  Drawing  from  these  results,  this  chapter  shows 
that  PMR  is  the  key  to  linking  AMR  and  PSL  sensor  networks,  which  have  tradition¬ 
ally  been  regarded  as  distinct.  Thus,  this  framework  provides  fundamental  insight 
into  the  natures  of  active  and  passive  distributed  RF  sensing.  Sec.  5.1  introduces 
the  unified  framework  by  showing  how  AMR,  PMR,  and  PSL  are  related  by  simple 
transformations  of  their  respective  signal  environments.  Based  on  these  relationships, 
signal  models  are  defined  for  AMR,  PMR  and  PSL  in  Sec.  5.2.  The  detectors  that 
result  from  these  signal  models  are  discussed  in  Sec.  5.3,  as  well  as  the  underlying 
signal  processing  operations  performed  by  each.  Conclusions  follow  in  Sec.  5.4. 

5.1  Transformations 

The  signal  environments  for  AMR,  PMR,  and  PSL  are  closely  related.  This  is 
shown  in  Fig.  40,  which  depicts  the  signal  environments  for  AMR,  PMR,  PMR  with¬ 
out  direct-path  references,  and  PSL  sensor  networks,  as  well  as  the  transformations 
that  relate  them.  These  transformations  are  discussed  in  turn.  First,  as  discussed  in 
Chapter  I,  AMR  and  PMR  are  distinguished  from  each  other  depending  on  whether 
the  transmitters  are  cooperative  (AMR)  or  non-cooperative  (PMR).  This  is  equiva¬ 
lent  to  whether  the  transmit  signals  are  a  priori  known  (AMR)  or  a  priori  unknown 
(PMR).  Although  AMR  may  have  co- located  transmitters  and  receivers,  while  in 
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Figure  39.  A  Unified  Theoretical  Framework  for  Detection  in  Active  and  Passive  Distributed  RF  Sensor  Networks 


PMR 


PMR  without  References 


Figure  40.  Transformations  between  AMR,  PMR,  and  PSL  Signal  Environments 


PMR  the  transmitters  and  receivers  are  necessarily  separate,  this  distinction  is  in¬ 
significant  from  a  detection  processing  perspective1. 

Although  the  transmit  signals  are  a  priori  unknown  in  PMR,  imperfect  estimates 
of  these  transmit  signals  are  provided  to  the  PMR  system  by  direct-path  signals. 
The  quality  of  these  estimates  is  quantified  by  the  average  DNR;  high-quality  and 
low-quality  estimates  are  provided  by  high-DNR  and  low-DNR  direct-path  signals, 
respectively.  This  defines  a  continuum  quantified  by  the  average  DNR.  At  the  high- 
DNR  extreme  the  transmit  signals  become  known  perfectly;  at  the  low-DNR  extreme 
nothing  is  known  about  the  transmit  signals.  This  low-DNR  extreme  exists  in  PMR 
if  direct-path  signals  are  unavailable  to  the  PMR  receivers,  which  might  occur  if 
there  is  there  is  physical  blockage  of  the  line-of-sight  paths  between  transmitters 

1  Although  this  distinction  will  affect  the  implementation  details  of  AMR  detection  processing, 
e.g.,  the  delay-Doppler  compensation,  it  does  not  fundamentally  alter  the  type  of  processing  (i.e., 
matched  filtering)  supported  by  the  AMR  signal  environment. 
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and  receivers,  or  if  illuminators  were  to  utilize  highly  directive  transmission  and/or 
null  the  receiver  directions.  The  resulting  sensor  network  is  termed  “PMR  without 
direct-path  references,”  which  is  abbreviated  “PMR  without  references”  in  Fig.  40 
and  throughout  the  rest  of  this  chapter. 

Finally,  PMR  without  references  and  PSL  are  closely  related  because  in  both 
only  target-path  signals  are  received  by  the  sensor  network.  They  are  distinguished 
according  to  whether  the  target-path  signals  are  scattered  (PMR  without  references) 
or  emitted  (PSL)  by  the  target.  Accordingly,  these  sensor  networks  are  distinguished 
by  the  type  of  propagation  channel  experienced  by  the  a  priori  unknown  target-path 
signals.  This  distinction  is  insignificant  from  a  detection  perspective.  Specifically, 
it  is  shown  in  Sec.  3.4  that  PMR  without  references  is  mathematically  equivalent  to 
PSL  when  the  effects  of  the  first  leg  of  the  bistatic  target-path  channel  are  absorbed 
into  the  unknown  transmitted  signal. 

5.2  Signal  Models 

In  this  section,  signal  models  are  given  for  AMR,  PMR,  and  PSL.  Due  to  the 
relationship  between  AMR,  PMR,  and  PSL  identified  in  Sec.  5.1,  these  signal  models 
are  derived  directly  from  the  PMR  signal  model  presented  in  Chapter  II. 

5.2.1  Active  MIMO  Radar. 

As  discussed  in  Sec.  5.1,  AMR  is  distinguished  from  PMR  by  the  use  of  cooper¬ 
ative  transmitters  that  transmit  known  signals.  In  addition,  the  transmitters  within 
AMR  are  typically  co-located  with  the  receivers,  although  this  is  not  necessarily  true. 
Consequently,  the  signal  model  is  identical  to  that  of  PMR  except  that  the  transmit 
signals  {u1}  are  assumed  known,  and  the  ith  baseband  signal  at  the  nth  element  of 
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the  jth  AMR  receiver  array  is  given  by  Eq.  (50)  in  Chapter  II  as 

=  ~'j  V'j  u 1  +  <■">"  ^  V\j  u'  +  n>3  (202) 

where  u*  is  the  known  transmit  signal  from  the  ith  transmitter,  T>lJ  =  T)L,(£lfj ,  i/3) 
and  T>ltJ  =  T>Li(£l3,  vlt3)  for  brevity,  and  all  the  remaining  variables  are  defined  as  in 
(50).  Note  that  (202)  is  a  CW  formulation  for  AMR,  whereas  AMR  is  often  studied 
in  the  context  of  pulsed  operation.  However,  (202)  also  applies  to  pulsed  operation 
if  u*  is  interpreted  as  a  pulse  train  that  represents  an  entire  CPI.  Note  also  that  the 
direct-path  signal  is  often  excluded  in  AMR  signal  model  formulations  because  it  is 
implicitly  assumed  that  the  receivers  are  blanked  while  any  pulsed  direct-path  signal 
is  incident. 


5.2.2  Passive  MIMO  Radar. 

The  PMR  signal  model  is  derived  in  detail  in  Chapter  II,  where  the  ith  baseband 
signal  at  the  nth  element  of  the  jth  PMR  receiver  array  is  given  by  Eq.  (50)  as 


4  =  4/  <3"‘"',(i'  Vj  u'  +  e**" (t)  V\j 


u'  +  n 


(203) 


where  u*  is  the  unknown  transmit  signal  from  the  ith  transmitter,  ,  z/)J) 

and  T>\3  =  T)Li[t3 ,v\3)  for  brevity,  and  all  the  remaining  variables  are  defined  as  in 
(50).  The  signal  model  for  PMR  without  references  is  similarly  given  by 

4=~l?^Viiu‘  +  n «  (204) 


which  is  equivalent  to  (203)  except  the  direct-path  signal  is  absent. 
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5.2.3  Passive  Source  Localization. 


As  discussed  in  Sec.  5.1,  PMR  without  references  and  PSL  are  distinguished  by 
whether  the  target-path  signals  are  scattered  (PMR  without  references)  or  emitted 
(PSL)  by  the  target.  This  distinction  may  be  defined  mathematically  by  dividing 
the  bistatic  target-path  channel  coefficient  and  delay-Doppler  operator  T>lt]  into 
components  that  represent  each  leg  of  the  bistatic  propagation  channel,  i.e.,  = 

and  'E’t’  —  i  as  discussed  in  Secs.  3.1  and  3.4  of  Chapter  III.  Then,  the 

ith  baseband  signal  at  the  nth  element  of  the  jth  PSL  receiver  array  is  given  by 

43'  =  7,‘i^"(t)A2U,  +  n«  (205) 


where  u*  is  the  ith  unknown  signal  emitted  by  the  target,  and  7 is  the  channel 
coefficient  between  the  target  and  the  jth  receiver,  dehned  by  (53)  in  Sec.  3.1  as 


h(r-^nJ2(  t)/c) 


\ 


tu*)  yigw 

(4vrRJ2(t))2 


(206) 


where  the  terms  in  (206)  are  dehned  in  Sec.  3.1.  Note  that  in  (53)  there  is  an  a y' 
that  represents  the  bistatic  rehectivity  with  respect  to  the  ij th  bistatic  target-path 
channel.  In  (206),  this  term  has  been  replaced  by  yj the  effective  radiated 
power  in  the  direction  of  the  jth  receiver,  to  rehect  the  difference  between  PSL  and 
PMR  without  references. 


5.3  Detectors 

When  the  centralized  target  detection  problem  is  formulated  for  AMR,  PMR, 
PMR  without  references,  and  PSL  sensor  networks  using  the  signal  models  presented 
in  Sec.  5.2,  the  resulting  detectors  exhibit  similarities  that  rehect  the  similarities  in 
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their  respective  signal  models.  These  detectors,  depicted  in  Fig.  39,  are  the  MF- 
GLRT,  RS-GLRT,  and  SS-GLRT  detectors  presented  in  Chapters  III  and  IV.  These 
are  defined  and  discussed  in  the  following  sections. 

5.3.1  Matched  Filter  GLRT. 

First,  the  matched  filter  GLRT  (MF-GLRT)  is  the  GLRT  for  centralized  detection 
in  AMR  sensor  networks  that  results  from  the  AMR  signal  model  in  (202).  This  GLRT 
is  given  by  Eq.  (150)  in  Sec.  4.2.3  as 

,  N k  Nr  H 

=  i «-/  <207) 

i.i  j.i 

where  s*-7  =  is  the  surveillance  signal  slJ  after  removal  of  the  delay  and 

Doppler  that  results  from  propagation  along  the  ij th  hypothesized  target-path  chan¬ 

nel  with  respect  to  the  cell  under  test.  As  shown  in  (207),  the  test  statistic  is 
formed  by  delay-Doppler  compensating  each  surveillance  signal,  matched  filtering, 
and  then  non-coherently  integrating  the  matched  filter  outputs  across  all  bistatic 
transmitter-receiver  pairs.  This  signal  processing  structure  is  common  in  AMR  de¬ 
tection  [24,  43].  The  test  statistic  £m/  follows  the  non-central  and  central  chi-squared 
distributions  under  "Hi  and  "Ho,  respectively  (Sec.  4.3.1).  The  integration  gain  real¬ 
ized  by  the  MF-GLRT  grows  coherently  with  increasing  signal  length,  as  expected  for 
matched  filtering  (Sec.  4.5.1).  In  addition,  the  ambiguity  properties  of  the  MF-GLRT 
can  be  explained  in  terms  of  bistatic  range,  bistatic  Doppler,  and  (for  array  receivers) 
AOA  (Sec.  4.5.2).  Note  that  matched  filtering  is  only  possible  because  the  transmit 
signals  (u'}  are  known  in  the  AMR  detection  problem  formulation. 
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5.3.2  Surveillance-Surveillance  GLRT. 

Next,  consider  the  surveillance-surveillance  GLRT  (SS-GLRT).  As  discussed  in 
Sec.  3.2,  this  detector  is  the  GLRT  for  centralized  detection  in  PMR  networks  without 
direct-path  references.  Due  to  the  mathematical  equivalence  between  PMR  networks 
without  references  and  PSL  networks,  it  is  also  the  GLRT  for  centralized  detection 
in  PSL  networks.  The  SS-GRLT  is  given  by  Eq.  (79)  in  Sec.  3.2  as 

1  ^  rti 

^  =  -EAl(G-)  ^  ««  (208) 

tr  *° 

where  Ai(-)  denotes  the  largest  eigenvalue  of  its  matrix  argument,  and  G*s  is  a  Gram 
matrix  with  jkt\i  element  defined  by 

[GLR  =  (209) 

It  is  shown  in  Sec.  4.4.1  that  the  entries  of  G*s  may  be  found  by  sampling  the  CAFs  be¬ 
tween  each  pair  of  surveillance  signals.  This  is  referred  to  as  surveillance-surveillance 
processing.  The  Gram  matrix  G*s  is  a  complex  Wishart  matrix  that  is  non-central 
under  'H\  and  central  under  770,  and  the  exact  distribution  of  the  test  statistic  £ ss 
is  found  in  the  recent  random  matrix  theory  literature  (Sec.  3.3).  In  comparison  to 
the  MF-GLRT,  the  SS-GLRT  achieves  non-coherent  rather  than  coherent  integration 
gain  with  increasing  signal  length  (Sec.  3.5. 1.2  and  4.5.1).  In  addition,  the  ambiguity 

properties  of  the  SS-GLRT  can  be  explained  in  terms  of  TDOA,  FDOA,  and  (for  array 

receivers)  AOA  (Secs.  3.5.2  and  4.5.2).  Finally,  transmitters  and  receivers  contribute 
asymmetrically  to  detection  sensitivity  in  that  adding  an  additional  receiver  increases 
detection  sensitivity  more  than  adding  an  additional  transmitter  (Sec.  3.5. 1.1).  This 
contrasts  with  the  MF-GLRT,  which  benefits  from  transmitters  and  receivers  equally. 
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5.3.3  Reference-Surveillance  GLRT. 

Finally,  consider  the  RS-GLRT.  As  shown  in  Sec.  4.2.1,  this  detector  is  the  GLRT 
for  centralized  detection  in  PMR  networks  with  direct-path  references.  It  is  given  by 
Eq.  (148)  in  Sec.  4.2.1  as 

1  Nt  y 

«™=-E(MG‘l)-MGt))  S'  «r.  (2W) 

a  J  n0 

where  G)  is  a  Gram  matrix  with  block  structure 

Gi  r^i 
**  vj 

G\  =  (211) 

Qi  Qi 

G*s  is  defined  by  (209),  G*s  =  (G®r)H,  G*r  and  G*r  are  dehned 

[G;r],t  =  (s«)"s“  (212) 

[GtR  =  (§?)"§“  (213) 

and  slJ  =  (Vf)Hslf  is  the  reference  signal  slf  after  removal  of  the  delay  and  Doppler 
shift  due  to  direct-path  propagation  from  the  ith  transmitter  to  the  jth  receiver.  It 
is  shown  in  Sec.  4.4.2  that  the  entries  of  G*r  may  be  found  by  sampling  CAFs  be¬ 
tween  reference  and  surveillance  signals.  The  calculation  of  such  reference-surveillance 
CAFs,  referred  to  as  reference-surveillance  processing,  can  be  interpreted  as  matched 
filtering  with  a  noisy  reference.  This  operation  is  the  basis  for  the  conventional  ap¬ 
proach  to  PMR  detection,  described  in  Sec.  1.4.1. 

The  Gram  matrix  G*  is  a  non-central  complex  Wishart  matrix  under  both  "Hi  and 
"Ho,  however  the  random  variables  Ai(GG  and  Ai(G*r)  in  (210)  are  not  independent 
due  to  the  dependence  of  G^  on  G*r.  This  case  does  not  appear  to  have  been  addressed 
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in  the  random  matrix  literature;  consequently,  the  probability  distribution  of  the 
statistic  £rs  is  not  known.  However,  the  distribution  of  £ rs  can  be  approximated 
under  certain  conditions,  described  shortly. 

The  appearance  of  both  Glss  and  G*r  in  (211)  suggests  that  the  RS-GLRT  exploits 
both  reference-surveillance  and  surveillance-surveillance  processing.  Comparison  of 
the  RS-GRLT,  MF-GRLT,  and  SS-GLRT  detectors  in  Sec.  4.5  shows  that  this  is 
true.  Specifically,  RS-GLRT  detection  performance,  in  terms  of  both  sensitivity  and 
ambiguity,  varies  between  that  of  the  MF-GLRT  and  SS-GRLT  as  a  function  of 
the  average  target-path-to-noise  ratio  (SNR)  and  the  average  direct-path-to-noise 
ratio  (DNR).  It  is  shown  in  Sec.  4.5  that  the  performance  of  the  RS-GLRT  can  be 
divided  into  three  distinct  SNR-DNR  regions:  a  high-DNR  region,  in  which  reference- 
surveillance  processing  is  dominant;  a  low- jo  region  (where  jo  denotes  the  average 
direct-path  to  target-path  power  ratio),  in  which  surveillance-surveillance  processing 
is  dominant;  and  a  transition  region  between  the  high-DNR  and  low-p  regions  where 
both  types  of  processing  are  significant.  These  regions  are  now  discussed  in  turn. 

In  the  high-DNR  region,  reference-surveillance  processing  (i.e.,  noisy  matched  fil¬ 
tering)  is  dominant,  and  RS-GLRT  performance  approaches  MF-GLRT  performance. 
Consequently,  the  RS-GLRT  realizes  coherent  integration  gain  with  signal  length 
(Sec.  4.5.1),  and  its  ambiguity  can  be  explained  in  terms  of  bistatic  range,  bistatic 
Doppler,  and  (for  array  receivers)  AOA  (Sec.  4.5.2).  Furthermore,  assuming  high- 
DNR  and  high- jo,  the  RS-GLRT  statistic,  £rs,  is  approximately  equal  to  a  scaled  ver¬ 
sion  of  the  MF-GLRT  statistic,  £m/  (Sec.  4.3.2).  Therefore,  £rs  approximately  follows 
the  non-central  and  central  chi-squared  distributions  under  "Hi  and  "Ho,  respectively, 
under  high-DNR  and  high-jo  conditions. 

In  the  low-p  region,  surveillance-surveillance  processing  is  dominant,  and  RS- 
GLRT  performance  approaches  SS-GLRT  performance.  Thus,  the  RS-GLRT  realizes 
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non-coherent  integration  gain  with  signal  length  (Sec.  4.5.1),  and  its  ambiguity  can 
be  explained  in  terms  of  TDOA,  FDOA,  and  (for  array  receivers)  AOA  (Sec.  4.5.2). 

Finally,  the  RS-GLRT  detector  exhibits  properties  of  both  surveillance-surveillance 
and  reference-surveillance  processing  in  the  transition  region.  In  particular,  its  detec¬ 
tion  sensitivity  is  intermediate  to  that  of  the  MF-GLRT  and  SS-GLRT  (Sec.  4.5.1), 
and  its  ambiguity  properties  are  influenced  by  bistatic  range,  bistatic  Doppler,  TDOA, 
FDOA,  and  (for  array  receivers)  AOA  (Sec.  4.5.2).  Note  that  in  the  transition  region 
the  sensitivity  and  ambiguity  performance  of  the  RS-GLRT  vary  smoothly  between 
the  sensitivity  and  ambiguity  of  the  SS-GLRT  and  MF-GLRT  with  increasing  DNR. 

5.4  Conclusion 

This  chapter  demonstrates  that  PMR  links  PSL  and  AMR  sensor  networks  within 
a  common  theoretical  framework.  This  framework  is  summarized  in  Fig.  39.  As 
discussed  in  Sec.  5.1,  imperfect  estimates  of  the  unknown  transmit  signals  are  provided 
to  the  PMR  system  by  direct-path  signals.  The  quality  of  this  estimate  is  quantified  by 
the  DNR;  high-quality  and  low-quality  estimates  are  provided  by  high-DNR  and  low- 
DNR  direct-path  signals.  This  defines  a  continuum  that  conceptually  spans  the  gap 
between  AMR  and  PSL  sensor  networks.  At  the  high-DNR  extreme,  PMR  approaches 
AMR  in  that  the  a  priori  unknown  signals  may  be  estimated  accurately  via  reception 
of  high-DNR  direct-path  signals.  At  the  low-DNR  extreme,  PMR  approaches  PSL  in 
that  there  is  no  auxiliary  source  of  information  about  the  a  priori  unknown  transmit 
signals  available  to  the  PMR  system  apart  from  the  target-path  signal  itself.  In  this 
way,  PMR  links  PSL  and  AMR  sensor  networks,  thereby  unifying  them  within  a 
common  theoretical  framework.  This  result  provides  insight  into  the  fundamental 
natures  of  active  and  passive  distributed  RF  sensing. 
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VI.  Conclusion 


This  dissertation  advances  the  theory  of  passive  MIMO  radar  detection.  The  main 
contributions  of  this  research  are  summarized  in  Sec.  6.1.  Promising  theoretical  and 
experimental  extensions  to  this  research  are  discussed  in  Sec.  6.2. 

6.1  Summary  of  Contributions 

As  stated  in  Sec.  1.6,  the  main  contributions  of  this  research  are: 

1.  The  derivation  and  statistical  characterization  of  a  novel  detector  for  the  cen¬ 
tralized  PMR  detection  without  direct-path  references  problem  (Chapter  111) 

2.  The  derivation  and  statistical  characterization  of  a  novel  detector  for  the  cen¬ 
tralized  PMR  detection  with  direct-path  references  problem  (Chapter  IV) 

3.  The  introduction  of  a  unified  framework  for  detection  in  active  and  passive 
distributed  RF  sensor  networks  that  encompasses  active  MIMO  radar,  passive 
MIMO  radar,  and  passive  source  localization  sensor  networks  (Chapter  V) 

6.1.1  PMR  Detection  without  References. 

Chapter  III  introduced  the  SS-GLRT,  a  novel  GLRT  detector  for  centralized  PMR 
detection  without  direct-path  references.  The  formulation  of  this  detector  extends 
the  formulations  of  previous  PMR  detectors  to  account  for  multiple  non-cooperative 
transmitters,  multichannel  (array)  receivers,  non-isotropic  target  scattering,  and  non¬ 
coherent  receivers.  Previous  PMR  detector  formulations  have  considered,  at  most, 
two  out  of  four  of  these  conditions.  These  extensions  are  significant  because  they 
enhance  detection  sensitivity  and  reduce  detection  ambiguity.  The  exact  distribu¬ 
tions  of  the  SS-GLRT  test  statistic  are  identified  under  both  hypotheses  using  recent 
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results  from  RMT,  which  represents  the  first  appearance  of  these  distributions  in  the 
context  of  radar  detection.  These  distributions  show  that  detection  sensitivity  is  only 
a  function  of  the  number  of  received  signal  samples,  the  number  of  transmitters  and 
receivers  in  the  PMR  network,  and  the  average  target-path  SNR.  Furthermore,  nu¬ 
merical  examples  demonstrate  important  properties  of  the  detector,  namely,  that  (a) 
receivers  and  transmitters  contribute  asymmetrically  to  detection  sensitivity  (which 
is  uncommon  in  radar  detection),  and  (b)  integration  gain  grows  non-coherently  with 
increasing  signal  length.  The  SS-GLRT  has  been  submitted  for  publication  in  [38]. 

6.1.2  PMR  Detection  with  References. 

Chapter  IV  introduced  the  RS-GLRT,  a  novel  GLRT  detector  for  centralized  PMR 
detection  with  direct-path  references.  This  detector  extends  the  formulation  of  the 
SS-GLRT  to  include  direct-path  signals  with  arbitrary  DNRs.  This  extension  is 
significant  because  it  quantifies,  for  the  first  time,  the  dependence  of  PMR  detection 
performance  on  the  quality  of  the  direct-path  reference  signals.  Analysis  of  the  RS- 
GLRT  test  statistic  distributions  shows  that  detection  sensitivity  depends  only  on 
the  number  of  received  signal  samples,  the  number  of  transmitters  and  receivers  in 
the  PMR  network,  the  average  target-path  input  SNR,  and  the  average  input  DNR. 
Numerical  simulations  reveal  that  PMR  detection  sensitivity  and  ambiguity  may  be 
divided  into  three  SNR/DNR  regimes:  first,  a  high-DNR  region,  in  which  reference- 
surveillance  processing  is  dominant;  second,  a  low-p  region,  in  which  surveillance- 
surveillance  processing  is  dominant;  and  third,  a  transition  region  between  the  high- 
DNR  and  low-p  regions,  in  which  both  types  of  processing  contribute  significantly  to 
sensitivity  and  ambiguity.  These  regimes  characterize,  for  the  first  time,  how  PMR 
sensitivity  and  ambiguity  depends  on  both  SNR  and  DNR.  The  RS-GLRT  has  been 
submitted  for  publication  in  [37]. 
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6.1.3  Unified  Detection  Framework. 


Finally,  a  novel  theoretical  framework  for  detection  in  active  and  passive  dis¬ 
tributed  RF  sensor  networks  is  proposed.  This  framework,  summarized  in  Fig.  39  of 
Chapter  V,  encompasses  AMR,  PMR,  and  PSL  sensor  networks.  It  identifies  how 
AMR,  PMR,  and  PSL  are  related  by  simple  transformations  of  their  respective  signal 
environments.  These  transformations  reveal  that  PMR  is  the  key  to  linking  AMR  and 
PSL,  which  have  traditionally  been  regarded  as  distinct.  These  relationships  also  sug¬ 
gest  that  PMR  detection  might  exhibit  properties  of  detection  in  both  AMR  and  PSL 
sensor  networks  under  the  appropriate  DNR  conditions.  Analysis  of  the  RS-GLRT 
and  SS-GLRT  detectors  shows  that  this  is  true,  i.e. ,  PMR  detection  performance, 
in  terms  of  both  sensitivity  and  ambiguity,  varies  between  that  of  AMR  and  PSL 
as  a  function  of  the  average  DNR.  In  high-DNR  scenarios,  RS-GLRT  performance 
approaches  AMR  performance.  In  such  scenarios,  reference-surveillance  processing 
(i.e.,  noisy  matched  filtering)  is  dominant.  In  low-p  scenarios  (where  p  is  the  average 
direct-path  to  target-path  signal  power  ratio),  RS-GLRT  performance  approaches  SS- 
GLRT  performance,  which  is  the  GLRT  detector  for  both  PMR  without  references 
and  PSL  sensor  networks.  In  such  scenarios,  surveillance-surveillance  processing  is 
dominant.  Finally,  at  intermediate  DNRs,  RS-GLRT  performance  exhibits  properties 
of  both  AMR  and  PSL.  These  results  are  explained  by  interpreting  PSL  and  AMR  as 
two  extremes  in  terms  of  knowledge  about  the  transmit  signals,  i.e.,  the  signals  are 
entirely  unknown  in  PSL  and  entirely  known  in  AMR.  These  extremes  are  spanned  by 
PMR,  for  which  direct-path  signals  provide  a  varying  degree  of  knowledge  about  the 
unknown  transmit  signals  that  is  quantified  by  the  DNR.  In  this  way,  PMR  may  be 
interpreted  as  generalizing  PSL  and  AMR.  This  result  provides  fundamental  insight 
into  the  natures  of  active  and  passive  distributed  RF  sensing. 
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6.2  Future  Research 


This  work  may  be  extended  in  several  ways.  First,  the  RS-GLRT  detector  could 
be  extended  to  explicitly  account  for  clutter-path  interference.  In  this  work,  the  RS- 
GLRT  formulation  considers  only  receiver  noise  and  direct-path  signals.  Direct-path 
interference  is  mitigated  via  deterministic  nulling  of  the  transmitter  direction  dur¬ 
ing  surveillance  channel  formation.  Consequently,  the  RS-GLRT  results  presented 
in  Chapter  IV  represent  noise- limited  performance,  which  may  be  regarded  as  an 
upper  bound  on  performance.  In  practice,  clutter-path  interference  often  limits  per¬ 
formance  by  masking  target  returns  [19].  This  problem  is  exacerbated  for  GW  signals 
in  comparison  to  pulsed  signals  because  a  GW  signal’s  ambiguity  function  fills  all  of 
range-Dopplcr  space.  Consequently,  the  ambiguity  function  floor  of  strong  near-range 
clutter  returns  may  mask  far-range  targets  [53,  51].  This  limits  the  applicability  of 
pulsed-radar  clutter  mitigation  techniques  such  as  bistatic  space-time  adaptive  pro¬ 
cessing  (STAP),  which  are  formulated  assuming  only  intra-range  cell  clutter.  It  is 
expected  that  explicitly  incorporating  clutter-path  signals  into  the  RS-GLRT  detec¬ 
tor  formulation  (as  in,  e.g.,  [19])  will  result  in  a  detector  that  implicitly  mitigates 
clutter-path  interference.  This  will  extend  the  applicability  of  the  RS-GLRT  to  sce¬ 
narios  with  challenging  interference  environments. 

Second,  the  unified  framework  presented  in  this  work  could  be  extended  to  in¬ 
clude  localization  in  addition  to  detection.  Maximizing  the  RS-GLRT  statistic  over 
a  hypothesized  target  position-velocity  space,  flsurv,  gives  the  MLE  of  the  unknown 
target  state,  i.e., 

(t,  t)  =  arg  max  £rs(p,  p).  (214) 

(p.p)en  surv 

The  corresponding  Cramer- Rao  Lower  Bound  (CRLB)  for  localization  in  PAIR  sensor 
networks  does  not  currently  exist.  Note  that  (214)  represents  a  single-step  approach 
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that  directly  estimates  the  MLE  of  the  unknown  target  state.  This  is  in  contrast  to 
traditional  two-step  localization  approaches  that  indirectly  estimate  the  target  state 
by  measuring  and  then  fusing  intermediate  quantities  such  as  TDOA,  FDOA,  and 
AO  A.  A  number  of  recent  papers  have  investigated  direct  single-step  localization 
performance  for  PSL  sensor  networks,  an  approach  termed  direct  position  determi¬ 
nation  [87,  5,  12,  86,  82],  It  is  expected  that  the  CRLB  for  PMR  localization  may 
be  similarly  derived.  It  is  also  expected  from  the  ambiguity  results  in  Sec.  4.5.2  that 
this  CRLB  would  be  a  function  of  both  SNR  and  DNR ,  and  would  implicitly  ex¬ 
ploit  AOA,  TDOA,  FDOA,  bistatic  range,  and  bistatic  Doppler  information,  thereby 
unifying  localization  in  AMR  and  PSL  sensor  networks  within  a  common  framework. 

Finally,  this  work  requires  experimental  verification.  Specifically,  it  would  be 
interesting  to  verify  the  ambiguity  characteristics  of  the  RS-GLRT  detector  in  the 
high-DNR,  transition,  and  low-p  regions.  A  minimal  experimental  setup  consists  of 
one  isotropic  target,  one  pulsed  transmitter,  and  multiple  geographically-distributed 
single-channel  receivers  in  an  anechoic  RF  environment.  Achieving  the  SNR  and  DNR 
required  for  each  region  requires  careful  control  of  the  transmitter  power,  target  radar 
cross  section,  and  scenario  geometry.  Use  of  pulsed  rather  than  CW  transmission  is 
convenient  because  it  potentially  enables  isolation  of  the  direct-path  and  target-path 
signals  via  time  gating,  provided  the  pulse  duration  is  sufficiently  short  relative  to  the 
scenario  geometry.  This  enables  the  use  of  single-channel  rather  than  multichannel 
receivers,  because  multiple  channels  would  not  be  required  to  spatially  isolate  the 
direct-path  and  target-path  signals.  An  anechoic  RF  environment  is  required  to  min¬ 
imize  clutter-path  interference,  which  is  not  addressed  in  the  RS-GLRT  formulation. 
Verification  of  the  ambiguity  characteristics  associated  with  each  of  the  three  perfor¬ 
mance  regions  would  confirm  the  relationship  between  AMR,  PMR,  and  PSL  sensor 
networks  within  the  unified  framework  of  Chapter  V. 


Ill 


Appendix  A.  Far-Field  Differential  Range  Approximation 

Let  the  nth  element  of  an  array  antenna  be  at  location  rn  =  r  +  Sn,  where  r  is  the 
location  of  the  array  reference  element,  and  dn  represents  the  offset  of  the  nth  element 
with  respect  to  the  reference  element.  Let  the  location  x  be  in  the  far-held  of  the 
array,  where  the  far-held  is  defined  as  fig  =  {x:  ||x  —  r||  3>  2  ||<5n||2  /  A}.  Then,  for 
x  e  Off,  the  range  between  x  and  rn,  Rn  =  [|x  —  rn|| ,  can  be  approximated  as  the 
sum  of  the  range  to  the  array  reference  element,  i?0  =  ||x  —  r||,  and  a  differential  term 
that  accounts  for  the  offset  of  the  nth  element  with  respect  to  the  reference  element. 
To  show  this,  rn  can  be  expanded  and  simplified  as  follows, 


Rn  =  1 1  x  —  r. 


(x  -  r)  - 


I x  -  r || 2  -  2(x  -  r)T<5„,  +  ||<>,n 


=  llx  -  r ||  Wl  - 


2(x  -  r)T8n  +  ||<h, 
||x-r||2 

\Ts 


=  llx-rlUl-^-^  +  'W1  /2+H.O.T. 


x  -  r 


lx  -  r ||  1  - 


(x  -  r)TSn 


x  -  r 


=  Rq  —  k(x)  •  Sr 


(215) 

(216) 

(217) 

(218) 

(219) 

(220) 
(221) 


where  k(x)  is  the  unit  vector  pointing  from  r  to  x.  Equation  (219)  results  from  the 
following  Taylor  series  expansion  about  x  =  0,  which  converges  for  |x|  <  1, 


Vl  +  x  = 


V  (~l)w(2n)! 

(1  —  2n)(n!)2(4n) 


8  16 


(222) 


Equation  (220)  then  follows  by  ignoring  the  higher  order  terms  due  to  the  large  mag¬ 
nitude  of  1 1 x  —  r 1 1 ,  and  observing  that  ||x  —  r||  >>  ||  A  ||2  due  to  the  far-held  condition. 
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Appendix  B.  Reference  and  Surveillance  Channel  Formation 


Given  an  arbitrary  unitary  matrix  Ay  e  C N^LxN^L ^  pie  sca]ar  ||s*i  _  Mfu*||2  in 
(135)  may  be  expressed  as 

||s^'  -  M)V:||2 

=  ||A<J'(s<J'-MiJu<)||2  (223) 

=  1 1 Aijsij 1 1 2  -  2Re { s^Tvi^u* }  +  1 1  A'AI'/u' 1 1 2  (224) 

We  consider  unitary  matrices  of  the  form  A*-7  =  (B,J  ®  IL)ff ,  where  BiJ  G  CNeXNe  is 
a  unitary  beamforming  matrix  that  consists  of  Ne  orthonormal  column  vectors, 


B^  =  [b^b?,biv..,b 


(225) 


and  bf7  and  by  denote  surveillance  and  reference  beamformers,  respectively,  dehned 
according  to  the  approach  depicted  in  Fig.  25.  Specifically,  by  steers  in  the  ith  trans¬ 
mitter  direction, 

(226) 


by  = 

r  imoi 


and  bf7  steers  in  the  direction  of  the  position  under  test,  p,  while  placing  a  null  in 
the  transmitter  direction, 

po  aij 

btj  =  ft  p ....  (227) 


ip v  | 
I  %±dp  I 


where  =  Ijve  —  bybyH  is  the  projection  matrix  into  the  orthogonal  complement  of 
b^.  The  remaining  columns  (b^  :  k  —  3 . . .  Ae}  are  dehned  so  that  (B’J)/,By  =  1^. 

A  geometric  interpretation  of  BtJ  is  given  in  Fig.  41  [72],  For  the  ?jth  bistatic  pair, 
the  direct-path  and  hypothesized  target-path  steering  vectors,  a(]  and  a*7,  are  linearly 
independent  but  not  necessarily  orthogonal,  and  they  span  a  rank- 2  subspace  of  CNe, 
denoted  ( [aj  a*5"]).  The  reference  and  surveillance  beamformers,  by  and  bf7,  are  an 
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<[*j?  ay'])-1 


Figure  41.  Decomposition  of  Cff'  into  Subspaces  Spanned  by  the  Columns  of  B,?  ;  the 
Rank-2  Subspace  ( !a)j  a'/ ] }  is  Depicted  by  the  Gray  Plane 

orthonormal  basis  for  a^]).  Similarly,  {b^  :  k  =  3 . . .  Ne }  is  an  orthonormal 
basis  for  the  orthogonal  complement  of  ( [aj  a*J]),  denoted  ( [a'J  a^])^.  Note  that 


B  B  " 


Ne 

b%  f  +  bWH  +  V  Whf  =  I 


=p' 


Ne 


(228) 


where  Py[s,  Py7,  and  P^±  are  orthogonal  projections  into  (b*-7),  (b*-7),  and  ([a^  a^])-1, 
respectively. 

The  terms  of  (224)  are  considered  in  turn.  First,  ||  A'-'s*-7 1|2  can  be  expanded 


Aijsij 


(B*'7  ®  Il)H  s®-7"  1 1 2 

(229) 

Ne 

s?ir+ n^'ir +EW  ii2 

(230) 

k= 3 

=Eij  , 

(rs)-*- 
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where  s*J  and  denote  surveillance  and  reference  channel  signals,  respectively, 


Nt 


n=  1 


(231) 


and  Wfr  is  given  by 


Ne 


11 

w/  = 


wk 


n=  1 


(232) 


The  notation  [x]„  denotes  the  nth  element  of  the  vector  x.  Note  that  ±  represents 
the  total  energy  within  the  projection  of  s into  ([a^J  a^])-1. 

Next,  the  value  s^Mj^ib  in  the  second  term  of  (224)  can  be  expressed  using 
(228)  as 

=  s'-'"  ((!’;)'  +  p;is  +  l%.)  X  I,)M',V  (233) 


The  term  in  (233)  involving  can  be  expanded  using  (121),  (226),  (228),  and  the 
mixed-product  property  of  the  Kronecker  product  as 


S«"(P,«  ®  IJMjV  =  s«"(P«®  Ii)  (' tf(a«®2>« )  +T?(a«®X>«))ui  (234) 

=  s«*  (t?(P?  a«)  ®  V‘j  +  7«(P«  a«)  ®  C«)u‘  (235) 

=  (j'j a-d  ®  +  %’P'ip  aj  ®  u‘  (236) 


where  /3^p  quantihes  the  mismatch  between  a)/  and  a^ , 


1 1  O  *•?  1 1  2 


(237) 


The  second  term  in  (236),  which  represents  leakage  of  the  target-path  signal  into  the 
reference  channel,  can  be  ignored  due  to  the  large  power  ratio  between  the  direct-path 
and  target-path  signals  typical  in  PMR,  i.e.,  \YJ\2  3>  and  noting  that  \/3iJp\  <  1 
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by  the  Cauchy- Schwarz  inequality  for  a'j  ^  a*U  Continuing  with  this  approximation, 

s'-'":  p;j  X  I/.  iM'.bi'  «  7Jsij'H(aJ'  X  •p'Mu'  (238) 

/  iVe  \  H 

=  (239) 

'  n=l  ' 

=  Tjv/At(S?)''ui  (240) 


where  v^e  —  ||a^||,  and  s*J  is  the  reference  signal  in  (231)  after  delay-Dopplcr 
compensation  that  removes  the  delay  and  Doppler  of  the  direct-path  channel, 

Ne 

%  =  vf  ]T[a =  cf  s«  (241) 

n—  1 

Similarly,  the  term  in  (233)  involving  P*(_  can  be  expanded  using  (121)  and  (228), 
and  simplified,  giving 


5hH(ph 


H 


U 


(242) 


where  Ne(l  —  \(3ljp\2)  =  ||P)'la))||,  and  §*/  is  the  surveillance  signal  s\j  in  (231)  after 
dclay-Doppler  compensation  that  removes  the  delay  and  Doppler  of  the  target-path 
channel, 

Ne 

%  =  Vf  £  [by]  *n  4  =  vf  s?  (243) 

n= 1 

Unlike  (240),  (242)  is  exact  rather  than  approximate. 

Finally,  the  term  in  (233)  involving  P^x  equals  zero  because  a ^  and  a))  are  in 
the  null  space  of  P^x , 


J3HrD'li 


(rs)J 


Il)M>1  =  8«Jr(7J'(P^aJ')®  Z>?'+7?'(P, 


(rs)J 


U 


(244) 


=  o 


=  o 
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Substituting  (240),  (242),  and  (244)  into  (233)  gives 


s'''"'  m  (245) 

where  ^  and  are  scale  factors  that  account  for  the  composite  scaling  of  the 
reference  and  surveillance  channels  resulting  from  channel  effects  and  beamforming, 
defined  as 


K?  =  (246) 

-  Kf)  (247) 

11  ‘112 

Finally,  the  term  || A'-'M'/u' ||  in  (224)  can  be  expanded  using  (121)  and  (225), 
and  simplified,  to  give 


||A*JM;V||2  «  L(\$\2  +  |^|2)  (248) 

The  approximation  in  (248)  amounts  to  ignoring  the  target-path  signal  leakage  in  the 
reference  channel,  discussed  earlier  in  relation  to  (236). 

Substituting  (230),  (245),  and  (248)  into  (224),  and  noting  that  ||s*J'||2  =  ||s^||2, 
Ill'll2  =  ||s^'||2,  and  ||iT||2  =  L,  the  scalar  ||su  —  MlJ’u*||2  can  be  expressed  as 

1 1 sij  -  ||2«||s^  ^V||2+|| sj*  -  /^u* 1 1 2  +  E^s)±  (249) 
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